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ABSTRACT 


This  thesis  describes  the  analysis  and  development  of  an  acoustic  vorticity  meter  to 
measure  shear  in  ocean-boundary  layers  over  smaller  measurement  volumes  than 
previously  possible.  A  nonintrusive  measurement  of  vorticity  would  filter  out  irrotational 
motion  such  as  surface  waves  and  currents  that  can  swamp  small  scale  measurements  of 
shear.  The  thesis  describes  the  desired  geophysical  measurements  and  translates  this 
oceanographic  context  into  design  goals. 

The  instrument  was  designed,  built,  tested,  and  deployed.  It  measures  three-axis 
vorticity  at  0.83  and  2.45  meters  below  the  ocean  surface  with  measurement  volumes  of 
0.45  meters  on  a  side.  The  instrument  forms  a  buoy  that  is  inertially  instrumented  to 
calculate  and  remove  buoy  motion  from  the  measurements.  The  instrument  uses  a 
complementary  filter  algorithm  to  estimate  attitude  and  motion  from  low-power, 
inexpensive,  strapdown  rate  gyros,  accelerometers,  and  fluxgate  magnetometers.  The 
instrument  performance  has  been  measured  to  have  a  vorticity  bias  of  not  more  than  1  x 
10'2  per  second  in  a  mean  flow  of  0.7  meters  per  second,  a  bias  of  not  more  than  1  x  10'2 
per  second  in  the  down- wave  and  vertical  directions  in  typical  ocean  waves,  and  a  30 
decibel  spectral  rejection  of  surface  wave  velocity. 

Two  instrument  deployments  are  described  to  show  the  potential  of  the  system. 
The  instrument  has  measured  shear  in  the  upper-ocean-boundary  layer,  and  these 
measurements  are  compared  to  concurrently  measured  wind  stress  and  stratification.  The 
instrument  was  also  deployed,  tethered  in  the  thermocline,  in  an  area  of  high  internal  wave 
activity.  Richardson-number  time  series  were  measured  and  compared  favorably  to 
concurrently  measured  Richardson  numbers  made  over  a  larger  spatial  scale. 


Thesis  Supervisor:  Albert  J.  Williams  III. 
Title:  Senior  Scientist 
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CHAPTER  1.  INTRODUCTION 


This  thesis  documents  the  development  and  testing  of  an  instrument  system  to 
measure  ocean  vorticity  and  shear  in  the  upper,  internal,  and  bottom-boundary  layers.  The 
purpose  of  this  project  has  been  to  provide  ocean  scientists  with  a  tool  capable  of 
measuring  ocean  shear  with  finer  resolution  and  or  closer  to  the  boundaries.  The 
deployment  data  are  presented  to  demonstrate  the  instrument's  potential  rather  than  to 
contribute  directly  to  understanding  marine  boundary  layers. 

This  chapter  describes  the  geophysical  environment  to  be  measured,  why  vorticity 
is  used  to  measure  shear,  how  shear  can  be  used  to  measure  vertical  diffiisivity,  and 
reviews  the  tradeoff  between  resolution  and  accuracy  in  a  shear  or  vorticity  measurement. 
The  upper-boundary  layer,  internal-boundary  layer,  and  bottom-boundary  layer  are 
reviewed  to  motivate  their  measurement  and  to  define  the  instrument  design  goals. 

Chapter  two  describes  the  mechanical  design,  electrical  design,  the  inertial  measurement 
unit  and  its  inertial  processing,  and  the  signal  processing  done  on  the  raw  data.  Chapter 
three  models  and  measures  the  sensor  performance.  Chapter  four  describes  two 
deployments  of  the  instrument:  one  in  the  upper  boundary  layer  to  measure  shear,  and  one 
in  the  thermocline  to  measure  gradient  Richardson  number.  Chapter  five  concludes  the 
thesis  and  discusses  future  work  and  improvements  to  the  system. 

l.A.  OCEAN  BOUNDARY  LAYERS 

Boundary  layers  mediate  the  turbulent  fluxes  of  heat,  momentum  and  chemical 
species.  Turbulent  transport  dominates  molecular  diffusion  everywhere  except  in  diffusive 
sublayers  on  the  order  of  one  millimeter  thick.  The  water  side  of  the  air-sea  interface 
controls  the  air-sea  transfer  of  most  gasses  (Kitaigorodskii  and  Donelan,  1984).  Gas 
transfer  resistance  is  a  function  of  turbulence  in  the  boundary  layer.  Many  quantities  in 
ocean  boundary  layers  such  as  Reynolds  stress,  chemical  flux,  or  shear  are  difficult  to 
measure  in  the  presence  of  gravity-wave  velocities  that  swamp  turbulent  velocities. 
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Typical  boundary-layer  flows  have  wave  velocities  of  order  0.5  m/s  and  shear  of  order  2  * 
1  O'2  per  second.  Over  large  measurement  separations  ie.  ten  meters,  this  shear  is  readily 
measured.  Yet  over  0.5  meters,  if  the  shear  is  to  be  measured  with  current  meters,  the 
meter's  accuracy  would  need  to  be  better  than  one  percent.  In  the  presence  of  waves  this 
accuracy  is  not  readily  achievable.  This  difficulty  results  in  a  shortage  of  near-surface 
measurements  of  shear  and  poorly-calibrated-closure  relations  for  ocean  and  climate 
models.  Vorticity  can  be  measured  as  a  surrogate  for  shear  because  time-average,  area- 
average,  horizontal  vorticity  equals  time-average,  area-average,  vertical  shear  in  ocean¬ 
boundary  layers.  A  nonintrusive  measurement  of  vorticity  would  measure  wind-driven 
shear  and  automatically  remove  irrotational  surface-gravity-wave  and  current  velocities, 
allowing  measurement  of  shear. 

This  thesis  will  emphasize  measurements  for  studies  of  vertical  diffusivity.  For 
readers  not  too  familiar  with  geophysical  flows,  vertical  (diapycnal)  diffusivity  is  orders  of 
magnitude  smaller  than  horizontal  (isopycnal)  diffusivity.  In  the  open  ocean,  far  away 
from  boundaries,  a  typical  horizontal  diffusivity  is  3  m2  /s  while  the  vertical  diffusivity  is 
just  1*10'5  m2  /s  (Ledwell,  Watson,  and  Law,  1993).  Density  structure  clearly  plays  a 
large  role  in  ocean  diffusivity. 

UPPER-BOUNDARY  LAYER 

The  measurement  of  upper-boundary-layer  shear  is  the  focus  of  this  instrument 
development.  The  ocean-upper-boundary  layer  is  often  compared  to  the  well  studied  and 
relatively  well  understood  unstratified  turbulent  flow  over  a  rigid  wall.  Factors  that  can 
complicate  the  ocean  surface  are  waves,  wave  breaking,  stratification,  rotation,  and 
organized  motions  such  as  Langmuir  cells. 

Shear  in  a  constant  stress  layer  in  an  unstratified  turbulent  flow  next  to  a  rigid  wall, 

zu 

is  given  by  equation  (1-1)  for  distances  greater  than  — >50  from  the  wall,  of  order  5  mm 

V  2 

in  the  ocean  (Monin_and  Yaglom,  1987).  In  this  equation,  —  is  Eulerian  shear,  w.is  the 
.  .  .  t  dz 

friction  velocity  —  the  square  root  of  the  shear  stress  divided  by  density, 

N  P 
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du  =  ^ 
dz  K 2 


(1-1) 


k  is  von  Karmon's  constant  usually  assumed  to  be  0.4,  and  z  is  the  distance  from  the  wall. 
This  shear  will  give  a  logarithmic  velocity  profile,  equation  (1-2),  and  the  dissipation  will 

u 

u(z)  =  —  ln(z)  +  const.  (1-2) 

K 

3 

be  —  (Tennekes  and  Lumley,  1989).  The  assumptions  for  this  log  layer  are  constant 

.  KZ 

shear  stress,  constant  density  and  a  Reynolds  number  high  enough  that  viscous  forces  are 
negligible  compared  to  turbulent  Reynolds  stress.  By  comparing  measured  ocean  shear  to 
measured  windstress,  the  effects  on  vertical  mixing  effectiveness  of  stratification,  wave 
breaking,  surfactants,  and  Langmuir  cells  can  be  measured. 


Stratification  can  inhibit  turbulence.  To  show  how  this  happens,  the  turbulent 
kinetic  energy  budget  equation  will  be  reviewed,  equation  (1-3)  (Stull,  1988).  In  this 
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equation,  term  a  is  the  material  derivative  of  turbulent  kinetic  energy  q  2  =u  2  +v  2  +w  2 , 
term  b  is  creation  of  turbulent  energy  by  shear,  term  c  is  pressure  diffusion,  term  d  is 
energy  diffusion,  term  e  is  viscous  dissipation,  and  term  /is  destruction  by  buoyancy.  In  a 
steady-state,  horizontally-uniform,  boundary  layer,  term  a  is  assumed  zero  and  the 
transport  terms  c  and  d  are  assumed  very  small  leaving  equation  (1-4).  The  energy 
destruction-by-buoyancy  term  takes  energy  away  from  turbulence. 


15 


(1-4) 


/  /  R  /  / 

-M  JV  -  «  6  +  — O  JV 

dz  p 

The  wind  driven  current  as  a  function  of  depth  is  a  function  of  the  vertical-eddy 
diffiisivity,  which  in  turn  is  strongly  influenced  by  stratification  (Price,  Weller,  and 
Schudlich,  1987).  If  a  constant  wind  stress  and  a  constant  vertical-eddy  diffusivity  with 
depth  is  assumed,  the  momentum  balance,  including  Coriolis  force,  generates  an  Ekman 
spiral.  If  eddy  diffiisivity  is  assumed  to  be  proportional  to  depth,  a  different  result  is 
obtained  which  is  closer  to  most  ocean  measurements  of  current  shear  (Madsen,  1977). 
The  solution  is  a  function  of  the  vertical  eddy  diffiisivity,  which  in  turn  is  a  function  of 
stratification  and  any  other  modifications  to  the  turbulence  such  as  wave  breaking  energy 
addition,  or  Langmuir  cells.  The  instrument  developed  in  this  thesis  measures  shears  close 
to  the  ocean  surface,  in  the  upper  three  meters,  where  there  are  few  shear  measurements. 

FREE  SURFACE  CONSIDERA  TIONS 

Free-surface  considerations  that  can  be  important  to  understanding  the  ocean- 
upper-boundary  layer  include  Langmuir  cells,  a  reduced-shear  layer  (as  compared  to  the 
turbulent  wall  layer  model),  gas  transfer,  wave  breaking  and  their  resultant  bubbles  and 
droplets,  and  surfactants. 

Langmuir  cells  are  three-dimensional  drift  currents  that  form  counter-rotating 
helical  vortices  parallel  to  the  wind.  Fig  1-1  (Langmuir,  1938).  Langmuir  cells  are  often 
identified  by  windrows  of  debris  floating  parallel  to  the  wind,  and  arranged  by 
convergence  zones  of  the  vortices.  When  these  cells  exist,  they  dominate  vertical  mixing 
over  their  extent  (Gordon,  1970).  While  cell-averaged  shears  are  of  order  1  x  10'2/s,  shear 
is  concentrated  on  the  edges  of  the  down  welling  jets  and  can  be  significant  compared  to 
average  current  shear.  Langmuir  cells  are  important  in  understanding  the  upper-boundary- 
layer  dynamics  and  complicate  measurements  in  this  environment  by  making  time-average 
shear  and  vorticity  horizontally  nonuniform. 


16 


Figure  1-1.  Artists  conception  of  Langmuir  cells 

The  transfer  of  momentum  from  wind  to  waves  and  current  is  not  well  understood 
and  is  an  area  of  active  research.  While  the  wind's  momentum  is  transferred  initially 
largely  to  waves,  and  most  of  this  to  short  wavelets,  most  of  this  momentum  is 
transferred  through  breaking  to  Eulerian  shear  and  currents.  Less  than  six  percent  of  the 
transferred  momentum  eventually  is  radiated  away  as  waves  (Mitsuyasu,  1985).  Many 
researchers  have  found  a  layer  of  reduced  shear  and  enhanced  dissipation  (some  by  a 
factor  of  100)  as  compared  to  a  wall  layer  (Csanady,  1983  and  1984,  Cheung  and  Street, 
1988  a  and  b,  Agrawal  et  al,  1992).  Santala  (1991)  reports  a  zone  of  no  shear  in  the 
direction  of  windstress  to  a  depth  of  ~  1.2  *  1 05  and  significant  shear  at  right 
angles  to  the  windstress  .  One  of  the  gcfals  in  developing  the  vorticity  meter  is  to  be  able 
to  make  better  open-ocean  measurements  of  shear  close  enough  to  the  surface  to  study 
this  reduced  shear  phenomenon. 
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Research  on  gas  transfer  between  the  ocean  and  atmosphere  has  received  increased 
interest  with  the  concern  over  atmospheric  carbon  dioxide  buildup.  The  transfer  of  low- 
solubility  gases  is  controlled  by  the  water-side  diffusive  sublayer,  which  is  of  order  one 
millimeter  thick  (Kitaigorodskii  and  Donelan,  1984).  Straining  and  renewal  of  the 
diffusive  sublayer  by  upper-boundary-layer  turbulence  is  a  major  contributor  to  gas 
transfer  (Brumley  and  Jirka,  1988).  Bubbles  from  wave  breaking  can  significantly  increase 
gas  transfer  when  the  wind  exceeds  some  velocity  (Broecker  and  Siems,  1984).  In 
addition  to  their  increased  effective-surface  area  for  diffusion,  bubbles  add  turbulence  and 
strain  the  diffusive  sublayer. 

The  presence  of  surfactants  affect  the  surface  boundary  condition,  dampen 
capillary  waves,  and  reduce  gas  transfer  (Hunt,  1984).  While  surfactants  reduce  the 
surface  drag  coefficient,  Wu  (1983)  reports  that  surfactants  can  still  increase  surface  drift 
currents. 

SENSOR-  WA  VE-CORRELA  TION  BIA  S 

Sensor-wave-correlation  bias  of  velocity  and  velocity-derived-shear  measurements 
is  caused  by  a  correlation  between  sensor  motion  and  the  wave-field-velocity  gradient 
(Pollard,  1973).  It  is  different  from  Stokes  drift  and  Stokes-drift  shear  because  sensors 
and  buoys  do  not  exactly  follow  water  motion.  In  open-ocean  fetch  conditions,  this  bias  is 
greater  than  wind-driven  Eulerian  shear  (Wu,  1975).  Santala  (1991)  and  Santala  and 
Terray  (1992)  derived  an  algorithm  for  removing  this  bias  from  measurements  taken  from 
an  inertially  instrumented  buoy,  but  the  accuracy  of  this  bias  removal  is  limited  by  the 
accuracy  of  the  measurements  of  the  directional-wave  spectrum.  A  nonintrusive 
measurement  of  vorticity,  on  the  other  hand,  is  not  effected  by  any  motion-correlation 
bias. 


REASONS  TO  MEASURE  EULERIAN  SHEAR 

In  addition  to  measuring  Eulerian  shear  to  understand  the  structure  of  surface 
boundary  layers,  one  reason  for  measuring  Eulerian  shear  is  to  be  able  to  determine  the 
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effective  turbulent  diffusivity.  Comparing  wind  stress  to  Eulerian  shear  results  in  a 
measurement  of  the  vertical-eddy  diffusivity,  and  by  the  Reynolds  analogy  for  turbulent 
diffusion,  is  also  a  measurement  of  the  effectiveness  of  turbulent  transfer  of  heat  and 
solutes  (Rohsenow  and  Choi,  1961).  This  measurement  of  turbulence,  although  indirect, 
is  better  conditioned  than  direct  measurements  of  Reynolds  stress  near  the  ocean  surface 
because  of  the  very  large  quadrature  components  from  wave  velocities  that  do  not 
contribute  to  the  vertical  flux  of  momentum,  heat  or  solutes.  Any  wave  reflections 
(standing  waves)  or  phase  lag  between  velocity  measurement  axes,  dooms  a  measurement 
of  Reynolds  flux  in  the  wave  field. 

MEASUREMENT  TYPES 

Different  methods  for  measuring  drift  currents  and  shear  include  Lagrangian 
drifters,  Eulerian  measurements,  and  surface-referenced  buoys.  For  measuring  mass 
transport  velocities,  Lagrangian  drifters  seem  like  an  obvious  inexpensive  option.  The 
problems  with  drifters  include  uncalibrated  drift  relative  to  their  target  depth  (Geyer, 

1989),  and  their  tendency  to  get  stuck  in  convergence  zones  of  three-dimensional  flow 
structures  such  as  Langmuir  cells,  which  can  have  unrepresentative  drift  velocities.  If 
measurement  of  Eulerian  shear  is  desired,  the  correlation  bias  of  sensor-wave  motion  of 
the  drifters  would  have  to  be  compensated  for,  which  would  require  inertial 
instrumentation  and  measuring  relative  velocity,  significantly  increasing  the  drifter's  cost 
and  complexity.  An  Eulerian  (fixed  in  space)  measurement  of  velocity  or  shear  requires  a 
tower  which  is  expensive,  has  a  large  flow  obstruction,  and  cannot  be  used  in  deep  water. 
Three-point  moorings  and  taut  moorings  are  not  stiff  enough  to  make  Eulerian-velocity 
measurements  free  of  motion-correlation  bias.  Surface-referenced  buoys  are  often  used  to 
measure  ocean  currents  but  buoy  motion  must  be  known  in  order  to  remove  wave  bias 
(Santala,  1991  and  Santala  and  Terray,  1992).  Buoys  ride  with  swell,  giving  a  non- 
Eulerian  reference  frame.  This  swell-based  reference  frame  can  actually  assist  near-surface 
measurements,  allowing  velocity  measurements  within  a  waveheight  of  the  surface 
(Cheung  and  Street,  1988).  Measuring  vorticity  from  a  surface  buoy  avoids  motion- 
correlation  bias,  but  the  buoy  must  be  inertially  instrumented  to  remove  buoy  motion  from 
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the  buoy's  measurements.  If  the  buoy  motion  is  not  known  and  not  compensated  for, 
motions  such  as  coning  motions  could  result  in  time  average  sensor  measurements  that  are 
caused  by  the  buoy  motion  and  do  not  exist  in  the  fluid.  Coning  motions  will  be  explained 
in  the  inertial  processing  section. 

INTERNAL  BOUNDARY  LAYER 

Density-stratified  layers  in  the  ocean  inhibit  turbulent-vertical  mixing,  support 
internal  waves,  and  can  be  treated  as  a  boundary  layer  (Salmon,  1990).  The  gradient 
Richardson  number  equation  (1-5)  is  a  ratio  of  the  relative  strengths  of  stratification. 


N2  P„  dz 
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which  inhibits  mixing,  and  shear,  which  encourages  mixing  (Turner,  1973).  In  this 

equation,  N  is  the  Brunt-Vaisala  frequency,  is  the  vertical  shear,  g  is  gravity,  p0  is 
dpr  .  .  dz 

average  density,  and  — -  is  the  density  gradient.  A  gradient  Richardson  number  less  than 
dz 

one  quarter  is  a  necessary,  and  in  practice  usually  sufficient,  condition  for  turbulent 
vertical  mixing  to  occur.  Internal  waves  have  vorticity  in  a  density-stratified  layer. 
Intemal-wave-shear  instability  is  thought  to  be  the  major  source  of  vertical  mixing  in 
density-stratified  layers  of  the  ocean  away  from  boundaries;  therefore  accumulating  more 
gradient-Richardson-number  statistics  in  the  ocean  over  different  spatial  scales  will  help  in 
further  understanding  vertical  mixing  in  the  ocean  (Gargett  et  al,  1981,  and  Gargett  and 
Holloway,  1984). 


BOTTOM-BOUNDARY  LAYER 

Understanding  turbulent  transport  and  stress  in  the  bottom-boundary  layer  is 
necessary  to  understand  processes  of  sediment  transport.  Measurement  of  bottom- 
boundary-layer  turbulence  near  shore  is  complicated  by  surface  wave  swell  that  reaches 
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the  bottom,  is  often  much  more  energetic  than  the  turbulence,  and  occupies  the  same 
frequency  band  (Grant  and  Madsen,  1986).  The  ability  to  measure  turbulence,  free  of 
surface  swell,  should  assist  ocean  scientists  in  studying  nonlinear  wave-current  interaction. 
By  measuring  vorticity,  the  surface  swell  is  filtered  out  leaving  the  turbulence  to  be  more 
readily  measured. 

I.B.  PREVIOUS  WORK 

Other  researchers  have  made  single-axis  vorticity  measurements  in  laboratories  and 
in  the  ocean.  Rossby  (1975)  proposed  measuring  ocean  vorticity  and  made  a  laboratory 
demonstration  of  single  axis  vorticity  measurement  by  measuring  the  circulation  around  a 
closed  triangle  by  measuring  the  difference  in  acoustic  travel  time.  The  lab  demonstration 
used  acoustic  mirrors  to  form  the  circulation  triangle.  He  proposed  measuring  circulation 
around  a  large  triangle  of  from  3  kilometers  on  a  side  to  ocean-basin  size.  Tsinober,  Kit, 
and  Teitel  (1986)  measured  single-axis  vorticity  electromagnetically  in  a  lab  over  very 
small  scale.  They  made  a  seven  electrode  probe  to  make  a  central  difference 
approximation  of  the  divergence  of  the  electric  potential  which  is  proportional  to  fluid 
vorticity.  Their  probe  was  only  2  millimeters  across.  Muller,  Lien,  and  Williams  (1988) 
estimated  relative  vorticity  in  the  ocean  from  current  meters  that  were  trimoored  in  the 
ocean  thermocline.  The  current  meters  formed  horizontal  triangles  from  8.5  to  1600 
meters  on  a  side  and  the  measurements  were  not  compensated  for  mooring  motion. 
Menemenlis  and  Farmer  (1992)  measured  vertical  vorticity  8  and  20  meters  beneath  the 
arctic  ice  sheet  acoustically.  They  measured  circulation  around  a  triangle  200  meters  on  a 
side.  Tom  Sanford  (APLUW,  personal  communication)  has  measured  small-scale,  single¬ 
axis  vorticity  in  ocean  boundary  layers  by  measuring  the  divergence  of  the  electric  field. 
The  instrument  whose  development  is  described  in  this  thesis,  has  made  the  first  three-axis 
vorticity  measurements  in  ocean  boundary  layers. 

l.C.  INSTRUMENT  DESIGN  GOALS 

The  primary  goal  in  developing  this  instrument  is  to  have  the  ability  to  measure 
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Eulerian  shear  in  the  ocean-upper-boundary-layer,  over  finer  resolution  and  closer  to  the 
surface,  than  previously  practical.  The  system  should  be  deployable  by  a  moderately  sized 
oceanographic  vessel  such  as  the  R.  V.  Asterias  (a  15  meter  workboat),  and  the 
measurement  volumes  should  be  scaled  and  at  a  depth  to  measure  the  reduced  shear  layer. 
Below  a  depth  of  about  5  meters,  moored  current  meters  are  adequate  for  measuring 
shear.  This  instrument  measures  shear  in  the  upper  5  meters  of  the  ocean  and  can  measure 
shear  in  the  upper  meter  of  the  ocean  in  order  to  measure  the  reduced  shear  layer.  A 
secondary  goal  for  the  system  is  that  it  be  able  to  measure  internal-wave  shear  and 
stratification  in  the  internal-boundary  layer  (thermocline)  over  small  scales  in  the  presence 
of  surface  swell.  An  ancillary  use  for  these  sensors  was  to  measure  turbulence  in  the 
bottom-boundary  layer  in  the  presence  of  surface  swell.  This  later  application  of  the 
vorticity  meter  is  not  covered  in  detail  in  this  thesis,  but  data  from  this  application  is 
presented  to  show  swell-spectral  rejection  in  a  coastal  deployment. 
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CHAPTER  2.  INSTRUMENT  DESIGN 


Chapter  two  describes  several  aspects  of  the  mechanical  design,  electrical  design, 
and  signal  processing.  The  signal  processing  section  reviews  the  calculation  of  instrument 
attitude  and  motion  in  inertial  space,  develops  an  algorithm  that  can  use  inexpensive  low- 
power  sensors  to  measure  instrument  motion,  describes  the  data  processing  done  on  all 
the  raw  signals,  and  simulates  the  inertial  processing.  The  inertial  simulations  show  that 
when  using  this  inertial  algorithm,  errors  in  measuring  vorticity  and  shear  resulting  from 
imperfect  knowledge  of  buoy  motion,  are  much  smaller  than  errors  resulting  from  flow 
disturbance.  Readers  not  interested  in  inertial  processing  can  skip  the  signal  processing 
subsections,  except  the  data  processing  subsection,  and  still  follow  the  rest  of  the  thesis. 

2.A.  MECHANICAL  DESIGN  OVERVIEW 

In  this  instrument,  circulation  around  a  closed  square  path  is  measured,  which  via 
Stokes  theorem,  is  the  area-integrated  vorticity  over  the  surrounded  area,  equation  (2-1). 


jjAd-ndA  =  jv-ds  =  T 
65  =  Vxv 


(2-1) 


In  this  equation  vis  velocity,  (5  is  vorticity  and  T  is  circulation.  Vorticity  is  the  curl  of  the 
velocity  field.  The  z-direction  vertical  shear  of  a  horizontal  x-directed  velocity 
corresponds  to  a  vorticity  in  the  y-direction.  A  square  was  chosen  for  the  circulation  path 
to  minimize  disturbance  to  circulation  in  a  constant  mean  flow.  A  triangle  has  fewer 
paths,  but  as  Fig.  2-1  shows,  a  mean  current  could  produce  a  wake  on  one  side  that  would 
generate  a  significant  measured  circulation  where  one  did  not  exist  in  the  undisturbed 
flow.  The  wakes  resulting  from  a  symmetrical  square  of  transducers  largely  cancel  out, 
minimizing  the  vorticity  error  due  to  wakes. 


The  water  velocity  in  each  path  around  the  square  is  measured  acoustically. 
Differential  acoustic  travel  times  on  each  path  are  measured  by  modified  Benthic  Acoustic 
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Fig.  2-1  Circulation  error  from  transducer  wakes  largely  cancel  out  with  a  square 
geometry  but  not  with  a  triangular  geometry. 


Stress  Sensor  (BASS)  electronics  (Williams  et  al,  1987).  The  water  velocity  v  parallel  to 
each  acoustic  path  is  given  by  equation  (2-2)  where  c  is  the  speed  of  sound.  At  is  the 
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differential  travel  time,  and  L  is  a  single,  acoustic-path  length  making  up  one  side  of  a 
square.  BASS  electronics  were  chosen  to  measure  velocity  because  of  their  speed, 
accuracy,  and  low  noise;  single-pulse  time  noise  is  forty  picoseconds.  The  option  of  using 
acoustic  mirrors  and  fewer  transducers  was  dismissed  due  to  problems  in  distinguishing 
the  signal  from  reflections  off  the  structure. 


The  velocity  of  each  path  around  the  square  is  added  to  give  the  circulation 
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divided  by  path  length  (Fig.  2-2).  Circulation  divided  by  path  length  is  plotted  instead  of 
true  circulation  or  vorticity,  because  it  has  the  same  units  as  velocity  so  that  cancellation 


Velocities  Around  the  Square 


Fig.  2-2.  Results  from  a  typical  tow  tank  run  showing  measured  velocity  from  each 
acoustic  path  around  the  square,  circulation  divided  by  path  length,  and  vorticity. 
This  sensor  was  at  an  angle  to  the  flow. 


of  antiparallel  velocities  can  be  evaluated.  The  prototype  lab  vorticity  meter,  when 
oriented  with  two  of  the  acoustic  paths  parallel  to  the  flow,  (the  worst  orientation  for 
velocity  measurement)  measured  ten  percent  less  velocity  than  the  undisturbed  flow. 
Parallel  paths,  however,  have  the  same  decrement  to  within  one  percent,  resulting  in  small 
vorticity  error. 

A  three-axis  vorticity  sensor  was  designed  to  minimize  wake-related  errors  when 
buoy  deployed.  For  an  instrument  to  be  deployed  on  an  ocean  buoy,  all  three  axes  of 
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vorticity  have  to  be  measured  because  buoy  attitude  can  change.  Buoy  motion  can  then  be 
compensated  for,  in  signal  processing,  if  the  buoy's  motion  is  known.  A  prototype 
vorticity  meter  with  15  centimeter  acoustic  paths  was  made  and  is  shown  in  Fig.  2-3.  The 


Fig.  2-3.  Photo  of  the  prototype  three-axis  vorticity  meter  with  15  cm  acoustic 
paths. 

1 5  centimeter  path  length  was  chosen  for  convenience  in  the  lab  and  compatibility  with 
existing  electronics.  Twelve  acoustic  paths  form  the  edges  of  a  regular  octahedron  and 
form  three,  orthogonal,  circulation  path  squares.  Each  sphere  contains  four  piezoelectric 
transducers  that  form  the  corner  of  two  circulation  path  squares.  Because  a  flow  sensor 
mounted  on  a  buoy  may  measure  flow  in  any  direction,  but  cannot  be  streamlined  in  all 
directions,  the  design  philosophy  was  to  minimize  wake  asymmetries  and  errors  in 
circulation.  This  geometry  maximized  symmetry  and  minimized  measured  circulation  flow 
disturbance.  If  some  of  the  parts  are  streamlined,  such  as  the  cylinders,  they  would 
become  lifting  surfaces  when  the  flow  angle  of  attack  changes  and  create  lift-related 
circulation,  ie.,  error.  Velocity  defects  symmetric  with  respect  to  acoustic  circulation 
paths,  do  not  contribute  significantly  to  error  in  measured  circulation.  This  prototype 
sensor  was  tested  for  bias  in  uniform  flow  in  tow  tanks  and  wave  rejection  in  a  wave  tank 
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with  the  results  of  these  tests  justifying  our  building  several  larger  versions. 

The  ratio  of  circulation  signal  to  wake  induced  noise  is  expected  to  increase  with 
longer  path  length  so  the  ocean  deployed  instruments  were  scaled  up  from  the  prototype. 
Two  vorticity  sensors  with  45  centimeter  acoustic  paths,  were  built  into  an  inertially 
instrumented  buoy  to  measure  shear  Fig.  2-4.  A  vorticity  sensor  with  a  1.5  meter  acoustic 


path  length,  was  built  as  its  own  tripod  for  measuring  the  turbulent  vorticity  and  shear  in 
the  bottom-boundary  layer;  this  sensor  will  only  be  briefly  commented  on  in  the 
performance  section  where  it  is  used  to  show  spectral  rejection  of  waves  during  a  coastal 
deployment  (Williams,  Terray,  Thwaites,  and  Trowbridge,  1994). 

The  sensor  geometry  uses  a  center  structural  stalk  to  reduce  circulation 
disturbance  as  compared  to  using  a  braced  frame.  A  braced  frame  such  as  one  that  BASS 
uses,  has  a  much  better  strength  and  stiffness  to  weight  ratio,  and  can  utilize  smaller 
structural  members.  Its  members  are  then,  however,  at  a  radius  from  the  measurement 
volume  center,  and  any  asymmetry  between  structural  member  wakes  and  acoustic  paths 
would  create  large  circulation  disturbances.  Tube  bending  strength  is  proportional  to  the 
tube  diameter  cubed  and  bending  stiffness  to  the  diameter  to  the  fourth  power,  while  drag 
and  consequent  wake  size  is  linear  to  the  tube  diameter.  If  one  is  looking  for  strength  to 
drag  ratio,  a  nonbraced  larger  center  tube  is  not  at  a  disadvantage.  An  additional  benefit 
of  this  design  is  that  all  the  wires  are  strung  inside  the  tubing  and  not  exposed  to 
hydrostatic  pressure.  Hydrostatic  pressure  applied  to  cables  squeeze  the  cables  and  change 
the  cables’  capacitance,  which  can  cause  instrument  zero  drift. 

The  sensor  structure  was  kept  smooth  to  avoid  tripping  its  boundary  layer.  The 
Reynolds  number  for  a  two-centimeter  tube  or  sphere  at  one  meter  per  second  velocity  is 
about  2  *  104  which  is  below  the  laminar-to-turbulent,  boundary-layer  transition  zone. 

For  a  sphere,  this  transition  zone  Reynolds  number  ranges  from  7  *  104  to  5  *  105 
depending  on  such  variables  as  surface  roughness,  ambient  turbulence,  and  structural 
vibration  (Potter  and  Foss,  1982). 

One  drawback  of  this  octahedron  design  is  a  lack  of  redundancy.  If  any  of  the 
twelve  acoustic  paths  fails  on  a  buoy  deployment,  there  will  be  no  earth-referenced 
vorticity  without  using  the  questionable  assumption  of  zero  correlation  between  buoy 
attitude  and  water  motion.  This  lack  of  redundancy  could  only  be  solved  with  a 
substantial  increase  in  flow  disturbance. 
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VORTEX  SHEDDING 

Any  bluff  body  (at  these  Reynolds  numbers)  moving  with  respect  to  a  fluid  has 
oscillatory  forces  exerted  on  it  by  the  fluid  which  are  associated  with  the  body's  von 
Karmon  vortex  street.  Acoustic  current  meters  are  sensitive  to  instrument  strumming. 
The  largest  oscillatory  force  on  a  tube  is  perpendicular  to  both  the  tube  and  the  flow 
direction,  at  a  frequency  given  by  equation  (2-3)  (Blevins,  1977).  In  this  equation/,  is  the 

SU 

D  (2-3) 

S  «  0.2 

frequency  of  the  periodic  forcing,  S  is  the  Strouhal  number,  U  is  the  relative  velocity,  and 
D  is  the  body  diameter.  Vibration  can  have  a  strong  organizing  effect  on  the  von  Karmon 
wake,  which  increases  the  lateral  force  exciting  the  structure,  and  can  cause  lock-in  of 
lightly-damped  structures.  In  addition  to  the  large  lateral  exciting  force  just  described, 
there  is  a  smaller  oscillating  force  parallel  to  the  flow  at  twice  the  frequency  given  above 
that  can  cause  lightly  damped  structures  to  lock-in  at  this  higher  Strouhal  number 
(Crandall,  Vigander,  and  March,  1975).  The  braised  stainless  steel  structure  of  the 
vorticity  meter  has  very  little  damping,  so  that  strumming  has  to  be  considered. 

The  prototype  with  the  1 5-centimeter  path,  showed  no  problems  with  strumming. 
The  first  natural  frequency  of  the  pod  arms  was  measured  to  be  123  Hz  and  that  of  the 
center  stalk  to  be  28  Hz.  The  pod  arms  had  polyurethane  injected  into  the  annulus 
between  the  wires  and  the  tube  interiors.  While  this  did  not  effectively  add  damping,  it  did 
keep  the  cables  from  rattling  around  and  acted  as  a  pressure  block.  Smaller  structures 
naturally  have  higher  natural  frequencies  than  larger  structures. 

The  vorticity  sensors  with  the  45  centimeter  acoustic  paths,  were  designed  with 
thin  pod  arms  for  minimum  flow  disturbance.  The  pod  arms  vibrate  at  30  Hz  with  a 
damping  ratio  of  only  0.00038,  measured  with  the  logarithmic  decrement  method  over 
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1,800  cycles  (Meirovitch,  1975).  In  a  constant  flow,  the  first  strumming  occurs  at  70  cm/s 
and  is  in  the  double  frequency  parallel  mode.  As  mentioned  previously,  the  buoy  that 
these  sensors  form,  was  designed  to  ride  with  ocean  swell.  In  the  open  ocean 
deployments  the  relative  velocity  between  the  sensors  and  water  has  never  been  this  large 
so  strumming  has  not  been  a  problem. 

As  an  additional  project,  two  more  45-centimeter  path,  vorticity  sensors  were  built 
to  measure  bottom-boundary-layer  turbulence  in  high  tidal  currents.  The  design  changes 
for  these  sensors  will  be  discussed  as  well.  The  solution  we  chose  was  to  raise  the  pod 
arms'  natural  frequency  and  damping  by  adding  pipes  that  were  half  as  long  as  the  pod 
arms,  screwing  them  into  the  center  stalk,  and  filling  the  annulus  between  the  pod  arms 
and  pipe  with  polyurethane  resin.  Polyurethanes  have  high  damping  and  are  often  used  to 
absorb  vibration  energy.  The  resulting  natural  frequency  of  the  pod  arms  was  doubled  and 
had  high  damping.  These  sensors  were  tested  and  observed  in  Vineyard  Sound,  MA.  in  a 
two  knot  tidal  current  and  showed  no  signs  of  strumming.  In  designing  this  modification, 
physical  dynamic  models  were  made  and  tested  with  polyurethane  by  itself,  fiberglass  in 
epoxy,  carbon  in  epoxy,  and  fiberglass  in  polyurethane,  in  addition  to  the  double  tube  filled 
with  resin.  None  of  the  other  damping  designs  were  effective  at  the  lower  structural 
modes  of  vibration. 

MATERIAL  SELECTION 

The  sensor  frames  are  silver-braised,  316  stainless  steel  while  the  instrument  case 
is  anodized  6061  aluminum.  The  silver  braise  did  not  corrode  in  fresh  water,  but  did  have 
to  be  protected  (in  this  case  by  polyurethane)  in  seawater.  Other  materials  considered  but 
not  chosen  included  3 16L  stainless  steel,  titanium,  and  carbon  fiber.  Titanium  and  3 16L 
stainless  steel  were  not  available  in  the  shapes  needed,  and  titanium  cannot  be  braised 
satisfactorily.  Gun-barrel  drilling  all  the  tubes  in  this  structure  would  have  been  too 
expensive.  The  instrument,  in  general,  was  designed  for  "one  off'  manufacturability  with 
all  shapes  machinable  using  standard  cutters,  and  all  material  available  from  stock  sizes. 
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Carbon  fiber  was  not  used  because  of  its  cost  and  in  this  application  provided 
marginal  advantages.  Carbon  fiber  structures,  when  at  least  vacuum  bagged  during 
curing,  have  better  strength-to-weight  ratios  and  much  better  stiffhess-to-weight  ratios 
than  metal;  however,  to  achieve  this  level  of  performance  requires  all  the  excess  resin  to 
be  squeezed  out.  This  necessitates  the  molds  to  be  "two  and  a  half'  dimensional  in  that 
they  must  have  one  surface  that  moves  as  the  excess  resin  is  removed.  This  reduces 
design  flexibility.  High-strength  carbon  structures  have  to  be  carefully  designed  and  laid 
up  because  of  carbon's  brittleness.  The  primary  motive  in  considering  the  use  of  carbon 
was  to  increase  structural  natural  frequencies,  to  push  up  the  minimum  velocity  of 
strumming.  In  air  wherein  added  mass  is  negligible,  carbon  fiber  structures  can  be  very 
competitive;  but  in  water  the  large  added  mass  lowers  the  natural  frequencies  to  where 
they  would  be  if  the  structure  was  made  of  metal.  The  one  structural  quality  of  composite 
structures  that  is  still  desirable  is  their  larger  structural  damping  than  metal.  However,  as 
the  strength  of  a  composite  material  improves,  its  damping  is  reduced.  Carbon's  small 
additional  structural  damping  did  not  justify  its  significantly  higher  cost. 


TRANSDUCERS  AND  THEIR  MOUNTING 

The  vorticity  meter  uses  piezoceramic  transducers  to  transmit  and  receive  sound 
waves  at  1.75  MHZ.  The  piezoceramic  used  was  Transducer  Products  LTZ-2. 
Specifications  for  alignment  accuracy  can  be  derived  from  modeling  the  transducer  as  a 
vibrating  piston  in  an  infinite  plane  wall  (Dowling  and  Ffowcs  Williams,  1983).  This 
model  gives  a  pressure  at  a  point  given  by  equation  (2-4)  where  p'  is  the  acoustic  pressure, 


p'M  = 
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Pa  is  the  density,  w  is  the  frequency,  a  is  the  transducer  radius,  U0  is  the  velocity  of 
transducer  surface,  R  is  the  distance  to  point  where  the  pressure  is  calculated,  Jx  is  the 
Bessel  function  of  the  order  one,  and  c  is  the  speed  of  sound.  The  relative  beam  patterns 
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for  0.25  inch  (0.63  cm)  and  0.375  inch  (0.953  cm)  diameter  transducers  are  shown  in  Fig. 
2-5.  The  half-power  beam  width  for  the  0.375  inch  transducers  is  2.7  degrees  and  for  the 
0.25  inch  transducers  is  4. 1  degrees.  These  angles  correspond  to  mounting  the 
transducers  to  the  same  tolerance  of  0.018  inch  (0.46  mm).  Because  of  concern  for 
ganging  of  errors,  the  final  mounting  accuracy  specification  was  0.01  inches  (0.25  mm). 
Alignment  of  the  frame  was  achieved  during  braising  by  attaching  a  jig  consisting  of  steel 
rods  that  were  screwed  in  each  acoustic  path  and  fitting  all  the  joints  loosely.  This 
technique  avoided  ganging  of  errors  from  all  the  joints. 

The  transducers  were  mounted  in  polyurethane  resin  on  a  stainless  steel  backbone 
as  shown  in  Fig.  2-6.  The  polyurethane  used  was  Conap  brand  EN-4;  we  chose  it  because 
it  is  acoustically  transparent  in  water.  Insulating  the  transducers  from  the  steel  is  a  spacer 


Fig.  2-5.  Relative  beam  patterns  for  0.25  inch  (0.635  cm)  and  0.375  inch 
(0.953  cm)  transducers  at  1.75  MHZ  calculated  from  eqn  (2-4). 
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of  Dacron  in  polyurethane,  which  was  found  to  meet  the  alignment  accuracy  needed.  The 
thin  Dacron  cloth  was  never  fully  saturated  by  the  resin,  leaving  microscopic  air  bubbles 
behind  the  transducers,  but  this  did  not  cause  any  measurable  problems.  The  only  stainless 
steel-to-urethane  primer  found  to  be  effective  in  seawater  was  Conap  brand  AD-6,  which 
has  an  acid  base.  If  water  diffuses  through  the  polyurethane  jacket,  dissolves  the  acid 
residue,  and  migrates  to  the  transducers;  the  frit  adhering  the  silver  electrodes  to  the 
ceramic  could  be  attacked  by  the  acid  causing  transducer  failure.  The  sensors  have  been 
left  in  freshwater  for  several  months  at  a  time  with  no  degradation,  but  this  potential 
problem  should  be  monitored  on  long-term  deployments.  The  entire  pod  is  encapsulated 
in  one-sixteenth  of  an  inch  of  polyurethane. 


electrode 


transducer 


spacer  of  dacron  saturated 
with  polyurethane 

■molded  polyurethane 
•stainless  steel  backbone 


Fig.  2-6.  Transducer  mounting  on  stainless  steel  backbone 

Piezoceramic  transducers  have  a  high  acoustic  impedance  relative  to  water, 
causing  a  low  acoustic  transmission  coefficient  of  about  0.04.  Quarter  wave  plates  with 
an  acoustic  impedance  equal  to  the  geometric  mean  between  that  of  water  and  ceramic 
were  considered  to  improve  the  transducer  effectiveness,  but  their  additional  cost  of  over 
one  hundred  dollars  extra  per  transducer  for  the  twenty  four  transducers  per  sensor,  made 
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this  option  unattractive.  Instead,  the  BASS  electronics  take  advantage  of  the  high 
transducer  Q  and  transmit  fifteen  cycles  of  the  carrier  wave  at  the  transducer  natural 
frequency.  The  BASS  circuit  times  off  the  fifteenth  cycle  of  the  receiver  transducer, 
whose  energy  is  coherently  summed  from  the  fifteen  cycles.  This  technique  largely 
overcomes  the  low  transducer  acoustic  match  with  water. 

The  minimum  pod  size  was  set  by  transducer  considerations.  Transducers  smaller 
than  0.25  inch  in  diameter  (0.64  cm)  lose  much  of  their  effectiveness  due  to  the  solder 
mass  connecting  the  electrode,  thereby  detuning  that  area  of  the  transducer  and  the 
depolarized  zone  underneath  the  solder.  BASS  electronics  use  a  low-input-impedance 
cascode  receiver  to  reduce  zero  drift  from  cable  capacitance  change,  but  this  low-input 
impedance  coupled  with  the  high  transducer  output  impedance  requires  a  strong  output 
signal  from  the  receiving  transducer. 

The  15  centimeter  path  prototype  vorticity  meter  used  0.25  inch  (0.64  cm) 
diameter  transducers;  these  are  the  smallest  available  with  both  wires  attached  to  one  side 
of  the  transducer.  The  transducer  output  gain  scales  as  the  diameter  to  the  fourth  power, 
requiring  that  its  transmitter  voltage  be  raised  relative  to  a  regular  BASS.  The  prototype 
also  used  EN-4  polyurethane  as  a  pressure  block.  This  is  not  a  hard  polyurethane  and 
could  cold  flow  on  a  long,  deep  deployment.  These  pods  ended  up  having  a  1.0  inch  (2.54 
cm)  outer  diameter. 

The  45  centimeter  path  vorticity  sensors  used  0.375  inch  (0.953  cm)  diameter 
transducers  to  avoid  excessive  transmitter  voltages.  The  acoustic  energy  from  the  receiver 
transducer  scales  as  one  over  the  path  length  squared.  The  electronics  of  this  instrument 
used  a  three-to-one  turns  ratio  in  the  toroidal  core  transformers  that  power  the  transducers 
(Williams  et  al,  1987).  In  these  sensors,  thermistors  (Yellow  Springs  Instrument  Co. 
44030)  were  mounted  inside  the  pod  arm  tubes  at  a  distance  of  2  inches  (5  cm)  from  the 
pods  on  each  side  of  each  sensor.  The  thermistors  measure  temperature  and  temperature 
stratification.  For  greater  depth  capability,  the  transducer  pods  use  glass-filled  epoxy 
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pressure  blocks.  These  pods  have  a  1.25  inch  (3.17  cm)  outer  diameter. 

The  1.5  meter  path  vorticity  meter  briefly  mentioned  uses  0.75  inch  (1.9  cm) 
diameter  transducers  and  operates  at  the  lower  frequency  of  875  KHz  to  keep  transducer 
alignment  from  becoming  too  difficult;  it  is  difficult  to  align  objects  to  less  than  one 
degree. 


BUOY  DESIGN 

A  shear  measuring  buoy  was  built  of  two  45-centimeter  path,  vorticity  sensors, 
their  electronics,  an  inertial  measurement  unit,  and  a  float  with  recovery  gear.  Fig.  2-4. 

The  buoy  was  designed  to  minimize  flow  disturbance  in  the  measurement  volumes  by 
moving  with  ocean  swell,  minimize  surface  expression  (the  float),  minimize  sensor  cylinder 
size,  and  locating  the  electronics  package  well  below  the  measurement  volumes. 
Measurement  volumes  are  centered  at  0.83  meters  and  2.45  meters  below  the  water 
surface,  and  the  total  buoy  height  is  5.06  meters.  The  buoy  has  a  strapdown,  inertial 
measurement  unit  to  measure  buoy  motion  and  remove  it  from  relative  flow,  and  to  rotate 
the  measurements  into  an  Earth  reference  system.  The  inertial  measurement  unit  is 
described  in  the  electrical  design  section  and  the  inertial  processing  is  described  in  the 
signal  processing  section. 

The  buoy  is  modular  and  is  bolted  together  with  interchangeable  sensors.  A 
bumper  of  rolled,  welded  stainless  steel  tubing  is  mounted  on  the  instrument  case  so  that 
the  buoy  can  be  laid  on  its  side  on  a  flat  surface  between  this  bumper  and  the  top  spherical 
float,  and  not  bend  the  pod  arms.  The  top  float  is  made  of  expanded  PVC  foam  with  a 
fiberglass  skin.  This  system  is  light,  relatively  incompressible,  failsafe,  shapable,  and 
corrosion  free. 

The  buoy  is  deployed  freely,  drifting  with  the  currents,  requiring  reliable  recovery 
aids.  The  recovery  systems  include  an  ARGOS  transmitter  that  transmits  to  a  satellite 
which  measures  and  relays  buoy  location,  a  strobe  light,  a  VHF  radio  transmitter,  and 
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when  deployed  below  the  ocean  surface,  an  acoustic  transponder  that  is  compatible  with  a 
diver-operable,  subsea  direction  finder. 

To  measure  internal  wave  shear  in  the  thermocline,  the  buoy  can  be  deployed  using 
an  elastic  tether  between  the  top  float  and  the  rest  of  the  instrument.  Fig.  2-7.  The  drifter 
uses  a  bungy  cord  strung  inside  a  hollow  sleeve  rope  that  is  slack,  to  isolate  float  heave 
from  the  instrument.  The  chain  on  the  float  bottom  and  poly  fishnet  floats  at  the 
instrument’s  top  keep  tension  on  the  elastic  tether  below  two  kilograms,  allowing  a  low- 
spring-rate  soft  bungy  to  be  used.  Observations  of  the  instrument  during  dives  showed 
very  effective  heave  suppression. 


Fig.  2-7.  Acoustic  vorticity  meter  configured  as  a  thermocline  shear  measuring 
drifter. 
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2.B.  ELECTRICAL  DESIGN 


ELECTRICAL  OVERVIEW 

Modified  BASS  electronics  (Williams  et  al,  1987)  form  the  heart  of  the  vorticity 
meter.  The  physical  electronics,  inertial  measurement  unit  and  data  logger  are  shown  in 
Fig.  2-8.  A  block  diagram  of  the  electrical  system  is  shown  in  Fig.  2-9.  The  BASS 


Fig.  2-8.  Photo  of  electronics  on  right,  inertial  measurement  unit  in  middle, 
battery,  and  data  logger  at  left  with  the  instrument  case  above. 


electronics  are  controlled  by  a  Tattletale  5  computer  made  by  Onset  computers,  Pocasset, 
MA..  The  Tattletale  5  computer  directs  each  pair  of  acoustic  transducers  to  measure  water 
velocity  along  their  acoustic  path,  polls  all  the  analog  inputs  from  the:  thermistors, 
accelerometers,  rate  gyros,  and  magnetometers;  and  sends  this  data  serially  to  a  data 
logger.  The  data  logger  is  a  Tattletale  6  computer  that  records  the  data  on  a  hard  disk 
drive  after  buffering  the  data  in  semiconductor  memory.  This  “buffering  to  memory”  is 
done  to  save  power;  the  hard  disk  is  only  spun  up  and  written  to  when  the  memory  is  full. 
The  analog  signal  conditioning  done  to  the  inertial  sensors'  signals  is  scaling  voltages  and 
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BASS 


Fig.  2-9.  Block  diagram  of  electrical  system 

low-pass,  anti-alias  filtering.  Each  of  the  recovery  devices:  the  ARGOS  transmitter,  VHF 
transmitter,  and  strobe,  has  its  own  battery  for  redundancy.  The  recovery  devices  are  not 
shown  in  this  diagram. 

One  problem  with  the  data  logger  system  is  the  high  current  required  to  spin  up  the 
hard  disk  before  a  software  timeout  is  reached.  The  main  batteries  of  this  instrument  are 
alkaline  which  have  a  high  output  impedance  when  cold,  and  are  unable  to  start  the  hard 
disk  below  a  certain  temperature.  To  resolve  this  problem,  a  lead  acid  gel  cell  with  lower 
output  impedance  at  low  temperature,  is  used  as  a  capacitor  to  start  the  disk.  This  gel  cell 
is  then  recharged  by  the  alkaline  main  batteries  between  disk  writes. 

The  longer  acoustic  path  length  required  modifying  the  transmit  voltage  and  timing 
of  the  BASS  electronics.  To  increase  the  transmit  voltage,  the  transformer  cores  that 
drive  the  acoustic  transducers  were  wound  with  a  three-to-one  ratio  instead  of  the  one-to- 
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one  ratio  in  a  regular  BASS.  The  fast  timing  to  control  the  acoustic  velocity  measurement 
is  done  in  analog  on  the  timing  and  burst  generator  card  to  save  power.  The  instrument 
transmits  fifteen  cycles  at  1.75  MHz.  for  8.57  ps.  The  receiver  is  turned  on  282.5  ps  after 
the  start  of  the  transmit;  this  corresponds  to  the  fastest  measurable  sound  speed  of  1590 
m/s.  The  receiver  effectively  turns  off  344.5  ps  after  the  transmit  start;  this  corresponds 
to  the  slowest  measurable  sound  speed  of  1370  m/s.  The  acoustic  path  velocity  is  then 
remeasured  383.5  ps  later  with  the  electronics  reversed  to  cancel  out  electronic  drift.  The 
measurement  cycle  consists  of  a  time  stamp  of  four  bytes,  24  forward  and  reverse  velocity 
measurements  that  take  18.4  ms,  13  analog  measurements  of  the  thermistors  and  inertial 
sensors,  and  then  sending  the  eighty  byte  measurement  sequence  at  9600  baud  to  the 
logger,  taking  83  ms.  The  whole  measurement  sequence  is  repeated  every  150  ms. 


The  data  format  is  shown  in  Table  2-1 .  All  the  data  is  transmitted  and  stored  in 
unsigned  binary  integers  except  the  velocity  measurements  which  are  in  two's  complement. 
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hexadecimal  record  header 
time  stamp  hr:min:s:counter 
four  acoustically  measured  velocities  making 
up  a  circulation  square 


thermistor  output 
accelerometer 
rate  gyro 
magnetometer 


Table  2-1.  Vorticity  meter  data  format.  The  numbers  in  brackets  indicate  the 
length  of  each  variable  in  bytes 
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INERTIAL  MEASUREMENT  UNIT 

A  strapdown  inertial  measurement  unit  was  built  into  the  shear  measuring  buoy  to 
measure  buoy  motion  allowing  the  motion  to  be  removed,  in  processing,  from  flow 
measurements.  The  Inertial  Processing  Background  section  shows  that  it  is  necessary  to 
measure  and  compensate  for  buoy  motion  even  if  only  time  average  measurements  from  a 
buoy  are  desired.  A  strapdown  inertial  measurement  unit  was  chosen  over  a  gimballed 
inertial  measurement  unit  to  save  power  and  size.  A  gimballed  inertial  measurement  unit, 
consisting  of  accelerometers  on  a  gyro-stabilized  platform  that  does  not  rotate  with 
respect  to  inertial  space,  is  more  accurate  than  a  strapdown  system  consisting  of  the 
accelerometers,  rate  gyros  and  magnetometers  mounted  "strapped  down"  to  the  buoy 
(VanBronkhorst,  1978).  Gimballed  systems  are  expensive,  delicate,  and  consume  much 
power;  whereas  strapdown  systems  are  inexpensive,  more  robust  mechanically,  and 
consume  less  power.  Strapdown  systems  do,  however,  require  more  computation. 
Gimballed  systems  use  rate  gyros  as  nulling  sensors,  do  not  expose  the  rate  gyros  to  the 
angular  rates  of  the  buoy,  and  therefore  can  measure  angular  rates  more  accurately. 
Strapdown  systems,  on  the  other  hand,  expose  their  rate  gyros  to  the  frill  angular  rates  of 
the  buoy,  requiring  rate  gyro  accuracy  and  linearity  over  a  wide  range.  Gimballed  systems 
are  able  to  cancel  out  some  error  terms,  such  as  some  small  accelerometer  misalignment, 
that  strapdown  systems  can  not  cancel  out  (Schmidt,  1978).  Some  wavebuoys  use  a  low- 
power  version  of  a  gimballed  system  by  floating  a  large  sphere  with  accelerometers,  in  oil, 
and  limiting  the  system's  righting  response  to  frequencies  lower  than  the  wind-wave 
spectrum.  These  wavebuoy  systems  are  however,  large  and  heavy.  The  shear  measuring 
buoy  uses  a  strapdown  inertial  measurement  unit  to  save  power  (battery  weight,  size,  and 
flow  disturbance),  to  avoid  the  bulk  and  weight  of  some  of  the  wavebuoy  inertial 
measurement  platforms,  and  to  save  money. 

The  strapdown  inertial  measurement  unit  consists  of  a  three-axis  accelerometer, 
three  single-axis  rate  gyros,  and  a  three-axis  magnetometer.  These  sensors  and  some  of 
their  specifications  are  listed  in  Table  2-2.  It  is  worthwhile  to  note  that  these  are  not 
inertial-navigation-grade  sensors.  For  example,  good  inertial  grade  rate  gyros  are  106 
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Three-axis  accelerometer  Columbia  Triaxial  Accelerometer 

model  SA-307  HPTV 

X  -Ito+IG 

Y  -Ito+IG 

Z  0  to+2G 

case  alignment  +/-  0.5° 

Three  single-axis  rate  gyros  Systron  Donner  Gyrochip  Angular  Rate  Sensor 
X,  Y,  and  Z  -50  to  +50  deg/s 

bandwidth  >  60  Hz. 
scale  calibration  1  % 
linearity  <0.05%  of  full  scale 
input  power  noise  requirement  <  0.01  v  rms  and 
<  0.001  v  rms  at  8.7  Khz.  +/-  500  Hz. 

Three-axis  magnetometer  Develco  model  9200  fluxgate  magnetometer 
X,  Y,  and  Z  -600  to  +600  mGauss 
alignment  +/-  1° 
output  ripple  0.4%  of  full  scale 
temp  stability  <3  %  from  0°C  to  60°C  sensitivity 
<1%  zero 

linearity  +/-  0.5%  of  full  scale 

Table  2-2.  Description  of  the  inertial  measurement  unit  sensors 


times  more  accurate  than  the  low-power  rate  gyros  used  in  the  vorticity  meter.  The 
alignment  of  the  sensors  used,  inside  their  sensor  cages,  is  only  good  to  about  one  degree. 
Techniques  such  as  Schuler  tuning  are  not  useful  when  the  rate  gyro  noise,  even  when 
averaged  over  ten  minutes,  and  drift  are  ten  times  the  earth's  rotation  rate.  However,  with 
appropriate  processing,  as  is  shown  in  the  Inertial  Processing  section,  these  sensors  are 
adequate  for  this  application. 

The  inertial  measurement  system  requirements  will  now  be  specified.  Errors  in 
estimated  buoy  attitude  and  motion  should  not  cause  more  than  five  percent  error  in  flow 
measurements.  This  translates  to  an  angle  error  specification  of  three  degrees.  Angle 
rates  in  sensor  coordinates  should  be  better  than  one  percent  to  compensate  for  buoy 
motion.  The  Inertial  Processing  Background  section  shows  that  the  heading  specification 
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of  three  degrees  actually  requires  pitch  and  roll  error  be  less  than  one  degree.  The  inertial 
system  should  be  self-aligning  in  a  seaway,  because  boats  that  the  buoy  will  be  launched 
from  do  not  have  three-axis  inertial  reference  for  a  conventional  system  alignment.  Errors 
from  coning  motions  should  not  be  significant.  The  Inertial  Processing  Background 
section  discusses  coning  motions  and  their  importance.  While  these  specifications  seem 
lax  by  inertial  navigation  standards,  they  are  adequate  for  removing  buoy  motion  from  the 
buoy  measurements  of  vorticity.  With  the  sensor  specifications  given,  it  would  not  be 
easy  to  do  significantly  better. 

CALIBRATION  OF  THE  INERTIAL  MEASUREMENT  UNIT 

All  calibrations  of  the  inertial  sensors  were  performed  with  the  sensors  inside  the 
buoy  instrument  case,  with  all  batteries,  the  same  electronics,  and  A/D  converter  that  were 
used  when  the  instrument  was  deployed  at  sea.  The  rate  gyros  were  calibrated  by  rotating 
the  instrument  case  in  a  rotating  welding  fixture  and  on  a  rotating  table,  at  different 
speeds,  and  around  each  axis.  The  accelerometers  and  magnetometers  were  calibrated 
statically  by  holding  the  case  at  different  roll,  pitch,  and  heading  angles  in  the  middle  of  an 
open  field.  The  magnetometer  was  calibrated  in  an  open  field  to  avoid  the  magnetic 
anomalies  in  the  laboratory.  The  angular  positions  were  measured  on  a  piece  of  leveled, 
one-inch  plywood  for  stability;  rested  on  measured  wedges  and  were  calibrated  with 
respect  to  a  single  axis  magnetometer  and  to  a  level.  In  all  these  calibrations,  the 
averaging  time  for  each  data  point  was  greater  than  five  minutes  to  ensure  reliable 
statistics.  The  calibration  results  are  shown  in  Figs.  2-10  and  show  good  linearity  of  the 
sensors,  except  at  the  very  ends  of  their  measurement  limits.  The  noise  from  all  except 
one  sensor  channel  was  white  to  the  background  levels  of  vibration  in  the  Bigelow 
building  basement.  The  Y  channel  of  the  rate  gyro  however,  had  a  17  second  oscillation 
that  corresponds  to  a  0.5  degree  rocking  motion  (Fig.  2-11).  The  cause  of  the  noise  is 
unknown,  but  is  within  the  manufacturer's  sensor  specification.  The  standard  deviations  of 
the  sensor  noise  are  listed  in  Table  2-3. 
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Fig.  2-10  c.  Calibration  of  Magnetometers.  The  top,  middle,  and  lower  graphs  are 
the  calibrations  of  the  x-axis,  y-axis,  and  z-axis  magnetometer  channels. 


AXES 

SENSOR  X  Y  5 

accelerometer  (m/s2)  0.013  0.012  C 

rate  gyro  (rad/s)  0.0014  0.0025  0 

magnetometer  (mG)  0.77  0.91  0 


Table  2-3.  Standard  deviation  of  the  inertial  sensor  noise 


2.C.  SIGNAL  PROCESSING 

This  section  reviews  the  kinematics  of  rotation,  develops  an  algorithm  to  calculate 
buoy  motion,  describes  the  total  data  processing,  and  simulates  this  inertial  signal 
processing  to  estimate  its  accuracy  . 

INERTIAL  PROCESSING  BACKGROUND 

COORDINATE  TRANSFORMATIONS 

This  subsection  reviews  the  kinematics  of  rotation  and  introduces  some  of  the 
mathematical  conventions  used  in  this  thesis.  The  rotation  of  coordinate  system  xyz  into 
reference  frame  x’y’z’  is  shown  in  Fig  2-12.  The  transformation  of  vector  coordinates 
from  one  reference  frame  to  the  other  can  be  described  by  a  matrix  multiplication  (2-5) 
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Fig.  2-12.  Rotation  of  coordinate  system  xyz  to  x'y'z' 
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(Crandall  et  al,  1968).  In  this  equation  the  vector  components  x,y,z  are  the  coordinates  of 
vector  fin  the  xyz  frame  and  the  primed  x'y'z'  vector  components  are  the  coordinates  of 
the  vector  in  the  primed  frame.  This  rotational-transformation  matrix  A  is  made  of 
direction  cosines  (2-6).  The  rotation  matrix  has  six-fold  redundancy  in  direction  cosines 

Ax'y  =  Ux‘  fy  (2-6) 

as  there  are  six,  independent  relations  among  the  nine  elements.  Some  properties  of  the 
matrix  are  that  its  inverse  is  its  transpose  and  its  determinant  is  plus  one  (2-7).  The  same 


[A]"1  =  [A]' 
\A\  =  -1 


(2-7) 


rotation  of  the  vector  f  referenced  in  a  constant  reference  frame  can  be  described  by  the 
same  equation  (2-5)  but  using  the  inverse  of  A.  Successive  rotations  are  not  commutative; 
a  rotation  A  about  the  x  axis  followed  by  a  rotation  B  about  the  y  axis  does  not  in  general 
equal  the  rotation  B  about  the  y  axis  followed  by  the  rotation  A  about  the  x  axis. 
Successive  rotations  can  be  described  as  matrix  multiplication  of  the  individual  rotation 
matrices,  and  matrix  multiplications  are  not  commutative  (2-8).  Because  of  the 

IAMB]  *  [B]*[A]  (2-8) 


noncommutivity  of  rotations,  finite  rotations  cannot  be  represented  by  vectors  that  follow 
the  rules  of  vector  addition. 

Infinitesimal  rotations  can,  however,  be  treated  as  vectors.  The  error  in  ignoring 
the  order  of  small  finite  rotations  is  proportional  to  the  angle  of  rotation  squared.  Angular 
velocities  can  be  treated  as  vectors  with  no  error. 
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EULER  ANGLE  CHOICE 

This  instrument  development  uses  a  3,2,1  Euler  angle  system  to  calculate  buoy 
rotations  (Goldstien,  1981).  In  a  3,2,1  Euler  angle  system,  the  object  or  reference  frame  is 
rotated  firstly  about  its  z  axis  by  an  angle  4>  (the  three  index),  secondly  about  its  new  y 
axis  by  an  angle  0  (the  two  index),  and  thirdly  about  its  new  x  axis  by  an  angle  ijr  (the  one 
index).  Fig.  2-13.  These  sequential  rotations  can  be  combined  to  form  a  single 


Fig.  2-13.  Order  of  rotations  for  a  3,2,1  Euler  angle  system. 
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rotational-transformation  matrix  (2-9).  The  body-centered  rotation  rates  can  be  computed 


cos((>  sin(J>  0 
D  =  -siruj)  coscj)  0 
0  0  1 

cos0  0  -sin0 
C  =  0  1  0 

sin0  0  cos0 
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1  0  0 
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from  the  Euler  angles  and  Euler  angle  rates  by  (2-10).  This  matrix  can  be  inverted  to  give 
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(2-10) 


the  Euler  angle  rates  in  terms  of  the  Euler  angles  and  the  body-referenced  angle  rates  (2- 
1 1).  This  new  matrix  will  be  called  the  Euler-angle  update  matrix  in  this  thesis  .  Knowing 
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the  initial  attitude,  successive  attitudes  can  be  stepped  through  with  measured  body- 
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centered  rotation  rates. 


The  Euler-angle  update  matrix  has  some  special  properties.  It  should  be  noted  that 
the  F  and  F"1  matrices  are  not  rotational-transformation  matrices;  their  determinants  are 
not  equal  to  one  and  their  inverse  is  not  their  transpose.  There  is  a  singularity  in  the 
Euler-angle  update  matrix  when  0  equals  plus  or  minus  tr/2,  when  the  buoy  is  on  its  side. 
There  is  no  heading  (j>  in  the  Euler-angle  update  matrix.  In  normal  operation,  the  heading 
angle  can  be  much  larger  than  either  the  pitch  or  the  roll  angle,  and  the  calculation  of 
heading  is  not  as  well  conditioned  (as  will  be  shown  in  the  magnetic  compass  section)  as 
the  calculation  of  pitch  or  roll.  This  3,2,1  Euler  angle  system  was  chosen  to  minimize 
error  propagation  and  to  have  its  singularity  at  an  attitude  that  the  buoy  has  not  been  seen 
to  approach  in  service.  Also,  if  the  Euler-angle  update  matrix  is  linearized,  the  error  in 
calculating  Euler  angle  rates  is  smaller  than  with  other  Euler  angle  systems. 

Real  buoy  motions  are  smooth  and  only  approximated  by  the  discrete  calculation 
steps  in  this  algorithm.  Commutivity  error  results  from  the  difference  between  the 
smooth,  real  motion  of  the  buoy  and  the  discrete  steps  assumed  by  any  Euler  angle  system. 
As  mentioned  above  in  this  section,  the  difference  between  the  3,2,1  angle  system  and  a 
1,2,3  angle  system  is  proportional  to  the  rotation  angle  squared.  Average  body-centered 
rotation  rate  is  measured  at  each  time  step  and  used  in  the  algorithm.  The  error  in  each 
step  of  the  algorithm  should  also  be  proportional  to  the  angle  squared  in  that  time  step. 
Increasing  the  calculation  step  rate  increases  the  number  of  errors  with  the  number  of 
steps,  and  decreases  the  error  of  each  step  by  one  over  the  number  squared.  The  result  of 
both  is  a  total  commutivity  error  for  a  buoy  evolution,  that  is  proportional  to  the  angular 
displacement  of  each  time  step. 

OTHER  TECHNIQUES  TO  CALCULATE  ATTITUDE 

There  exist  several  other  choices  of  algorithms  to  calculate  the  attitude  of  a 
strapdown  inertial  system,  including:  direct  direction-cosine-matrix  updates,  four- 
parameter  techniques,  and  other  Euler  angle  systems  (McKern,  1968).  The  direct 
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direction-cosine-matrix  update  and  the  four-parameter  systems  such  as  the  Euler 
parameter  technique  have  the  advantages  of  no  singularities  and  fewer  calculations  with 
transcendental  functions.  Both  techniques  are  popular  in  real-time  inertial-navigation 
systems  that  require  an  unrestricted  operational  envelope,  and  that  may  have  limited 
processing  power.  As  mentioned  earlier,  the  singularity  of  the  chosen  3,2, 1  Euler  angle 
system  is  an  attitude  that  the  buoy  has  not  approached  in  the  ocean,  therefore  the 
existence  of  this  singularity  is  not  a  problem. 

The  signal  processing  for  this  instrument  was  performed  after  the  buoy  was 
recovered,  on  an  80486  computer  in  the  Matlab  environment.  Before  deciding  which 
inertial  attitude  computation  technique  to  use,  a  test  was  performed  to  measure  how  much 
transcendental  functions  would  slow  the  computation.  A  series  of  tangents  and  adds  of 
vectors  of  random  numbers  took  less  than  twice  as  long  as  the  same  number  of  multiplies 
and  adds  of  the  same  vectors.  I  do  not  know  if  this  is  more  a  function  of  Matlab's 
overhead  or  efficient  processing  of  transcendental  functions  by  the  floating  point 
processor.  There  is  no  significant  increase  in  processing  speed  when  using  these  less 
intuitive  algorithms  in  this  processing  environment.  Additionally,  Euler  angles  use  less 
memory.  On  a  more  general  inertial-measurement  system,  these  other  techniques  would 
have  to  be  considered. 

Two  other  Euler-angle  systems  that  I  considered  but  rejected  were  the  3, 1,3 
system  and  the  1,2,3  system  (Pio,  1966).  The  3,1,3  system  is  used  in  many  dynamics 
textbooks  and  is  convenient  for  describing  the  motion  of  a  spinning  top.  The  order  of 
rotations  is  about  the  z  axis,  the  new  x  axis,  and  then  the  new  z  axis.  The  3,1,3  Euler- 
angle  update  matrix  however,  has  a  singularity  at  the  operating  point  of  floating  right  side 
up,  and  was  therefore  dismissed.  The  1,2,3  Euler-angle  system  rotates  about  the  x  axis, 
then  about  the  new  y  axis,  and  then  about  the  new  z  axis.  This  system,  however,  has  the 
heading  angle  4>  in  the  Euler-angle  update  matrix.  The  heading  ranges  from  0  to  360 
degrees  while  pitch  and  roll  seldom  exceed  20  degrees.  Additionally,  as  the  next  section 
explains,  the  measurement  of  heading  is  poorly  conditioned  as  compared  to  the 
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measurement  of  pitch  and  roll,  resulting  in  greater  heading  error.  The  3,2,1  Euler-angle 
system  results  in  smaller  errors  in  buoy  attitude  processing  than  either  the  3,1,3  or  the 
1,2,3  systems,  and  was  therefore  chosen. 

MAGNETIC  COMPASSES 

Magnetic  compasses  work  by  projecting  the  magnetic  lines  of  flux  on  a  plane  that 
is  perpendicular  to  the  perceived  gravitational  field.  In  Woods  Hole  MA.,  the  magnetic 
lines  of  flux  are  inclined  from  the  horizontal  by  about  seventy  degrees,  making  this 
measurement  poorly  conditioned.  This  is  why  aircraft  need  gyrocompasses  to  measure 
their  heading  in  a  turn.  The  vorticity  meter  has  a  three  axis  fluxgate  magnetometer.  The 
inertial  processing  of  the  instrument,  described  in  the  Inertial  Algorithm  section,  rotates 
the  measured  magnetic  fluxes  into  a  horizontal  plane  and  then  computes  heading.  An 
error  in  pitch  or  roll  in  the  algorithm  used  is  equivalent  to  having  the  magnetometer  at  the 
attitude  of  the  pitch  or  roll  error,  which  can  result  in  a  heading  error  that  is  greater  than 
the  pitch  or  roll  error. 

This  section  calculates  an  example  of  what  happens  if  instrument  heading  is 
calculated  from  magnetometer  measurements  made  in  a  plane  that  is  tilted  and  the  tilt  is 
not  accounted  for.  The  following  derivation  uses  the  coordinate  system  shown  in  Fig  2- 
14.  With  this  coordinate  system,  the  heading  of  a  horizontal  plane  can  be  computed 
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Fig.  2-14.  Coordinate  system  used  in  compass  example 
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from  the  inverse  tangent  of  the  x-direction  flux  divided  by  the  y-direction  flux.  If  the  plane 
of  the  x  and  y  magnetometer  is  put  at  an  angle  to  horizontal,  the  measured  fluxes  will  be 
given  by  equation  (2-12).  In  this  equation,  Bx,  By,  and  Bz  are  the  magnetic  fluxes  in  the 


By  =  [A] j  540 mG  cos(-70°)  • 
Bz  540 mG  sin(-70°) 


(2-12) 


rotated  (buoy  attitude  direction)  x,y,z,  and  [A]  is  the  rotational-transformation  matrix 
described  earlier.  The  measured  heading  will  then  be  given  by  equation  (2-13).  In  this 
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equation,  4>t  is  the  true  heading,  (j)m  is  the  measured  heading.  The  instantaneous  heading 
error,  at  some  headings,  for  one  degree  of  pitch  or  roll  is  2.75°.  An  accuracy  specification 
of  three  degrees  in  heading  will  require  pitch  and  roll  accuracy  of  one  degree. 

Because  of  the  nonlinearities  in  equation  (2-13),  taking  the  mean  of  measured 
heading  of  a  rocking  buoy,  and  not  compensating  for  pitch  and  roll,  can  result  in 
significant  error.  Fig.  2-15  shows  the  results  of  a  fifteen  degree  buoy  oscillation  in  0  when 
the  true  heading  ((j>t )  is  55  degrees.  The  rocking  angle  0,  the  true  heading,  and  the 
measured  heading  are  shown  through  the  rocking  cycle.  The  error  in  mean-measured 
heading  is  over  eight  degrees.  Not  knowing  which  way  is  up  can  result  in  significant  error 
in  measured  heading,  even  in  the  mean. 
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Measured  heading  from  15  degree  buoy  oscillation  in  theta  with  heading  constant 
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Fig.  2-15.  Measured  heading  for  rocking  buoy  if  heading  is  calculated  from 

tan  'B/By  and  pitch  and  roll  are  not  accounted  for. 

CONING  MOTIONS 

One  motive  for  measuring  and  compensating  for  instrument  attitude  is  to  avoid  the 
measurement  errors  that  can  result  from  coning  motions  (Goodman  and  Robinson,  1958). 
Coning  motion  occurs  when  a  line  in  the  object  moves  in  space  as  to  follow  the  surface  of 
a  cone.  This  section  describes  one  such  coning  motion  and  calculates  the  error  that  results 
if  the  motion  is  not  accounted  for.  The  Euler  angles  in  this  example  have  a  time  history 
described  by  equation  (2-14),  where  ij/  is  a  sine  curve  and  0  is  a  cosine  curve.  The 
measured  body-referenced  angle  rates  are  given  by  equation  (2-15).  Plugging  in  the 
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(2-14) 
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values  given,  assuming  small  coning  motion  angles,  and  taking  the  small  angle 
approximation  of  sin(x)=x  and  cos(x)=l,  results  in  equation  (2-16).  The  body-referenced 


a)  =  au>  cos(to  t) 

<*y  =  -«wsin(co/)  (2-16) 

o  =  a2(s>  sin2(o)  t) 

z  axis  measures  a  time-averaged  rotation  rate  while  the  true  heading  (j)  is  zero  and 
constant.  If  the  inertial  measurement  system  does  not  keep  up  with  buoy  motion,  this 
error  can  result.  If  a  buoy  measures  vorticity,  ignores  instantaneous  attitudes  and  angle 
rates,  and  only  compensates  measured  relative  vorticity  for  average  drift  in  heading,  it 
would  calculate  a  earth-referenced  vorticity  of  minus  two  times  this  angle  rate  (coj  in  a 
still  fluid. 

ALGORITHM  ERROR  PROPAGA  TION 

Because  the  Euler  angle  rates  used  to  calculate  future  Euler  angles  depend  on  the 
current  Euler  angles,  an  error  in  the  current  estimate  of  Euler  angles  can  propagate  into 
the  future.  This  section  describes  a  perturbation  expansion  of  the  Euler  angle  update 
algorithm  to  estimate  how  errors  propagate.  Equation  (2-17)  is  the  continuous  limit  of  the 

*  =  [F"1('F)]{<d}  =  [G(Y)]{<o}  (2-17) 
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Euler-angle  update  algorithm.  The  Euler  angle  is  broken  into  the  true  Euler  angle  Y  and  a 
small  error  6  equation  (2- 1 8)  and  plugged  into  equation  (2- 1 7)  to  give  equation 
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(2-19).  This  equation  is  solved  for  how  the  error  changes  using  a  Taylor  expansion  for  G 


¥+8  =  [G('F+8)]{to}  (2-19) 

and  dropping  second  and  higher  order  terms  in  error  (2-20).  From  equation  (2-20)  we  can 
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estimate  how  a  small  error  propagates  during  a  specified  body  motion.  We  can  also 
conclude  that  errors  in  the  heading  Euler  angle  d ,  do  not  propagate  to  roll  or  pitch 
estimates  because  the  heading  <j>  is  not  in  the  Euler-angle  update  matrix. 

The  first  example  coning  motion  to  be  discussed  will  be  the  coning  motion 
example  from  Coning  Motions  section.  The  angles  and  angle  rates  are  reviewed  in 
equation  (2-21).  In  this  equation,  a  is  the  angular  amplitude  and  cd  is  the  frequency  of  the 


i| !  =  a  sin(cof) 
0  =  a  cos(cof) 
<j>  =  0 


vj;  =  dcocos(cDf) 
0  =  -ao)sin(o)f) 

<j>  =  0 


CD  =  tfCDCOS(cDf) 

CD,  =  -dCDsin(cDf)  +  o(a3)  (2-21) 

cd2  =  fl2CDsin2(tDf) 
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coning  motion.  The  last  relations  are  good  to  order  a2  in  that  a  small  angle  approximation 
was  made.  Using  these  parameters  in  equation  2-20,  making  the  small  angle 
approximation,  and  dropping  terms  that  have  a  time  average  of  zero  results  in  equation  (2- 
22)  showing  the  Euler  angle  errors  stay  constant.  Dropping  the  zero  time-average  terms 

b  =  0 

c  =  0  +  o(a3)  (2-22) 

d  =  o 

is  reasonable  because  under  these  assumed  motions,  the  errors  b,  c,  and  d  will  change 
much  slower  than  co  by  at  least  order  a3. 

The  second  example  of  coning  motion  is  one  with  zero  time-average,  buoy- 
referenced,  z-axis  rotation  rate.  The  angles  for  this  motion  are  shown  in  equation  (2-23). 

i|r  =  a  sin  (cut)  ijf  =  tfwco  s(o)t)  wx=  ^cos(ot)  +o(a3) 

Q  =  a  cos(wt)  0  =  -acosin(o)f)  (o^=  -«tosin(o)t)  +o(a3)  (2-23) 

<j>  =  fa  2a)sin2((Dt)  dt  <j>  =  ^2o)sin2(o)f)  co  =  0  +o(a 3) 

Plugging  this  motion  into  equation  (2-20),  making  the  small  angle  approximation,  and 
dropping  terms  of  order  a3  and  smaller  gives  the  angular  error  propagation  (2-24).  These 

b  =  -£tf2cosin2((of) 

c  =  foi2<osin2(cot)  (2-24) 

d  =  0 

equations  were  calculated  for  an  initial  error  of  0.05  radians  in  i|r  and  the  results  are  shown 
in  Fig.  2-16.  This  solution  shows  that  under  certain  buoy  motions,  slow  oscillating  errors 
in  Euler  angles  can  occur. 
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INERTIAL  ALGORITHM 

This  section  describes  the  algorithm  used  to  compute  the  buoy  attitude  and  motion 
that  are  in  turn  used  to  remove  buoy  motion  from  measured  fluid  motion.  Complementary 
filtering  allows  the  use  of  noisy  sensors  with  drift  and  still  meet  the  accuracy  specifications 
mentioned  in  the  Inertial  Measurement  Unit  section.  Using  more  accurate  sensors  would 
increase  cost  and  require  more  power.  The  rate  gyros  have  a  manufacturer  drift 
specification  of  0. 1  degrees  per  second.  To  meet  the  three-degree  accuracy  specification, 
a  conventional  strapdown  algorithm  that  just  integrates  rate  gyro  output,  would  only  meet 
the  accuracy  specification  for  thirty  seconds.  The  use  of  these  noisy,  drifty  rate  gyros 
requires  information  from  other  sensors  to  compensate  for  gyro  drift.  The  algorithm  I  use 
to  compute  buoy  attitude  and  motion  makes  assumptions  about  the  buoy  motion  that 
allow  the  use  of  accelerometers  and  magnetometers  to  compensate  for  gyro  drift.  For  a 
general  use  inertial  system,  such  as  one  for  an  aircraft,  these  motion  constraints  could  not 
be  made  and  this  algorithm  could  not  be  used. 

CONSTRAINTS  OF  BUOY  MOTION 

Constraints  on  buoy  motion  allow  the  inertial  system  to  align  itself  in  a  seaway. 
These  motion  assumptions  are  based  on  the  measured  information  that  there  is  almost  no 
wave  energy  below  a  twenty-five  second  period.  This  analysis  ignores  tides  because  tides 
would  not  tilt  the  buoy  and  the  inertial  measurement  unit  is  not  accurate  enough  to 
measure  accelerations  at  a  frequency  as  low  as  a  tidal  frequency.  The  motion  assumptions 
for  periods  longer  than  twenty-five  seconds  are  zero  average  heave  (vertical 
displacement),  zero  average  heave  velocity,  and  zero  average  horizontal  velocity  (surge 
and  sway).  Currents  on  the  order  of  several  tenths  of  a  meter  per  second  that  vary  over 
periods  of  hours  or  days  cannot  be  measured  by  this  inertial  measurement  unit,  and  are 
assumed  to  be  zero. 

COMPLEMENTARY  FIL  TERS 

The  algorithm  uses  complementary  filters  to  mix  redundant  estimates  of  Euler 
angles  from  rate  gyros,  accelerometers,  and  magnetometers.  A  block  diagram  of  a 


60 


complementary  filter  is  shown  in  Fig.  2-17.  In  this  figure,  s(t)  is  the  true  signal,  n,  or2(t) 


Fig.  2-17.  Complementary  filter  of  redundant  sensor  measurements 

are  noise  added  to  the  signal,  G(<o)  is  the  transfer  function  of  filter  one,  and  x(t)  is  the 
output.  Weiner  filter  theory  gives  the  result  equation  (2-25)  for  an  optimal  non-causal 
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GH  =  - —  (2-25) 

s  +  s 

nlnl  n2n2 

filter  (Brown,  1983).  In  this  equation,  Sn  is  the  power  spectral  density  of  noise  in  signal 
I.  This  solution  assumes  that  noise  from  the  two  sensors  n,  and  n2  are  uncorrelated.  Non- 
causal  filters  were  chosen  over  causal  filters  because  the  processing  does  not  need  to  be 
done  in  real  time,  to  avoid  phase  shifts,  and  to  reduce  mean  square  error.  Mean  square 
error  from  an  optimal  non-causal  filter  is  a  little  smaller  than  from  an  optimal  causal  filter. 
A  Kalman  filter  was  not  needed  because  one  can  say  a  priori  in  what  frequency  bands  each 
sensor  has  less  noise. 

A  conceptual  block  diagram  of  the  inertial  processing  is  shown  in  Fig.  2-18.  The 
high  pass  and  low  pass  filters  are  complementary  and  I  justify  their  form  in  the  next 
section.  The  rate  gyros  are  integrated  to  make  angle  estimates  of  pitch,  roll,  and  heading. 
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Lowpass 

filter 


Fig.  2-18.  Conceptual  block  diagram  of  inertial  processing. 

A  second  estimate  of  pitch  and  roll  is  made  by  taking  the  inverse  tangent  of  the  horizontal 
acceleration  divided  by  the  vertical  acceleration.  This  estimate  has  errors  when  the  buoy 
experiences  wave  acceleration.  A  second  heading  estimate  is  made  by  taking  the  inverse 
tangent  of  measurements  of  orthogonal,  horizontal,  magnetic-flux  measurements. 

Heading  estimates  must  have  their  branch  cuts  unwrapped  before  filtering.  The  next 
section  models  the  noise  of  the  angle  estimates,  in  a  seaway. 

MODEL  OF  WAVE  MOTION  CA  USED  ANGLE  ERROR 
This  section  models  noise  in  the  computation  of  buoy  angle  to  justify  the 
complementary  filter  form  used  in  the  inertial  processing.  The  model  is  of  a  buoy  in  a 
two-dimensional  wave  field,  assumes  that  the  buoy  horizontal  acceleration  follows  the 
water  wave  horizontal  acceleration,  assumes  that  the  buoy  angle  is  computed  from 
redundant  measurements  of  a  rate  gyro  and  horizontal  and  vertical  accelerometers,  and 
linearizes  some  of  the  computations.  The  rate  gyro  has  close  to  white  noise,  and  its  noise 
is  modeled  as  white  with  the  same  variance  as  the  measured  sensor  noise.  After 
integrating  the  rate  gyro  output  to  estimate  angle,  the 
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noise  spectrum  is  given  by  equation  (2-26).  In  this  equation,  Nrr  is  the  noise  power 


_  9.5*1 0~8/W2/sec 
(/o))2  to2 


(2-26) 


spectral  density  of  the  rate  gyro  and  i|/rr  is  the  noise  spectrum  of  the  angle  estimate  from 
the  integrated  rate  gyro.  In  this  model  every  angle  and  frequency  is  expressed  in  radians 
for  convenience.  The  accelerometer-derived  angle  is  the  inverse  tangent  of  buoy-axis 
horizontal  acceleration  divided  by  buoy-axis  vertical  acceleration  (2-27).  This  equation 


t  =  tan'1-  *  -  (2-27) 

£  g 


will  be  linearized  for  the  model  leaving  the  horizontal  acceleration  divided  by  the 
acceleration  due  to  gravity.  The  model's  assumed  wave  spectrum  is  a  Pierson-Moskowitz 
spectrum  whose  energy  is  given  by  equation  (2-28)  (Newman,  1977).  In  this  equation,  a 


S(  <•>)  = 


(2-28) 


and  P  are  empirical  constants  assumed  to  be  8. 1  x  10'3  and  0.74  respectively,  g  is  the 
acceleration  due  to  gravity,  U  is  the  windspeed  at  19.5  meters,  and  o>  is  the  radian 
frequency.  The  angle  noise  spectral  density  from  this  lateral  wave  acceleration  is  then 
given  by  equation  (2-29).  Accelerometer  electronic  noise  is  modelled  as  white  with  the 


(2-29) 

measured  variance  having  a  noise  spectral  density  of  4.2  x  10‘8  rad2  s.  These  noise  spectra 
are  plotted  as  a  function  of  the  radian  frequency  for  two  windspeeds  in  Fig.  2-19.  For  the 
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Fig.  2-19.  Angle  noise  for  buoy  in  a  two-dimensional  seaway  subjected  to  the 
horizontal  accelerations  of  a  Pierson-Moskowitz  spectrum. 

higher  windspeed,  the  lateral  acceleration  spectrum  moves  up  and  to  the  left.  This  model 
of  measured  angle  noise  justifies  lowpass  filtering  the  angle  estimate  from  the 
accelerometers  and  highpass  filtering  the  angle  estimate  from  the  integrated  rate  gyro. 
Returning  to  describing  the  inertial  processing  algorithm,  the  actual  inertial  processing 
uses  a  fourth-order,  butterworth,  highpass  filter  with  a  thirty-second  period  on  the 
integrated  rate-gyro  output,  and  its  complement  on  the  other  estimates  of  Euler  angles. 

HORIZONTAL  REFERENCE  FRAME  FOR  HEADING 
As  described  in  the  Magnetic  Compasses  section,  to  get  an  accurate  measurement 
of  heading  from  magnetic  flux  requires  that  the  magnetic  flux  measurements  be  in  a 
horizontal  plane.  This  is  done  by  first  calculating  pitch  and  roll  angles,  and  then  using 
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these  angles  to  rotate  the  buoy-referenced  magnetometer  measurements  into  an 
intermediate,  horizontal,  computational  reference  frame  (2-30).  In  this  equation,  xyz 
refers  to  the  buoy  axis  and  x'y'z'  refers  to  the  intermediate  computational  reference  frame. 
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INERTIAL  ALGORITHM 

The  complete  inertial  algorithm  is  shown  as  a  block  diagram  in  Fig.  2-20.  The 
simpler,  conceptual  block  diagram  of  Fig.  2-18  linearizes  the  Euler-angle-update  matrix  as 
the  identity  matrix,  but  gives  a  reasonable  first  estimate  of  the  angles.  The  first  angle 
estimate  is  then  used  in  the  nonlinear,  Euler-angle-update  matrix  and  the  calculation  is 
repeated.  In  simulation,  the  algorithm  is  stable  and  is  more  accurate  than  linearized 
versions.  In  the  figure,  the  A^,  rotates  the  buoy-based  magnetometer  measurements  into  a 
horizontal,  computational  frame.  The  slow  0,  measurement  uses  the  inverse  tangent  of  the 
x  acceleration  multiplied  by  cost]/  divided  by  the  vertical  acceleration.  This  is  exact  and 
results  from  a  nonlinearity  in  the  3,2,1  Euler  angle  system.  The  first  and  last  minute  of 
each  processed  record  is  ignored  because  of  filter  transients.  This  algorithm  meets  the 
accuracy  needed  while  using  low-power,  noisy  ,and  drifty  sensors. 
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DATA  PROCESSING 

The  data  processing  for  the  instrument  is  represented  in  a  block  diagram  in  Fig.  2- 
21 .  The  raw  data  is  transferred  from  the  data  logger  hard  disk  to  a  personal  computer 
hard  disk  and  then  converted  into  Matlab  readable  files  with  the  programs  offld6.exe  and 


Program  or  .m  file 

offld6.exe 

bunt.exe 

vortload 


n2cal 

n2t4 

n2v 

veltr45f 

vorttran 


buoyangb 

rel2abs4 

b2E 


VEL1  VORT1  tempo 
VEL2  VORT2 


Fig.  2-21.  Vorticity  meter  data  path 


bunt.exe.  Data  is  then  loaded  into  the  Matlab  environment  and  the  two's  complement 
numbers  converted  to  integers  with  the  file  vortload.m.  In  the  Matlab  environment,  xxx.m 
files  are  sequences  of  processing  commands  that  Matlab  interprets  and  executes.  The 
processing  files  are  listed  in  appendix  B  along  with  descriptions  of  all  intermediate 
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variables.  These  integers  are  then  converted  into  physical  units  with  the  files  n2cal.m, 
n2t4.m,  and  n2v.m.  The  buoy-relative  fluid  velocities  and  vorticities  are  rotated  from 
sensor  coordinates  to  buoy  coordinates  with  the  files  veltr45f.m  and  vorttran.m.  There  are 
redundant  velocities  that  are  weighted  by  their  direction  cosines  in  computing  the  buoy- 
referenced,  relative  velocity.  Euler  angles  are  computed  by  the  file  buoyangb.m.  These 
Euler  angles  are  then  used  to  rotate  buoy  relative  velocity  and  vorticity  into  Earth 
coordinates,  adding  buoy  velocity  and  rotation  with  the  files  rel2abs4.m  and  b2E.m. 

Twice  the  buoy  rotation  rate  is  added  to  vorticity  to  compensate  for  buoy  motion,  because 
the  vorticity  of  a  nondeforming,  rotating  block  of  fluid  is  twice  the  rotation  rate.  The  final 
output  consists  of  Earth-referenced  velocities  VEL1  and  VEL2  in  cm/s,  Earth-referenced 
vorticities  VORT1  and  VORT2  in  1/s,  and  temperature  at  the  top  and  bottom  of  each 
measurement  volume  tempc  in  degrees  C.  The  numbers  1  and  2  in  the  velocity  and 
vorticity  variables  refer  to  sensors  1  and  2;  sensor  1  is  adjacent  to  the  instrument  case  and 
sensor  2  is  adjacent  to  the  float.  The  earth-based  coordinate  system  is  direction  1  east, 
direction  2  north,  and  direction  3  up.  In  general,  variables  in  lowercase  letters  are  in  buoy 
or  sensor  coordinates  and  variables  in  capitals  are  in  earth-based  coordinates  of  east,  north 
and  vertical. 


The  last  processing  stage  involves  rotating  buoy  relative  velocity  and  adding  buoy 
motion,  using  the  rel2abs4.m  file.  Velocity  in  an  inertial  reference  frame  can  be  computed 
from  velocity  in  a  moving  and  rotating  reference  frame  by  equation  (2-3 1).  The  Earth- 
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v  =  — -  +  v  +  d ixr  (2-31) 

dt  rel 

dRo 

reference  velocity  vis  equal  to  the  velocity  of  the  moving  frame  — and  the  velocity 

dt 

relative  to  the  moving  frame  f  (  ,  and  the  cross  product  of  the  rotation  rate  of  the  moving 
frame  with  the  position  vector  to  the  point.  The  buoy  velocity  is  defined  at  the 
accelerometer  block,  so  f  is  the  vector  from  the  accelerometers  to  the  center  of  the 
measurement  volume.  The  processing  of  rel2abs4.m  is  shown  as  a  block  diagram  in  Fig. 
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2-22.  In  the  top  half  of  the  figure,  the  wxr  term  is  added  to  buoy-relative  flow  in  the 
buoy  frame  and  then  rotated  into  Earth  coordinates.  In  the  bottom  half  of  the  figure,  the 


VEL1 


-g  cosOcosiJ; 


Fig.  2-22.  Algorithm  to  compute  Earth-referenced  flow  velocity 


acceleration  due  to  gravity  in  the  buoy  reference  frame  is  subtracted  from  the  buoy 
reference  accelerations.  The  net  buoy  acceleration  is  then  rotated  into  earth  coordinates, 
integrated  into  buoy  velocity,  high-pass  filtered  in  accordance  to  the  zero-average  velocity 
assumption.  This  final  buoy  velocity  is  added  to  relative  buoy  flow  velocity  to  form  the 
Earth-referenced  absolute  velocity.  As  mentioned  earlier,  this  velocity  does  not  include 
drift  currents  that  have  a  period  slower  than  thirty  seconds;  and,  the  Earth  reference  frame 
is  considered  inertial. 
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SIMULATION  OF  INERTIAL  PROCESSING 
This  section  describes  numerical  simulations  of  the  inertial  processing.  I 
conducted  simulations  to  help  develop  the  algorithm  used  and  to  estimate  the  accuracy  of 
the  algorithm  in  wave  motion.  In  these  simulations,  I  also  evaluated  errors  due  to  possible 
coning  motions,  and  the  effectiveness  in  using  an  intermediate-computational,  reference 
frame  to  calculate  heading.  The  major  causes  of  algorithm  error  are  broken  down  into  the 
following  types  of  errors:  commutivity,  integration,  round-off,  and  sensor  lag.  Lastly, 
simpler  inertial  algorithms  are  simulated  to  evaluate  the  tradeoff  between  algorithm 
complexity  and  accuracy.  ' 

Simulation  of  the  inertial  processing  is  shown  as  a  block  diagram  in  Fig.  2-23.  The 
buoy  is  subjected  to  representative  wave  motions,  sensed  physical  quantities  are 
calculated,  sensor  noise  added,  and  then  these  simulated  measurements  are  processed  by 
the  normal  buoy  software.  These  simulated  processed  variables  can  then  be  compared  to 
the  input  simulated  buoy  motion  to  measure  the  inertial  system  accuracy.  In  the 
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Fig.  2-23.  Inertial  processing  simulation  block  diagram.  Imulc,  imu2c,  and  imu3c 
are  accelerometer,  rate  gyro,  and  magnetometer,  vortb  and  watvelr  are  buoy- 
relative  vorticity  and  velocity,  VEL1,  VEL2,  VORT1,  and  VORT2  are 
computed,  earth-referenced  velocity  and  vorticity. 
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simulations,  I  used  four  deep-water  waves,  two  in  the  +X  direction  (east)  and  two  in  the 
+Y  direction  (north).  The  simulated  heading  of  the  buoy  is  the  measured  buoy  heading 
from  a  deployment  South  of  Martha's  Vineyard.  The  buoy  is  assumed  to  be  a  surface 
heave  follower  and  that  the  horizontal  velocity  of  the  top  sphere  and  bottom  instrument 
case  follow  the  horizontal  water  velocity  at  their  average  depths.  Assuming  that  the  buoy 
top  and  instrument  case  follow  the  wave  horizontal  excursions  at  their  respective  depths, 
the  contributions  to  buoy  angles  are  given  by  —  ( 1  -  eh)  multiplied  by  wave  amplitude, 
for  each  wave.  In  this  expression,  L  is  the  distance  from  the  buoy  top  to  the  instrument 
case  and  k  is  the  wavenumber.  The  simulated  measured  water  velocity  relative  to  the  buoy 
is  zero.  This  zero  relative  velocity  exaggerates  what  the  actual  vertical  velocity  would  be 
at  depth,  because  it  ignores  the  ekz  depth  attenuation  of  vertical  velocity.  The  earth- 
referenced  vorticity  is  modeled  as  zero  resulting  in  a  relative  vorticity  of  minus  two  times 
the  buoy  angle  rates.  Although  the  simulated  motion  is  arbitrary,  all  the  calculations  from 
these  assumed  motions  use  the  full  nonlinear,  exact  equations. 

In  all  except  one  of  the  simulations,  white  noise  of  the  same  variance  as  the  actual 
sensors,  was  added  to  the  simulated  sensor  measurements.  As  described  in  the  Inertial 
Measurement  Unit  Calibration  section,  the  y  axis  rate  gyro  had  a  small,  1 7-second  sinusoid 
added  to  white  noise.  In  one  set  of  simulations,  the  same  size  sinusoid  was  added  to  white 
noise  for  this  channel,  to  measure  the  inertial  algorithm’s  response  to  this  realistic  noise. 

GENERAL  ACCURACY  RESULTS  OF  SIMULA  TIONS 

The  simulations  show  that  the  inertial  measurement  unit  and  the  inertial  processing 
meet  the  specifications  described  in  the  Inertial  Measurement  Unit  section,  for  measuring 
shear  and  vorticity.  The  results  of  two  sets  of  such  wave  simulations  and  their  system 
error  are  described.  A  second  order  statistic  of  wave  height  called  significant  wave  height 
is  reviewed  and  used  to  describe  these  two  simulations.  The  significant  wave  height  is  the 
average  wave  height  of  the  one-third  largest  waves.  For  a  normally  distributed  wave  field, 
this  height  is  four  times  the  square  root  of  the  variance  of  the  sea  surface  (Earle  and 
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Bishop,  1984).  In  these  simulations,  the  waves  are  not  normally  distributed,  but  this 
formula  is  used  as  a  description  of  sea  roughness.  One  set  of  simulations  (sim7,  sim7f, 
and  sim7s)  had  four  waves  at  periods  of  nine,  six,  four  and  two  and  a  half  seconds;  had  a 
1.25-meter  significant  wave  height;  and  had  white  sensor  noise  only.  The  nine  second  and 
four  second  waves  traveled  east  while  the  six  second  and  two  and  a  half  second  waves 
traveled  north.  In  these  simulations,  the  root  mean  square  (rms)  Euler  angle  error  in  pitch 
and  roll  was  0.25  degrees  each  and  1.0  degrees  in  heading.  A  second  set  of  simulations 
(sim7c)  had  larger  waves  at  periods  thirteen,  nine,  six,  and  four  seconds;  had  a  3.7-meter 
significant  wave  height;  and  had  the  more  realistic  y-axis  rate  gyro  noise  of  white  and 
sinusoidal  noise.  The  thirteen  second  and  six  second  waves  traveled  north  while  the  nine 
second  and  four  second  waves  traveled  east.  In  the  second  set  of  simulations,  the  rms 
Euler  angle  error  in  pitch  was  0.65  degrees,  roll  was  0.9  degrees  and  heading  was  1.5 
degrees. 

The  simulations  section  makes  several  references  to  a  deployment  of  the  vorticity 
measuring  buoy  south  of  Martha’s  Vineyard  in  significant  swell.  The  wave  conditions 
during  this  deployment  are  considered  typical  for  open-ocean  conditions.  This  second  set 
of  simulations  had  the  same  rms  pitch  and  roll  angles,  and  pitch  and  roll  rates,  as  the 
deployment  south  of  Martha's  Vineyard  with  substantial  swell.  These  simulated  errors  are 
within  the  specification  error  of  three  degrees. 

INTERMEDIA  TE  COMPUTA  TIONAL  REFERENCE  FRAME  FOR  HEADING 

The  effect  on  the  accuracy  of  the  simulated  heading  with  and  without  using  an 
intermediate,  horizontal,  computational  reference  frame  was  computed.  Using  the  higher 
sea  state  of  sim7c  with  a  significant  wave  height  of  3.7  meters,  the  raw  heading  with  the 
computed  intermediate  horizontal  reference  frame  for  magnetic  flux  had  a  rms  error  of  2.2 
degrees.  After  complementary  filtering  with  the  heading  estimate  from  the  rate  gyros,  this 
error  was  reduced  to  1 .5  degrees.  If  the  raw  heading  is  computed  directly  from  the 
inverse  tangent  of  the  x  and  y  magnetic  fluxes  in  the  buoy  axes,  the  rms  heading  error  was 
15.3  degrees. 
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REJECTION  OF  CONING  MOTION 

The  ability  of  a  strapdown,  inertial  measurement  system  to  not  bias  angle  rates 
when  undergoing  a  coning  motion  is  a  classic  strapdown  system  test.  A  simulation  was 
done  that  had  two  waves  causing  coning  in  the  presence  of  two  other  waves.  As  the 
Coning  Motion  section  developed,  if  ij;  is  a  sine  wave,  0  is  a  cosine  wave,  and  the  heading 
<j)  is  zero,  then  a  vertical  rate  gyro  on  the  instrument  will  measure  a  time-average  heading 
rate  of  one  half  the  sine  and  cosine  amplitude  squared  multiplied  by  the  radian  frequency 
(2-16).  The  simulation  exposed  the  buoy  to  4.8-degree  coning  motion,  which  should  give 
a  time-average,  buoy-measured,  vertical  rotation  rate  of  0.0037  rad/s.  In  the  simulation, 
the  calculated  Euler  heading  rate  rejected  95  %  of  this  bias.  This  resulting  error  in  earth- 
referenced  vorticity  is  below  the  1  x  10'2  Is  that  is  the  measurement  goal.  Coning  motion 
of  a  buoy  built  to  measure  this  scale  vorticity  would,  in  reality,  rarely  cause  error  in  excess 
of  1  x  1  O'2  / s  even  without  compensating  for  it.  This  simulation  is,  however,  a  good 
verification  of  the  inertial  system  and  shows  that  coning  motion  is  not  a  problem  for  this 
system. 


LOW  FREQUENCY  ERRORS 

While  the  inertial  measurement  system  measures  buoy  attitude  adequately  to 
measure  vorticity  and  shear,  it  does  produce  low  frequency  errors  in  calculated  buoy 
velocity.  Errors  in  calculated  buoy  velocity  are  summed  with  measured  buoy  relative 
velocity  producing  errors  in  calculated  earth-referenced  fluid  velocities.  Most  of  the  error 
in  computed  Euler  angles  occurs  at  low  frequencies.  These  low  frequency  errors  in  pitch 
and  roll  are  on  the  order  of  a  quarter  degree.  Pitch  and  roll  angle  errors  cause  errors  in 
calculated  horizontal  acceleration  when  removing  gravity  from  the  buoy  measured 
accelerations.  These  calculated  horizontal  accelerations  are  then  integrated  to  calculate 
horizontal  buoy  and  wave  velocities.  Integrating  low  frequency  errors  amplifies  the 
errors. 


The  size  of  these  low  frequency  errors  in  velocity  can  be  inferred  from  data  taken 
from  an  instrument  deployment  in  the  open  ocean.  The  calculated  power  spectral  densities 
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of  east  and  vertical  wave  velocity  from  a  deployment  South  of  Martha's  Vineyard  are 
shown  in  Fig.  2-24.  The  horizontal  velocities  below  a  period  of  15  seconds  and  the 
vertical  velocities  below  a  period  of  25  seconds  have  large  errors.  The  vertical  velocity 
estimate  is  much  less  affected  by  the  small,  low  frequency  errors  in  pitch  and  roll  because 
in  subtracting  gravity,  vertical  g  is  multiplied  by  the  cosine  of  the  pitch  and  roll  angles. 
The  pitch  and  roll  angles  rarely  exceeded  ten  degrees.  Horizontal  g  compensation. 


Fig.  2-24.  Calculated  east  (X  vel)  and  vertical  (Z  vel)  velocity  spectral  densities  of 
the  lower  sensor  from  deployment  south  of  Martha's  Vineyard. 


however,  is  multiplied  by  the  sines  of  these  angles.  The  Error  Propagation  section 
showed  that  the  Euler-angle  error  from  the  integrated-rate  gyro,  changed  very  slowly  and 
could  have  slow  oscillations.  Pitch,  roll  wave-measuring  buoys  are  not  very  accurate 
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either  at  low  frequencies,  because  these  longer  wavelength  waves  usually  have  small 
slopes. 


In  simulating  various  test  algorithms  while  developing  the  inertial  algorithm,  it  was 
helpful  to  have  a  single  number  to  represent  low  frequency  noise.  None  of  the  simulations 
had  waves  of  frequency  lower  than  a  13  second  period.  The  statistic  chosen 
to  represent  low-frequency  noise  was  the  velocity  variance  between  the  frequencies 
corresponding  to  a  30-second  wave  and  a  15-second  wave.  This  statistic  is  just  the 
velocity  spectrum  integrated  from  1/30  s  to  1/15  s.  In  these  simulations,  any  energy  in  this 
band  is  error.  Recall  that  in  processing,  the  velocities  are  high-pass  filtered  with  a  30- 
second  high-pass  filter  so  there  should  not  be  much  energy  below  the  low  frequency 
cutoff 


Much  of  the  inertial  measurement  system  error  comes  from  sensor  noise.  The 
relative  contribution  to  error  from  sensor  noise  is  shown  in  Fig.  2-25.  This  figure  shows 
the  power  spectral  densities  of  east  velocity  in  two  simulations.  The  solid  line  is  the 
spectrum  of  the  simulated  measured  velocity  from  sim7  and  has  a  variance  of  18.5  cm2/s2 
between  30  seconds  and  15  seconds.  The  dash-dot-dash  line  is  the  spectrum  of  the 
simulation  input  velocity.  The  dotted  line  is  the  spectrum  of  a  two-dimensional, 
complementary-filtered,  horizontal-velocity  estimate  from  just  the  white  noise  from  one 
accelerometer  channel  and  one  rate  gyro  channel.  The  velocity  output  was  high-pass 
filtered  and  has  a  variance  of  3.2  cm2/s2  between  30  seconds  and  15  seconds.  A  perfect 
three  dimensional  algorithm  with  these  sensors,  at  this  sample  rate,  could  not  do  better 
than  this  level.  Most  of  the  noise  is  from  the  rate  gyro.  The  integrated  accelerometer 
noise  that  is  high-pass  filtered,  has  a  variance  of  0.2  cm2/s2  in  the  30  second  to  15  second 
period  band.  A  high-pass  filtered,  velocity  estimate  from  just  the  rate-gyro  noise  has  3.6 
cm2/s2  of  variance  in  this  band.  This  is  more  than  the  complementary-filtered,  two- 
dimensional  model  showing  some  of  the  noise  reduction  with  complementary  processing. 
The  high  spectrum  values  in  the  lowest  one  or  two  frequency  bins  is  an  artifact  of  Matlab's 
detrending  and  can  be  ignored. 
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Fig.  2-25.  Simulated  east  velocity  spectra. 


The  noise  floor  in  this  plot  above  0.4  Hz  is  primarily  due  to  the  white,  rate-gyro 
noise.  The  buoy  coordinate  system  origin  is  at  the  accelerometers  which  is  not  collocated 
with  the  measurement  volumes.  As  described  in  the  Data  Processing  section,  the  o>  x  r 
term  is  added  to  the  calculated  velocity  of  the  buoy  coordinate  origin.  This  adds  a 
substantial  radius  multiplied  by  the  rate  gyro  output,  to  buoy  velocity  and  causes  this  noise 
floor  in  velocity.  Measured  vorticity  spectra  do  not  come  down  to  the  noise  floor  due  to 
the  rate  gyro  noise.  Recall  that  twice  the  rate-gyro  output  is  added  to  measured  relative 
vorticity  in  calculating  earth-referenced  vorticity. 

Logarithmic  graphs  of  spectra  can  exaggerate  the  low  frequency  energy  of  a  signal. 
In  Fig.  2-25,  the  plot  of  power  spectral  density  of  horizontal  velocity  looks  like  error  is 
very  significant.  However,  only  4  %  of  the  velocity  variance  is  error  and  96  %  is  the  input 
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velocity. 


Getting  back  to  Fig.  2-24  and  the  low  frequency  velocity  variance;  simulations 
(sim7c)  with  the  same  rms  Euler  angles,  Euler  angle  rates,  and  the  more  accurate  y-axis 
rate  gyro  noise,  had  the  same  horizontal  velocity  variance  as  the  deployment  south  of  the 
Vineyard.  If  these  simulations  did  not  have  this  low-frequency,  velocity  variance,  the 
simulations  would  have  been  inaccurate. 

ALGORITHM  ERROR  BREAKDOWN 

The  four  major  sources  of  error  in  a  strapdown  inertial  algorithm  are  commutivity 
error,  integration  error,  round-off  error,  and  sensor  lags.  The  least-count,  analog-to- 
digital  conversion  noise  in  this  system  is  well  below  sensor  noise  and  will  not  be 
considered  further.  This  instrument  uses  a  16  bit  analog-to-digital  converter  and  the 
analog  inputs  are  scaled  to  effectively  use  the  converter's  range.  I  will  discuss  each  of 
these  four  major  error  terms  separately  in  conjunction  with  the  matrix  of  simulations  in 
Table  2-4.  The  simulations  shown  in  this  table  varied  sample  rates  while  keeping  the 
waves  constant,  and  also  doubled  the  Euler  angles  and  angle  rates  at  the  0. 15  second 
sample  rate.  For  each  simulation,  the  rms  Euler  angle  error  ( t|r,  0,  and  4>)  and  east 
velocity  variance  between  a  30-second  period  and  a  15-second  period  are  listed. 

The  commutivity  error  results  from  the  difference  between  the  actual,  smooth 
continuous  rotation  of  an  object  and  the  discrete,  sequential,  orthogonal  motion  assumed 
by  the  Euler  angle  system.  This  error  is  proportional  to  the  square  of  the  angular 
displacement  in  each  calculation  step.  The  total  error  for  an  angular  evolution  is  then 
proportional  to  the  sampling  time  step  and  is  proportional  to  the  square  of  the  angle  rates 
during  that  evolution.  For  the  Euler  angle  rates  modeled,  the  maximum  Euler  angle  error 
from  commutivity  over  a  30-second  evolution,  is  2.6  degrees.  The  simulations,  however, 
show  angle  errors  of  order  a  quarter  degree  due  to  some  canceling  out  of  these  errors  in 
cyclic  motion.  Commutivity  error  is  inherent  in  any  strapdown  inertial  system  and  can  be 
reduced  by  sampling  faster.  The  instrument  is  close  to  its  hardware  speed  limit,  so  to 
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0.41°  0.45° 


0.51  second  sampling,  sim7s 
std.  angle  error  0.36° 

1/15 

(PXJF-  124  cm2/s2 

1/30 


0.15  second  sampling,  sim7 

std.  angle  error  0.23°  0.27°  1° 

1/15 

!PXJF-  18.5  cm2 Is2 

1/30 


0.05  second  sampling,  sim7f 

std.  angle  error  0.17°  0.19°  3.1° 

1/15 

fPXJF  =  S&cmV 

1/30 

0. 1 5  second  sampling  2D 

accelerometer  and  rate  gyro  noise  into  complementary  filter 
velocity  high  passed 

1/15 

f  PXv/F  =  3.2  an2 Is2 

1/30 


double  amplitude 

0. 1 5  second  sampling,  sim7b 

std.  angle  error  1°  1°  2.25° 

1/15 

f  PXv/F  =  39 8  an  2/s  2 

1/30 


Table  2-4.  Array  of  simulations  used  to  infer  error  terms. 


sample  significantly  faster  would  require  changes  in  hardware.  Referring  to  Table  2-4, 
sim7b  is  the  same  as  sim7  but  with  all  waves,  Euler  angles  and  Euler  angle  rates  doubled 
The  angle  errors  in  the  doubled  simulation  are  quadrupled,  consistent  with  commutivity 
error  dominating  algorithm  error.  The  velocity  variance  increased  by  more  than  sixteen 
due  to  doubling  the  angle  rates.  In  sim7,  over  one-sixth  of  the  velocity  variance  is 
accounted  for  by  the  two-dimensional  sensor  noise,  sim5d. 

The  inertial  algorithm  uses  trapezoidal  rule  integration.  The  integration  error 
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during  an  evolution  should  be  proportional  to  the  time  step  squared  and  proportional  to 
the  second  derivative  of  the  quantity  integrated  (Thomas,  1972).  The  angle  error  accrued 
in  30  seconds  of  processing  with  the  simulated  angle  rates  is  of  order  0.2°.  In  cyclic 
motions,  however,  much  of  this  error  cancels  out  and  this  error  term  is  much  smaller  than 
commutivity  error.  In  one  simulation,  I  tested  second-  and  third-order  Runge  Kutta 
integration  but  the  output  error  was  not  measurably  better,  so  trapezoidal  integration  was 
kept. 


Round-off  error  in  double-precision,  floating  point  math  can  occur  if  processing  is 
not  well  conditioned.  The  complementary  filter  processing  used  is  not  well  conditioned 
because  the  crossover  frequency  is  a  very  small  fraction  of  the  Nyquist  frequency. 
Referring  to  the  simulations  in  Table  4,  when  the  sampling  rate  is  slowed  from  0.15 
seconds  to  0.51  seconds,  the  angle  error  is  increased  a  little  and  the  horizontal  velocity 
variance  is  increased  much  more.  This  increase  in  error  is  largely  due  to  commutivity 
error.  As  the  sampling  rate  is  increased  to  0.05  seconds,  however,  the  error  does  not  drop 
as  much  as  would  be  expected  if  commutivity  error  still  dominated  error  at  the  faster 
sample  rate.  Both  the  commutivity  error  and  the  integration  error  should  drop 
significantly  at  this  faster  sampling  rate.  To  explain  why  the  faster  sample  rates  do  not 
decrease  error  more,  the  complementary  filters  were  tested.  Normal  random  vectors,  of 
standard  deviation  equal  to  one,  were  complementary  filtered  with  the  filters  from  the 
three  sampling  rates.  The  slow  sampling  rate  filter  output  had  a  standard  deviation  of 
error  of  1  x  10'11.  The  medium  sampling  rate  filter  for  the  0. 15  second  sampling, 
produced  an  output  with  a  standard  deviation  of  error  of  1  x  10'10,  which  is  still  small  but 
much  larger  than  the  least  count  scale  of  2.2  *  10'16.  When  the  same  random  vector  is 
filtered  by  the  0.05  second  filter  the  standard  deviation  of  error  jumped  to  0.024.  The 
pole  zero  diagram  of  the  high-  and  low-pass  filters  for  the  0.05  second  sampling  is  shown 
in  Fig.  2-26.  When  the  poles  and  zeros  are  this  close  together,  even  double  precision  math 
can  have  errors.  If  the  hardware  allowed  it,  sampling  at  0.05  second  would  have  less  error 
than  sampling  at  0. 15  seconds;  but  if  sampling  rates  were  desired  to  be  significantly  faster 
than  0.05  seconds,  this  complementary-filter  conditioning  would  become  a  problem.  The 
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high-pass  part  of  a  complementary  filter  cannot  be  downsampled.  If  significantly  faster 
sampling  was  done,  some  other  processing  algorithm  would  have  to  be  considered. 


Fig.  2-26.  Pole-zero  plot  for  0.05  second  sample  complementary  filter 

The  electronics  do  not  have  a  simultaneous  sample  and  hold,  so  error  from  sensor 
lags  and  nonsimultaneous  sampling  have  to  be  considered.  The  velocity  and  analog 
sampling  in  each  time  step  is  done  in  under  40  milliseconds.  The  analog  inertial  sensors 
are  anti-alias  filtered  by  a  single-pole  filter  and  at  the  frequencies  of  interest,  the  filter 
phase  lag  is  close  to  a  partial-time-step  delay.  Sampling  the  analog,  anti-alias  filtered, 
inertial  sensors  after  the  velocity  measurements,  and  considering  all  samples  to  have  been 
taken  at  the  same  time,  largely  cancels  out  the  anti-alias  filter  lag.  Measurements  of  the 
phase  of  the  cross  spectrum  between  sensor-referenced  vorticity  and  its  corresponding 
rate  gyro  signal,  show  a  small  lead  for  the  analog  measurements. 
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To  estimate  the  possible  errors  resulting  from  a  partial  period  sensor  lag,  a 
simulation  was  run  with  a  full  period  lag.  The  sensor  lag  error  from  a  fractional  period  lag 
should  be  at  most  that  fraction  multiplied  by  the  sensor  lag  error  from  the  full  period  lag 
simulation.  The  simulation  that  delayed  the  analog  measurement  a  whole  time  step 
increased  errors  over  simulations  with  no  lag.  The  standard  deviation  of  pitch  and  roll 
angle  error  increased  from  0.25  degrees  to  0.65  degrees,  heading  error  increased  from  1.4 
degrees  to  1.6  degrees,  and  the  low  frequency  velocity  variance  increased  by  1  cm2/s2.  The 
maximum  actual  sensor  lag  is  a  quarter  of  a  time  step.  The  actual  sensor  lag  errors  are  at 
most  0. 1  degrees  in  pitch,  roll,  or  heading  and  0.25  cm2/s2  in  low  frequency  velocity 
variance. 

SIMPLER  INERTIAL  ALGORITHMS 

Most  oceanographic  instruments  are  sample  rate  limited.  Most  oceanographic 
instruments  are  battery  powered  requiring  slow,  low  power  electronics.  Sample  rates  are 
usually  limited  by  the  slow  electronics  or  the  storage  space  required  for  the  months  of  data 
that  a  typical  oceanographic  instrument  is  deployed.  In  strapdown  inertial  processing, 
faster  sample  speed  for  the  same  algorithm  complexity  gives  better  accuracy.  This 
instrument  development  was  hardware  limited  in  speed  and  sample  rate  requiring,  by 
oceanographic  standards,  a  complex  inertial  algorithm.  In  a  more  general  instrument 
development,  the  designer  has  to  tradeoff  algorithm  complexity  and  sample  rate. 

Simpler  inertial  algorithms  were  simulated  to  evaluate  the  tradeoff  between 
algorithm  complexity  and  accuracy.  Three  of  the  simulations  are  summarized  in  Table  2- 
5.  The  first  simulation  listed  is  the  same  sim7  as  in  Table  2-4  for  comparison.  The  second 
simulation  listed,  sim71,  used  a  linear,  Euler-angle  update  matrix  and  the  same 
complementary  filtering  as  sim7.  The  standard  deviations  of  angle  error  increased  for 
pitch  and  roll  from  0.23°  and  0.27°  to  0.40°  and  0.47°  while  the  heading  error  standard 
deviation  stayed  the  same  at  one  degree.  The  low  frequency  velocity  variance,  however, 
jumped  from  18.5  cm2/s2to  63  cm2/s2.  If  only  buoy  attitude  was  desired  from  a  strapdown 
inertial  system,  this  simpler  algorithm  might  be  adequate. 
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A  simulation  was  run,  but  not  listed,  that  took  out  the  cosi|r  term  in  computing  the 

slow  estimate  of  0.  For  a  3,2, 1  Euler  angle  system,  the  slow  estimate  of  0  is  given  by 

tan1  *  cosl^ .  In  these  simulations,  the  pitch  and  roll  angles  never  exceeded  15°  so  the 
z 

simulation  results  were  the  same  with  and  without  the  cosijf  term.  This  term  was  left  in 
because  theoretically  it  should  be  there  and  if  the  pitch  and  roll  angles  were  greater,  should 
improve  accuracy. 


Full  nonlinear  processing  sim7 

std  angle  error  0.23°  0.27°  1.0° 

1/15 

f  PXv/F  =  18.5  an  2 Is 2 

1/30 

Linear  Euler-angle  update  matrix  sim71 

std  angle  error  0.40°  0.47°  1.0° 

1/15 

f  ^Xve^  =  63  cm2 Is2 

1/30 

Linear  update  matix,  first-order,  causal  complementary  filter 

std  angle  error  0.54°  0.59°  8.04° 

1/15 

f  Pa/F  -  8°  cmV 

1/30 

Table  2-5.  Simulations  of  simpler  inertial  algorithms 

The  last  simulation  listed  in  Table  2-5  is  of  an  algorithm  with  a  linear  Euler  angle 
update  algorithm  and  first-order  causal  high-  and  low-pass  complementary  filters.  Its 
pitch  and  roll  error  standard  deviations  increased  to  0.54°  and  0.59°  while  its  heading  error 
jumped  to  over  8°.  This  simpler  algorithm  gives  a  good  vertical-velocity  spectrum  but  has 
unacceptable  heading  errors.  If  in  a  future  instrument  development,  in-situ  processing  of 
attitude  was  desired,  one  of  these  simpler  algorithms  sampled  and  processed  at  a  faster 
rate  would  have  to  be  considered. 
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These  inertial  simulations  show  that  this  inertial  system  meets  the  specifications  for 
measuring  shear  and  vorticity.  The  simulations  of  estimated  error  in  measured  vorticity 
due  to  inertial  system  errors,  predict  mean  errors  of  5  *  10'5  /s  and  standard  deviation  of 
error  of  4  x  10'3  /s.  Most  of  the  time  varying  error  is  due  to  rate  gyro  noise.  The 
measured  standard  deviation  of  vorticity,  in  an  open  ocean  deployment,  was  fifteen  times 
this  magnitude  showing  that  most  of  the  measured  vorticity  variance  is  due  to  wakes  and 
turbulence  in  the  flow  and  not  IMU  errors.  If  a  system  is  desired  to  measure  horizontal 
velocity  spectra  to  periods  longer  than  15  seconds,  faster  sampling  and/or  better  rate  gyros 
are  needed. 


\ 
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CHAPTER  3.  SENSOR 
PERFORMANCE 


This  chapter  describes  the  vorticity  sensor’s  major  sources  of  error  and  also 
describes  measurements  of  the  instrument's  accuracy.  The  inertial  system  simulations  of 
the  last  chapter  showed  that  errors  in  instantaneous  and  mean  vorticity  resulting  from 
inertial  system  errors  was  much  less  than  1  *  10‘2  per  second.  Errors  from  wakes  of  the 
sensor  advecting  through  the  sensing  volume  are  the  major  error  sources.  This  chapter 
models  in  constant  flow:  potential  flow  disturbance  due  to  the  transducer  pods,  wake  error 
due  to  the  center  stalk  wake,  wake  error  due  to  the  transducer  pods,  circulation 
disturbance  from  lift  on  the  center  stalk,  and  sound  speed  changes.  This  chapter  then 
shows  through  measurements  that:  electronic  noise  and  bias  are  much  less  than  1  x  1  O'2 
per  second,  symmetry  of  the  sensor  design  succeeded  in  keeping  biases  in  constant  flow 
to  less  than  1  x  10'2  per  second,  down-wave  and  vertical  vorticity  biases  in  open-ocean 
waves  are  less  than  1  *  10'2  per  second,  and  cross-wave  vorticity  bias  in  open-ocean 
waves  is  about  1  *  10'2  per  second  and  scales  with  rms,  sensor-relative  velocity. 

3.A.  MODELING  SENSOR  ERRORS 

This  section  develops  models  of  sensor  flow  disturbance  and  error.  Potential  flow 
is  used  to  model  flow  around  the  transducer  pods  and  predicts  a  reduced  gain  (that  can  be 
calibrated)  but  no  bias  or  noise.  Wake  models  are  developed  to  describe  bias  in  measured 
circulation  due  to  the  center  stalk  and  the  transducer  pods.  The  circulation  bias  that 
would  result  from  lift  off  the  center  stalk  is  predicted  and  the  result  is  used  to  justify  not 
streamlining  the  cylinders  in  the  sensor. 

The  potential  flow  solution  of  flow  around  spheres  can  be  used  to  gain  some 
insight  into  flow  errors.  The  geometry  and  coordinate  system  modeled  are  shown  in  Fig. 
3-1.  This  model  only  considers  the  spheres;  it  ignores  pod  arms,  and  it  assumes  that  the 
spheres  are  moving  relative  to  a  still  fluid.  In  this  model,  L  is  the  distance  between 
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Fig.  3-1.  Geometry  of  potential  flow  model. 

spheres,  0  is  the  angle  between  the  flow  and  the  acoustic  path,  r  is  the  radius  from  the 
sphere  center,  and  rQ  is  the  sphere  radius.  The  sphere  diameter  is  much  smaller  than  the 
distance  between  the  spheres  so  there  is  little  interaction  between  the  sphere  disturbance 
fields.  This  model  calculates  the  line  integral  of  velocity  disturbance  from  one  sphere  over 
an  acoustic  length  L,  and  doubles  it.  The  flow  potential  field  O  for  the  disturbance  field 
from  one  sphere  is  given  by  equation  (3-1).  The  radial  velocity  line  integral  for  the 


2 


COS0 


(3-1) 


disturbance  from  one  sphere  is  then  equation  (3-2).  The  potential  flow  disturbance  retards 
radial  flow  both  upstream  and  downstream  of  the  sphere  causing  acoustic  velocity 
measurement  to  under-report  velocity.  To  calculate  the  velocity  speed  error,  the  radial 
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cos0  dr  = 


Ur  cos0 
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velocity  line  integral  is  doubled  and  divided  by  the  line  integral  of  undisturbed  velocity  (3- 
3).  For  the  1 5 -centimeter  path  vorticity  meter,  the  velocity  error  is  0.0741  cos  (0)  *  U  or 


velocity  error  d  _ 

- ^ -  =  — COS0 

velocity  2  L 


(3-3) 


about  7.5  %  of  the  flow,  which  is  a  significant  error.  However,  because  of  the  cos  0  term, 
this  error  just  changes  the  velocity-calibration  constant  of  the  sensor.  Potential  flow 
theory  predicts  that  if  the  flow  upstream  of  the  sensor  is  irrotational,  the  flow  around  the 
vorticity  sensor  will  stay  irrotational.  If  the  flow  around  the  sensor  is  irrotational,  there 
will  be  zero  measured  circulation,  ie,  no  vorticity  bias  in  irrotational  flow.  If  the  flow 
upstream  of  the  sensor  has  a  small  amount  of  shear  compared  to  the  shears  of  the  sphere’s 
disturbance  fields,  the  above  model  predicts  under-reporting  of  the  shear  or  vorticity  by 
the  same  ratio  that  velocity  is  under-reported.  Calibrating  the  velocity  gain  should  also 
calibrate  the  vorticity  gain.  This  potential  flow  model  should  apply  well  upstream  of  a 
sphere,  but  because  of  boundary  layer  separation,  the  model  is  not  realistic  downstream  of 
a  sphere. 


Models  were  developed  for  the  circulation  error  resulting  from  wakes  from  the 
center  stalk  and  transducer  pods.  Downstream  of  bluff  bodies,  wakes  resulting  from 
boundary  layer  separation  over  bluff  bodies  disturb  the  flow  more  than  the  potential  flow 
solution.  Fig.  3-2  represents  the  cross  sections  of  two  vorticity  meters,  showing  the 
center  stalk  in  the  middle  of  a  square  of  transducer  pods.  The  top  half  of  the  figure  shows 
a  circulation  square  path  in  the  worst  direction  for  velocity-path-wake  errors.  As  explained 
in  the  Mechanical  Design  Overview  section,  the  sensor  was  designed  with  symmetry  to 
minimize  error  when  measuring  circulation  or  vorticity.  The  lower  half  of  the  figure 
shows  the  same  sensor  at  an  unsymmetrical  angle  of  attack.  The  wake  errors  on  paths  2 
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Fig.  3-2.  Geometry  of  wake  error  models. 

and  4  are  the  largest  and  still  cancel  in  measuring  circulation.  The  wake  errors  in  paths  1 
and  3,  however,  do  not  cancel.  Two  models  were  developed  to  estimate  circulation  error 
resulting  from  path  wakes  from  the  center  tube  and  from  a  sphere. 


The  first  model  developed  was  for  the  two-dimensional  wake  of  the  center  stalk. 
The  model  assumes  the  far-field,  two-dimensional  turbulent  wake  described  in  Schlichting, 
(1979).  The  model  geometry  is  shown  in  Fig.  3-3  and  assumes  an  entrainment  relation 
given  by  equation  (3-4).  In  this  equation,  b  is  the  wake  width,  x  is  the  distance 


(3-4) 


downstream,  P  is  the  entrainment  constant,  U„  is  the  free  stream  velocity,  and  u,  is  the 
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Fig.  3-3.  Geometry  of  two-  and  three-dimensional  models  of  wake  circulation 
disturbance.  The  two-dimensional  model  stalk  diameter  is  D  and  the  three- 
dimensional  model  pod  diameter  is  d. 

wake  velocity  difference  from  the  undisturbed  fluid.  The  momentum  deficit  of  the  wake 
is  given  by  equation  (3-5)  and  is  constant  with  distance  downstream.  In  this  equation,  n  is 


/  =  p  fu(Um-u)dA 
J  =  -p  U^u/A 


(3-5) 


the  velocity  in  the  wake,  J  is  the  drag  per  unit  length  on  the  center  stalk,  p  is  density,  and 
the  area  over  which  the  wake  is  integrated  is  perpendicular  to  the  wake  and  downstream 
of  the  strut.  The  wake  width  scales  as  equation  (3-6)  and  the  wake  velocity  scales  as 


b  ~ 


(3-6) 


equation  (3-7). 
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The  circulation  disturbance  through  the  wake  is  the  line  integral  of  velocity  along  a 
circulation  path  and  is  given  by  equations  (3-8).  Only  the  line  integral  through  the  wake 


a  circulation  ~ 


/"• 

h 

j"  i/jtanOiy 


dl 

sin  Qdl 


(3-8) 


contributes  to  circulation  disturbance  because  outside  the  wake,  even  in  the  areas  where 
the  flow  is  accelerated,  the  flow  is  irrotational.  Using  the  scaling  expressions  in  equation 
(3-8)  gives  equation  (3-9)  where  K  is  an  as  yet  undetermined  constant.  The  far  field 

ai circulation  =  Ku  {b  (3-9) 


momentum  deficit  and  circulation  disturbance  equations  are  both  linear  in  u,  the  velocity 
anomaly  and  b  the  wake  width,  which  means  that  the  actual  distribution  of  u,  is 
unimportant.  For  this  model  a  top  hat  distribution  is  assumed  where  u,  is  constant  over 
the  wake  width.  The  undetermined  constant,  K,  can  be  determined  from  the  momentum 
equation  (3-5)  to  be  one  half.  The  result  for  the  circulation  and  vorticity  disturbance  is 
given  by  equation  (3-10).  Line  averaging  across  a  two-dimensional  wake  in  the  far  field. 


a  circulation  =  —U  C^D  tan0 

2  ~  D 


Avorticity 


J  UCdD  tan0 


(3-10) 


is  constant  with  distance  downstream,  is  independent  of  the  velocity  distribution  in  the 
wake,  and  is  independent  of  p  the  entrainment  constant.  The  measured  circulation 
depends  on  the  circulation  path  in  rotational  flows. 


This  bias  correction  from  the  center-stalk  wake  model  was  applied  to  data  taken 
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with  the  sensor  towed  through  still  water.  The  geometry  of  this  model  application  is 
shown  in  Fig.  3-4.  The  circulation  path  squares  look  like  the  figure  when  viewed  parallel 


u 


Fig.  3-4.  Center  stalk  model  used  on  tow  tank  data 


to  the  center  stalk.  Half  of  the  wake  is  assumed  to  flow  on  each  side  of  the  downstream 
pod  and  the  resulting  measured  wake  circulation  and  vorticity  errors  are  calculated  by 
equation  (3-1 1).  Vorticity  vectors  perpendicular  to  each  circulation  path  square  are 


circulation  error  =  U  — tan(60")  -  U  — 130(30°) 

“4  ”4 

0.288  UD  (3’u) 

vorticity  error  =  - - — 

L 2 


labeled  avort ,  bvort,  and  cvort.  The  vorticity  measurement  avort  was  measured  with  its 
circulation  path  as  shown  in  the  Fig.  3-4  and  cvort  was  measured  with  a  circulation  path 
that  was  the  mirror  image  of  the  figure.  In  the  tow  tests,  the  vorticity  was  rotated  into 
axes  parallel  to  the  tank  yvort,  across  the  tank  xvort,  and  vertical  zvort.  The  error  in  these 
coordinates  was  then  given  by  equation  (3-12).  This  is  a  simple,  far-field  model 
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xvort  = 
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xcvort 


—  xavort  -  — 

s[2  fl 

0.288  (7  D 

xvort  =  y2  x - - — 


(3-12) 


being  applied  in  the  near  field.  This  simple  model’s  assumption  of  the  wake  flowing  on 
each  side  of  the  pod  and  the  pod  not  effecting  the  measured  circulation  is  unrealistic. 

Some  of  the  stalk  wake  flows  around  the  downstream  pod  in  the  dimension  parallel  to  the 
stalk  and  out  of  the  plane  of  this  two-dimensional  model.  The  actual  correction  applied  to 
measurements  from  the  45 -centimeter  path  vorticity  sensor  was  68.7%  of  the  above 
model,  and  in  constant  flow  accounted  for  89%  of  the  measured  bias. 


The  second  model  developed  was  for  the  three-dimensional  wake  of  the  transducer 
pods.  The  geometry  of  the  model  is  shown  by  Fig.  3-3  again,  using  d  for  the  pod 
diameter.  The  entrainment  model  used  was  equation  (3-4)  as  in  the  two-dimensional  wake 
model.  The  width  is  proportional  to  the  distance  downstream  to  the  one-third  (3-13)  and 
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the  wake  relative  velocity  is  proportional  to  the  distance  downstream  to  the  minus  two- 
thirds.  The  measured  wake  circulation  error  was  then  given  by  equation  (3-14)  for  an 


a  circulation  ~  U 

O0 


p3*3 


tan0 


(3-14) 
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acoustic  path  through  the  wake  center.  The  measured  vorticity  error  assuming  the 
distance  downstream  is  L,  is  given  by  equation  (3-15).  This  result  shows  that  the 


a vorticity  ~ 


U(ct 


4 

3  tan0 


(3-15) 


measured  circulation  or  vorticity  error  changes  with  distance  downstream,  depends  on  the 
entrainment  constant  P,  and  depends  on  the  velocity  distribution  in  the  wake.  Although 
applying  this  far-field  wake  model  to  the  relatively  near-field  wake  in  this  sensor 
application  is  tenuous,  this  solution's  scaling  of  error  with  pod  size  and  path  length 
provides  a  framework  in  which  measurements  of  error  can  be  compared. 


The  circulation  bias  that  would  result  from  lift  from  the  center  stalk  will  now  be 
calculated  to  justify  not  streamlining  the  stalk.  If  the  center  stalk  were  not  round,  there 
could  have  been  lift  off  the  stalk  and  consequent  lift-related  circulation.  The  geometry  of 
a  lift-related  circulation  model  is  shown  in  Fig.  3-5.  The  relation  between  lift  and 


Fig.  3-5.  Geometry  of  lift  related  circulation  model. 
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circulation  is  given  by  equation  (3-16)  where  CL  is  the  coefficient  of  lift,  T  is  circulation. 


T  =  udx. 

L  =  pLT  =  ^pu2dcl 


(3-16) 


and  p  is  density.  The  change  in  circulation  and  vorticity  is  given  by  equations  (3-17).  As 


a circulation 
Avorticity  = 


=  0.5UDCl 
0 .5UDCl 

~L2 


(3-17) 


an  example,  with  the  dimensions  of  the  45-centimeter  path  vorticity  sensor  if  the 
coefficient  of  lift  was  one  (a  reasonable  lift  coefficient  for  a  streamlined  strut)  and  the 
constant  current  was  one  meter  per  second,  the  lift-related  vorticity  would  be  0.086  per 
second.  If  the  center  stalk  was  streamlined  and  had  an  angle  of  attack  relative  to  the  flow, 
lift  caused  by  asymmetrical  boundary  layer  separation  off  the  stalk  would  cause  significant 
bias.  Rather  than  trying  to  deal  with  this  bias,  the  center  stalk  was  kept  round  and  as  small 
as  structurally  sound,  and  should  have  zero  time-average  lift. 


3.B.  MEASUREMENTS  OF  PERFORMANCE 

This  section  measures  electronic  noise  and  bias  in  vorticity  to  show  that  they  are 
much  less  than  1  x  10'2  /s.  Sound  speed  changes  in  gain  are  predicted.  Bias  and  noise  in 
constant  flow  are  measured  to  be  less  than  1  *  10'2  Is.  Biases  in  simulated  wave  flow  are 
measured  to  show  that  down- wave  and  vertical  vorticity  biases  are  less  than  1  *  1 0"2  / s 
and  that  cross-wave  bias  scales  with  rms  relative  velocity.  Measurements  in  a  wave  tank 
and  in  the  ocean  are  used  to  show  that  vorticity-meter  bias  is  insensitive  to  heading.  A 
deployment  with  significant  swell  and  negligible  wind  stress  will  be  described  to  show  that 
cross-wave  vorticity  bias  in  open-ocean  waves  is  about  1  x  10'2  /s.  Vorticity  meter 
rejection  of  waves  is  measured  to  show  how  effective  the  vorticity  measurement  is  at 
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removing  irrotational  motion  from  its  measurements. 


ELECTRONIC  NOISE 

Flow  distortion  due  to  the  instrument  being  in  the  water  dominates  the  electronic 
instrument  noise  and  bias  in  all  but  the  most  still  flows.  The  two  significant  sources  of 
electrical  error  are  electronic  noise  and  zero  drift.  Current  measurements  made  with  a  1 5- 
centimeter  path  sensor  are  shown  in  Fig.  3-6.  The  top  graph  is  measured  velocity  in  a  still 


Fig.  3-6.  Current  measurement  in  a  bucket.  The  top  graph  shows  typical  bucket 
convective  water  velocities  summed  with  electronic  noise  and  the  bottom  graph 
plots  the  measured  velocities  in  a  gel,  that  is  still,  leaving  electronic  noise. 
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bucket  of  water  and  the  lower  graph  is  measured  velocity  in  a  gel  of  carrageenin,  to  stop 
convective  flow.  The  electronic  noise  of  this  measurement  has  a  standard  deviation  of 
0.0095  cm/s.  Normal  convective  velocities  in  a  still  bucket  are  typically  2  to  5  mm/s. 

ELECTRONIC  ZERO  DRIFT 

The  zero  of  this  instrument  can  change,  by  as  much  as  0.3  cm/s,  if  the  sensor 
cables  are  moved  or  reconnected,  causing  a  bias.  To  make  accurate  measurements  with 
this  instrument  requires  calibrating  the  zero  flow  in  a  bucket  of  still  water  each  time  the 
sensors  are  reconnected.  These  “bucket  zeros”  have  to  be  performed  carefully,  avoiding 
differential  bucket  heating  and  averaged  for  many  hours.  Typical  “bucket  zeros”  of 
velocity  before  and  after  a  deployment,  vary  by  less  than  a  millimeter  per  second  on  each 
acoustic  path.  Because  the  cables  are  under  atmospheric  pressure  all  the  time,  the 
instrument  zero  does  not  change  with  instrument  depth.  If  the  cables  were  exposed  to 
hydrostatic  pressure,  the  resulting  compression  would  change  the  cable  capacitance  and 
could  change  the  instrument's  zero.  To  put  these  zero  changes  in  perspective,  typical 
propeller  current  meters  do  not  turn  until  flows  reach  one  centimeter  per  second.  Except 
in  very  still  conditions  electrical  instrument  errors  can  be  ignored. 

SOUND  SPEED  CHANGES 

If  the  sound  speed  in  water  changes  due  to  changes  in  temperature  or  salinity,  and 
if  this  local  speed  of  sound  in  water  is  not  used  to  compute  the  water  velocity  from  the 
differential  sound  travel  time  (2-2),  an  error  in  sensor  gain  will  result.  This  gain  error  will 
be  (c2-c2)/c2  where  ca  is  the  assumed  speed  of  sound  and  cr  is  the  actual  speed  of  sound. 
A  one  degree  centigrade  difference  in  temperature  can  result  in  a  0.6  %  error  in  gain.  This 
potential  source  of  error  can  be  avoided  by  using  the  local  sound  speed. 

THERMISTOR  TIME  RESPONSE 

The  thermistor  time  constant  was  measured  by  plunging  the  instrument  into  cool 
water  from  warm  air  and  measuring  the  time  delay  for  the  measured  temperature  to  reach 
63%  of  the  final  temperature  change.  This  time  constant  model  assumes  a  first  order 
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thermal  system.  The  time  constant  was  effected  by  the  relative  flow  velocities  of  the 
water,  and  varied  from  12  seconds  to  25  seconds.  The  thermal  response  time  is  long 
because  the  thermistors  were  mounted  inside  the  fairly  thick  wall,  stainless  steel  tubing. 
The  purpose  of  the  thermistors  was  to  measure  time  average  thermal  stratification.  The 
thermistor  time  constant  was  sacrificed  for  protection. 

BIAS  IN  CONSTANT  FLOW 

WAKE  ERROR  CANCELLA  TION 

Wake  error  cancellation  was  measured  to  evaluate  the  effectiveness  of  the  sensor 
design,  in  reducing  bias  in  constant  flow.  Wake  error  cancellation  with  a  square 
circulation  path  is  shown  in  Fig.  3-7.  The  15-centimeter  path  sensor  was  towed  with  two 
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Fig.  3-7.  Measured  path  velocities  around  circulation-path  square  while  sensor  was 

towed  through  still  water. 
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of  its  acoustic  paths  parallel  to  the  carriage  velocity,  the  worst  orientation  for  measuring 
velocity.  The  acoustic  paths  parallel  to  the  flow  measured  about  10%  less  than  the  free 
stream  velocity,  but  these  two  decrements  were  the  same  to  within  1%  of  the  carriage 
velocity.  The  symmetry  of  the  sensor  cancels  90%  of  the  error  from  transducer  pod 
wakes.  Lines  were  drawn  at  the  carriage  velocity  and  show  initial  velocity  measurements 
with  no  or  reduced  wake  decrement.  Until  the  sensor  has  moved  a  path  length,  there  is 
still  undisturbed  fluid  in  at  least  part  of  the  acoustic  path.  In  a  wave  field,  water  is 
constantly  accelerating  in  different  directions  so  constant-flow  sensor  calibrations  cannot 
automatically  be  used. 

The  three-axis  vorticity  sensor  also  has  acoustic  paths  that  form  triangles  so 
vorticity  bias  in  constant  flow  for  a  triangular  path  can  be  computed.  For  a  tow  at  0.37 
m/s,  the  vorticity  bias  for  triangle  circulation  paths  on  the  top  and  bottom  of  the  sensor  is 
0.53  /s,  while  the  bias  for  the  square  path  is  0.027  /s,  Fig.  3-8.  The  vorticity  bias  of  the 
triangle  was  twenty  times  that  of  the  square  because  of  the  90%  error  cancellation  of  the 
square  and  because  the  area  of  the  triangle  was  half  the  area  of  the  square  path.  While  this 
sensor  greatly  reduces  vorticity  bias  in  constant  flow  by  designing  the  acoustic  circulation 
paths  to  be  symmetric  with  sensor  wakes,  it  does  not  reduce  the  intensity  or  size  of  the 
sensor  wakes.  In  a  wave  field,  these  strong  wakes  left  by  the  sensor  can  advect  back  into 
the  measurement  volume;  this  resulting  error  will  be  treated  separately  in  the  bias  in  wave 
flow  section. 
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Fig.  3-8.  Vorticity  bias  in  constant  flow  for  triangular  and  square  (avort) 
circulation  paths.  Tow  speed  0.37  m/s. 

MEASURED  BIAS  IN  CONSTANT  FLOW 

Vorticity  bias  in  constant  flow  was  measured  by  towing  both  15 -centimeter  path 
and  45-centimeter  path  vorticity  sensors  in  tanks  of  still  water.  Uncorrected  average 
vorticities  of  both  sensors  for  different  tow  carriage  speeds  are  shown  in  Fig.  3-9.  For  the 
45-centimeter  path  sensor,  the  x  direction  was  horizontal  and  perpendicular  to  the  flow, 
the  y  direction  was  parallel  to  the  flow,  and  the  z  direction  was  vertical.  These 
measurements  were  corrected  for  the  center  tube  bias  error  discussed  in  the  Modeling 
Sensor  Errors  section,  Fig.  3-10.  A  separate  plot  of  the  measured  biases  of  the  15- 
centimeter  path  vorticity  meter  with  its  axes  labeled  is  shown  in  Appendix  C.  Both 
sensors  have  the  same  ratio  of  center  stalk  diameter  to  path  length;  the  buoy  center  stalk 
had  to  be  this  large  to  be  strong  enough  to  be  picked  up.  The  pod  diameter  to  path  length 
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;ig.  3-9  Uncorrected  vorticity  bias  of  15-  and  45-centimer 
:onstant  flow. 


Fig.  3-10.  Vorticity  bias  of  15-  and  45-centimeter  sensor  in 


stalk-bias  correction  described  in  text. 


is  62.5%  smaller  for  the  larger  45  cm  path  sensor.  Correcting  for  center  tube  bias  reduced 
bias  to  below  0.01  /s  in  constant  flow  up  to  70  cm/s.  Above  70  cm/s,  the  pod  arms 
strummed  in  the  parallel,  double  frequency  mode  and  measurements  of  velocity  noise 
drastically  increased.  The  x-axis  (cross  channel)  bias  error  at  0.257  m/s  corresponds  to 
the  wave  speed  of  a  wavelength  just  longer  than  the  diameter  of  the  center  tube.  The 
depth-varying,  extra  drag  at  this  speed  is  a  Froude-number,  free-surface  effect  of  the  test, 
resulting  from  the  center  stalk  penetrating  the  water  surface.  In  an  actual  deployment,  all 
of  the  center  stalk  is  underwater  and  this  Froude  effect  should  not  exist  at  that  velocity. 


Correcting  for  the  center  stalk  bias  leaves  mostly  pod  wake  errors  which  scale  as 

the  pod  diameter  over  the  path  length.  The  biases  corrected  for  the  center  stalk  wake  for 

both  size  sensors  for  a  0.643  m/s  tow  are  listed  in  Table  3-1.  The  ratio  of  these  biases  is 

about  nine.  The  t|iree-dimensional  pod-wake  model  predicted  that  vorticity  bias  would 

U(  d\ 


scale  as  —  ^  —  j  .  Using  the  dimensions  of  the  15-centimeter  path  and  45-centimeter 
path  vorticity  meters  in  the  pod- wake  model ,  the  model  predicts  a  ratio  of  bias  of  nine. 


Corrected  bias 


X 

y 

z 

15 -cm  path  sensor 

-0.068 

0.057 

-0.064 

45-cm  path  sensor 

0.0067 

0.007 

-0.006 

ratio 

10.2 

8.2 

10.7 

Table  3-1.  Vorticity  sensor  bias  in  0.643  m/s  constant  flow  corrected  for  center 
stalk  bias. 

Comparing  error  data  from  the  tow  tank  and  strumming  onset  speed,  shows  one  of 
the  design  tradeoffs.  The  45-centimeter  path  sensor  was  designed  with  small  cross 
sections  to  reduce  wake  effects  which  reduced  error,  but  this  also  lowered  the  strumming 
onset  speed.  For  any  sensor  size,  there  is  a  tradeoff  between  the  error  within  an  operating 
envelope  and  the  size  of  that  operating  envelope. 
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VARIANCE  IN  CONSTANT  FLOW 

The  size  of  velocity  fluctuations  caused  by  vortex  shedding  was  measured.  The 
geometry  is  shown  in  Fig.  3-11  where  D  is  the  center  stalk  diameter  and  L  is  the  acoustic 
path  length.  The  expected  fluctuating  velocity  is  given  by  equation  (3-18). 

Ml  DU„ 

a velocity  =  — -  =  constant  * - -  (3-18) 

l 

The  velocity  of  the  vortices  shed  scale  with  the  free-stream  velocity  and  the  vortice’s  size 
should  scale  with  the  stalk  diameter.  Acoustic  measurement  of  velocity  line  averages  the 
fluctuating  velocity  over  the  acoustic  path.  The  constant  in  this  parameterization  was 
measured  by  towing  the  15-centimeter  path  sensor  through  still  water  Fig.  3-12.  The 
measured  constant  varied  from  1.2  to  1.8. 


Fig.  3-1 1.  Geometry  used  to  measure  size  of  shed-vortex,  measured  velocity 
variance. 
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Fig.  3-12.  Fluctuating  measured  velocity  perpendicular  to  and  behind  center  stalk 
of  15  cm  path  vorticity  meter  in  4.06  cm/s.  tow. 

Fluctuations  of  velocity  perpendicular  to  the  center  stalk  and  tow  velocity,  were 
measured  to  increase  with  tow  velocity.  A  power  spectral  density  of  this  measured 
velocity  at  one  tow  speed  is  shown  in  Fig.  3-13.  The  large  spectral  peak  in  the  center  of 
this  figure  corresponds  to  the  shedding  frequency  of  the  center  stalk.  Spectra  of  cross 
velocity,  such  as  shown  in  Fig.  3-13,  were  measured  at  different  tow  speeds  and  are 
compiled  together  in  Fig.  3-14.  As  the  carriage  moved  faster,  the  peak  frequency 
increased  and  velocity  variation  amplitude  increased.  The  peak  at  the  highest  speed  was 
aliased  and  was  actually  above  the  Nyquist  frequency.  At  the  carriage  speeds  run,  vortex 
shedding  variance  from  the  pod  arms  was  aliased. 
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Fig.  3-13.  Power  spectral  density  of  measured  cross  tank  X  and  parallel  Y 
velocities  of  a  45-cm  path  sensor  towed  at  0.26  m/s  showing  the  Strouhal  peak  in 
cross  velocity. 


carriage  speed  in  m/sec 


Fig.  3-14.  Power  spectral  density  of  measured  cross  tank  velocity  at  different 
carriage  speeds. 
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BIAS  IN  WAVE  FLOW 

Measurement  bias  of  mean  vorticity  in  a  wave  flow  is  difficult  to  accurately 
measure  or  model.  Vorticity  errors  in  a  wave  flow  are  primarily  from  the  sensor's  wake 
advecting  back  into  the  measurement  volume.  Subtle  differences  in  how  and  where  the 
wake  is  advected  back  could  significantly  affect  errors.  Constant  flow  models  where  the 
sensor  wake  is  only  advected  away  from  the  measurement  volume  cannot  be  directly 
applied  to  a  wave  flow. 

A  nondimensionalization  of  error  in  a  monochromatic  wave  field  is  equation  (3- 
19).  In  this  equation,  is  the  vorticity  bias,  L  is  the  path  length,  u  ’  is  the  rms  relative 


wbiasL  r/  Tu'  gL 

-  =  t  (geometry,  - ,  — — , 

u '  L  u'2 


h  z 
- 7  t),  - , 


(3-19) 


u 


u 


velocity,  T  is  the  wave  period,  g  is  the  acceleration  due  to  gravity,  h  is  the  water  depth,  0 
is  the  heading  angle  of  the  sensor  with  respect  to  the  waves,  z  is  the  sensor  depth,  u  is  the 
average  flow  past  the  sensor  parallel  to  the  waves,  and  v  is  the  average  flow  past  the 
sensor  perpendicular  to  the  waves.  The  Tu  7L  term  is  the  square  root  of  two  multiplied  by 
the  Keulegan-Carpenter  number  which  is  proportional  to  the  wave  excursion  divided  by 
the  sensor  diameter  (Faltinsen,  1990).  The  gL/u  ’2  term  is  an  inverse  Froude  number  for  the 
relative  wave  velocity.  This  nondimensionalization  is  convenient  to  apply  to  laboratory 
rotating  arm  tests  of  bias,  wave  tank  tests  of  bias,  and  is  also  convenient  for  estimating  the 
accuracy  of  a  measurement  after  it  has  been  made.  This  formulation  is  not,  however,  the 
most  convenient  to  apply  in  a  predictive  role  given  an  expected  sea-state  without  knowing 
the  buoy  response  function.  The  measured  buoy  response  functions  are  plotted  in 
Appendix  C. 


In  the  tests  performed  so  far  gL/u  ’2, 6,  and  mean  drift  u  were  varied.  To  keep 
geometry  constant,  only  data  taken  with  the  45-centimeter  path  sensor  is  discussed  here. 
Because  of  the  sensor  size,  the  water  depth  and  sensor  depth  could  not  be  significantly 
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varied  in  laboratory  tests.  Cross  wave  velocity  drift,  v ,  was  zero  for  all  the  data  taken  so 
its  effect  cannot  be  estimated  from  the  data.  The  parameters  that  I  varied  were  wave 
period,  wave  amplitude,  sensor  angle  with  respect  to  waves,  and  average  flow  past  the 
sensor  parallel  to  the  waves.  This  section  discusses  data  taken  with  the  sensor  on  a 
rotating  arm  in  still  water  to  measure  bias  and  estimate  how  bias  varies  with  relative  rms 
velocity  and  down-wave  drift.  Data  taken  in  a  wave  tank  and  during  an  open-ocean 
deployment  are  used  to  show  that  sensor  heading  does  not  significantly  affect  sensor  bias. 
I  use  data  from  a  buoy  deployment  south  of  Martha's  Vineyard,  MA.  with  no  wind  stress 
but  significant  swell,  to  estimate  bias  magnitude  in  waves. 

ROTATION  ARM  BIAS 

This  section  measures  the  vertical  and  down-wave  vorticity  biases  to  be  less  than  1 
x  10'2/s,  measures  bias  sensitivity  to  down- wave  drift,  and  shows  that  cross-wave  bias  is 
proportional  to  rms  relative  wave  velocity.  The  45-centimeter  path  sensor  was  mounted 
on  a  rotating  arm  at  the  David  Taylor  Model  Basin.  The  rotating  arm  traversed  the  sensor 
such  that  all  points  of  the  sensor  traveled  in  the  same,  but  displaced,  circles.  Any  line  in 
the  sensor  stayed  in  the  same  direction  as  the  arm  traversed. 

To  measure  vorticity  sensor  bias  in  wave-like  flow,  the  sensor  was  traversed  in 
circles  relative  to  the  carriage  at  different  rates  while  the  tow  carriage  traveled  at  different 
speeds.  A  typical  plot  of  measured  vorticity  means  for  the  same  arm  rotation  rate,  at 
different  carriage  speeds  is  shown  in  Fig.  3-15.  In  this  plot,  the  arm  was  rotated  at  a 
period  of  seven  seconds  while  the  carriage  speed  was  stepped  through  different  speeds. 
Each  point  is  a  measured  mean  of  at  least  two  minutes  of  data.  The  x-direction  vorticity  is 
across  the  tank,  parallel  to  the  axis  of  rotation  of  the  rotating  arm,  and  is  corrected  for  the 
center  stalk  bias  described  in  the  Modeling  Sensor  Errors  section.  The  y-direction 
vorticity  is  down  the  tank  parallel  to  tow  carriage  velocity  and  the  z-direction  is  vertical. 
Included  on  the  plot  are  lines  of  standard  error.  The  standard  error  of  a  mean  is  the 
standard  deviation  of  a  signal  divided  by  the  square  root  of  the  number  of  independent 
samples  of  the  signal  that  were  used  to  calculate  the  mean.  The  standard  error  is  a 
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Fig.  3-15.  Vorticity  means  of  the  45-cm  path  vorticity  meter  measured  on  the 
rotating  arm  with  a  seven  second  period.  The  standard  errors  of  each  axis  are  also 
plotted. 


measurement  of  the  variability  of  a  sample  mean.  The  standard  errors  in  these 
measurements  were  not  small  because  the  vorticity  variance  was  large.  The  Strouhal 
wakes  behind  bluff  bodies  contain  a  great  deal  of  time  varying  vorticity;  the  standard 
deviations  of  vertical  vorticity  for  the  data  taken  shown  in  Fig  3-15  were  about  0.5  per 
second.  Plots  of  vorticity  means  for  different  arm  rotation  rates  are  shown  in  Appendix 
C. 


To  measure  whether  the  down- wave  vorticity  and  vertical  vorticity  were  biased, 
the  measured  means  are  plotted  divided  by  their  standard  errors  and  compared  to  a  normal 
distribution.  The  down-wave,  y  axis,  vorticity  means  divided  by  their  standard  errors  for 
all  the  rotation  rates  and  carriage  speeds  are  plotted  in  Fig.  3-16,  and  the  vertical,  z  axis, 
results  in  Fig.  3-17.  In  both  these  figures,  the  number  of  independent  samples  was 
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Fig.  3-16.  Down- wave  vorticity  means  divided  by  their  standard  errors  compared 
to  a  normal  distribution. 


Fig.  3-17.  Vertical  vorticity  means  divided  by  their  standard  errors  compared  to  a 
normal  distribution. 
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assumed  to  be  half  the  number  of  samples  of  the  means.  An  autocorrelation  of  a  vorticity 
measurement  made  during  a  rotating  arm  test  is  shown  in  Fig.  3-18.  The  width  of  the 
autocorrelation  function  at  half  its  peak  value  is  two  samples.  Figures  3-16  and  17  show 
that  there  are  no  statistically  significant  biases  of  down-wave  or  vertical  vorticity.  If  there 
are  biases  in  either  of  these  measurements,  the  biases  are  less  than  1  *  10'2  /s. 


Fig.  3-18.  Autocorrelation  of  vorticity  measurement  during  a  rotating  arm  test 
showing  an  autocorrelation  time  scale  of  two  samples. 

To  measure  the  influence  of  the  down- wave  drift  velocity,  u ,  on  bias,  I  performed 
rotation-arm  tests  at  different  tow  carriage  speeds  and  rotation  rates.  The  data  from  these 
tests  are  plotted  in  Fig.  3-19.  The  abscissa  is  the  relative  drift  velocity  divided  by  the 
relative  wave  velocity,  and  the  ordinate  is  the  mean  vorticity  scaled  by  the  acoustic  path 
length  and  relative  velocity.  This  data  shows  that  relative,  down-wave  drift  velocity  u 
does  not  strongly  effect  bias. 

Vorticity  bias  was  found  to  scale  with  relative,  rms  velocity.  Bias  for  rotation  arm 
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tests  with  the  same  ratio  of  drift  speed  to  circular  wave  velocities  are  plotted  as  a  function 
of  rms  relative  velocity  in  Fig.  3-20.  That  errors  scale  with  relative  velocities  in  a  wave 
flow,  makes  sense  because  wake  velocities  scale  with  relative  flow  and  these  wakes 
advected  back  into  the  measurement  volume  cause  the  bias  error. 


Fig.  3-19.  Rotating-arm,  sensor  vorticity  means  plotted  with  respect  to  relative 
drift  velocity. 


Data  from  these  tests  have  to  be  used  carefully  to  infer  wave  bias  conclusions, 
because  a  rotating  arm  test  is  a  model  of  wave  flow,  and  the  rotating  arm  vibrated 
strongly.  A  rotating  arm  test  is  a  model  of  wave  flow  because  while  the  flow  moves 
relative  to  the  sensor  in  the  circular  paths  like  a  wave  flow,  the  pressure  gradients  on  the 
sensor  are  not  the  same  as  in  a  wave  field.  Pressure  gradients  that  are  parallel  to  flow, 
affect  boundary-layer  separation  and  therefore  affect  wakes  left  by  the  sensor.  Sensor 
wakes  advecting  back  into  the  measurement  volume  are  the  dominant  bias  errors. 
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Fig.  3-20.  Rotating-arm  vorticity  means  for  constant  drift  velocity  divided  by 
wave  velocity  plotted  with  respect  to  relative  rms  velocity. 

therefore  anything  that  affects  boundary  layer  separation  on  the  sensor  could  influence  the 
sensor  bias. 

The  rotation  arm  used  vibrated  strongly.  The  power  spectral  density  of  measured 
horizontal  acceleration  is  shown  in  Fig.  3-21;  it  indicates  aliasing  of  vibration  and  motion. 
These  spectra  can  be  compared  to  the  acceleration  spectra  from  an  ocean  deployment  Fig. 
3-22.  The  buoy's  most  effective  anti-alias  filter  is  the  second-order,  lightly-damped 
oscillator  that  is  the  buoy  response.  The  fastest  natural  frequency  of  the  buoy  is  in  heave 
with  a  period  of  about  a  second;  pitch  and  roll  frequencies  are  slower.  In  the  ocean 
acceleration  spectra,  the  spectral  magnitude  above  three  hertz  is  the  measured  noise  level 
of  the  accelerometers.  Sensor  vibration  can  also  affect  boundary-layer  separation  and 
wake  formation.  The  buoy,  when  deployed  in  the  ocean,  does  not  have  aliasing  problems, 
but  when  strapped  rigidly  to  the  vibrating  rotating  arm  did  have  aliasing. 


Ill 
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Fig.  3-21.  Horizontal  acceleration  spectra  from  rotating  arm  test  showing  strong 
vibration. 


Fig.  3-22.  Open-ocean  deployment  acceleration  spectra;  X  is  horizontal 
acceleration,  and  Z  is  vertical  acceleration. 


WA  VE  BIAS  SENSITIVITY  TO  ORIENTA  TION 

Sensor  bias  was  measured  to  not  vary  significantly  with  sensor  orientation  in  a 
wave  tank.  Wave  tanks  have  real  wave  flow  which  avoids  the  unmodeled  pressure 
gradients  of  a  rotating  arm  test.  The  actual  shear  drifts  can,  however,  be  non-zero.  In  a 
wave  tank,  the  sensor  is  held  still  while  waves  radiate  past.  The  buoy  in  the  ocean, 
however,  moves  with  the  longer  wavelength  swell  and  wave  flow  relative  to  the  buoy  is 
more  complicated.  Again,  I  inferred  the  sensitivity  of  an  error  parameter  from  a  test  that 
does  not  precisely  reproduce  flow  parameters. 

The  wave  tank  used  was  in  the  Canada  Center  for  Inland  Waters  and  when  used 
carefully,  was  stationary  and  measurements  were  repeatable.  Tank  flow  measurements 
had  to  be  averaged  for  at  least  ten  minutes  to  get  consistent  measurements.  Fig.  3-23 
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Fig.  3-23.  CCIW  wave  tank  velocity  measurements  made  with  an  acoustic  current 
meter  showing  low  frequency  variability.  The  velocity  trace  in  this  long  time  series 
of  a  rapidly  varying  signal,  is  smeared  and  aliased  by  the  printing  but  does  show 
the  time  varying  envelope  of  tank  velocity.  The  sampling  at  10  Hz.  of  a  wave  of 
a  several  second  period  is  not  aliased. 
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shows  a  velocity  measurement  that  had  significant  low  frequency  variability.  Ten-minute 
averages  provided  consistent  mean  flow  measurements  to  about  1  mm/s,  even  when  the 
elapsed  time  between  measurements  was  over  three  hours. 

The  wave  tank  testing  took  advantage  of  the  tank's  stationarity.  The  mean  flow 
pump  and  wave  maker  were  started  and  run  at  least  two  hours  for  tank  transients  to  die 
down  before  any  measurements  were  made.  With  the  vorticity  sensor  out  of  the  water, 
ten-minute,  mean  current  measurements  were  made  at  the  top  and  bottom  depths  of  the 
vorticity  measurement  volume.  BASS  and  Denshi  Kogyo  (a  Japanese  instrument  maker) 
acoustic  current  meters  were  used.  The  current  meters  were  removed  from  the  water,  the 
vorticity  sensor  was  put  in  the  water,  and  at  least  a  ten-minute  measurement  of  vorticity 
was  made  without  wake  disturbance  from  the  current  meters.  The  vorticity  meter  was 
removed  from  the  water,  the  current  meters  put  back  in,  and  ten-minute  current 
measurements  were  made  again  at  the  top  and  bottom  of  the  measurement  volume  to  see 
if  the  wave  tank  drift  currents  had  changed  over  the  elapsed  time.  The  wave  tank  was 
found  by  all  the  instruments  to  be  stationary  but  had  large  shears,  of  order  0. 1 1  /s.  This 
shear  is  larger  than  typical  for  unstratified  geophysical  flows.  Current  meters  have 
unknown  errors  in  wave  flows  and  the  BASS  and  Denshi  Kogyo  current  meter 
measurements  of  shear  often  differed  by  more  than  0.04  /s.  The  wave  tank  was  too  wide 
to  use  a  laser  Doppler  velocity  meter. 

The  consistency  of  current  meter  measurements  did  however  show  tank 
stationarity  that  was  used  to  measure  how  vorticity  meter  bias  changes  with  orientation. 
Under  one  set  of  wave  and  current  conditions,  the  vorticity  sensor  was  measured  at 
different  heading  angles  relative  to  the  waves.  The  mean  vorticity  measurements  were 
rotated  into  tank  coordinates  and  are  listed  in  Table  3-2.  Because  of  sensor  symmetry 
about  the  vertical  axis,  the  three  physical  angles  are  repeated  every  120  degrees  and 
mirrored  every  60  degrees.  Sampling  at  these  three  angles  is  equivalent  to  sampling  every 
15  degrees  about  the  vertical  axis.  No  changes  in  measured  horizontal  vorticity  greater 
than  0.014  /s  were  caused  by  sensor  angle.  The  measured  vertical  vorticity,  however. 
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changed  by  0.02  /s.  A  second  measurement  of  wave  bias  sensitivity  to  orientation,  made 
in  the  open  ocean,  showed  less  change  with  orientation,  and  will  be  described  in  the  next 
section. 


Sensor  Orientation 

Mean  Vorticity 

Parallel  to  Tank  Cross  Tank 

Vertical 

0° 

-0.0093 

0.1159 

0.0251 

-15° 

-0.0102 

0.1055 

0.0101 

-30° 

-0.006 

0.1064 

0.0048 

Table  3-2.  Measured  vorticity  means  in  tank  coordinates  at  three  sensor  rotation 

angles. 


OPEN-OCEAN  BIAS  MEASUREMENT 

The  buoy  was  deployed  in  the  open  ocean  with  negligible  wind  stress  and 
significant  swell  to  measure  instrument  bias  in  waves.  The  deployment  was  south  of 
Martha’s  Vineyard  MA.  after  a  storm.  An  actual  deployment  "simulates"  all  error  terms; 
the  only  shortcoming  is  that  the  true  shear  is  not  accurately  known.  Mean  vorticity 
measured  by  the  shallow  sensor  along  with  the  standard  errors  of  each  segment  are  shown 
in  Fig.  3-24.  Recall  that  vorticity  is  the  curl  of  the  velocity  field  and  is  therefore 
perpendicular  to  the  corresponding  velocities.  Mean  vorticity  and  corresponding  standard 
errors  measured  by  the  sensor  2.45  meter  deep  are  plotted  in  Fig.  3-25.  These 
measurements  will  be  used  with  some  assumptions  to  estimate  what  instrument  bias  could 
be. 


The  dominant  wave  direction  will  be  used  along  with  the  results  from  the  rotating 
arm  tests  to  estimate  what  part  of  the  measured  vorticity  is  real  and  what  part  could  be 
bias.  The  wave  directional  spectrum  of  velocity  was  computed  using  the  maximum 
entropy  method  and  is  shown  in  Fig.  3-26  (Longuet-Higgins,  Cartwright,  and  Smith, 
1963,  and  Lygre  and  Krogstad,  1986).  Most  of  the  waves  were  radiating  at  a  heading  of 
45  degrees.  The  measured  mean  vorticity  components  were  rotated  into  down-wave  and 
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Fig.  3-24.  Mean  vorticity  and  standard  errors  measured  by  the  sensor  0.83  meter 
deep  during  deployment  south  of  Martha’s  Vineyard. 


Fig.  3-25.  Mean  vorticity  and  corresponding  standard  errors  measured  2.45  meter 
deep  during  deployment  south  of  Martha’s  Vineyard. 
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Fig.  3-26.  Directional  wave  spectrum  of  velocity  measured  on  the  south  of 
Martha’s  Vineyard  deployment,  calculated  with  the  maximum  entropy  method. 

cross- wave  components  with  the  coordinates  defined  by  Fig.  3-27,  and  are  listed  along 


North 


Fig.  3-27.  Coordinate  rotation  from  east,  north  and  vertical  to  down-wave, 
cross-wave,  and  vertical. 

with  rms  relative  velocity  in  Table  3-3.  The  rotating  arm  tests  showed  no  statistically 
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Sensor  depth 

down-wave 

cross-wave 

vertical 

rms  relative  velocity 

0.83  m 

0.016 /s 

-0.01  /s 

-0.0042  /s 

28.1  cm/s 

2.45  m 

-0.015 /s 

0.0096  /s 

-0.0047  /s 

27.9  cm/s 

Table  3-3.  Mean  measured  vorticity  rotated  into  down-wave  and  cross-wave 
coordinates  from  deployment  south  of  Martha’s  Vineyard. 

significant  bias  in  the  down-wave  or  vertical  directions.  In  this  deployment,  the  down- 
wave  and  vertical  vorticities  are  assumed  to  be  real  and  the  cross-wave  vorticities  are 
assumed  to  be  error.  The  rotating  arm  tests  also  showed  that  bias  errors  were 
proportional  to  rms  relative  velocity.  The  error  bars  presented  in  the  applications  section 
are  these  cross-wave  vorticity  means  scaled  by  rms  relative  velocity. 

The  actual  shear  for  this  deployment  could  have  been  as  large  as  0.022  /s.  Current 
shears  in  a  stratified  flow  have  very  little  vertical  shear  stress  or  mixing  if  the  gradient 
Richardson  number  is  above  0.25.  The  measured  Brunt-Vaisala  frequency  was  0.01 1 
rad/s  giving  a  9.4-minute  buoyancy  period.  This  stratification  would  have  a  Richardson 
number  of  0.25  if  the  shear  was  0.022  /s.  If,  on  one  hand,  the  actual  shear  was  this  large 
and  in  the  direction  measured  by  the  buoy,  the  bias  could  have  been  zero.  If,  on  the  other 
hand,  the  actual  shear  was  this  large  in  the  opposite  direction  of  the  buoy  measured 
vorticities,  the  bias  could  have  been  larger. 

If  vorticity  is  assumed  to  be  constant  for  a  20  minute  period,  the  fact  that  the  buoy 
rotated  can  be  used  to  check  the  presumed  zero  bias  in  down-wave  and  vertical  vorticity. 
Measured  vertical  vorticity  is  looked  at  first.  A  20  minute  time  series  of  measurements 
was  broken  up  into  five  degree  heading  bins  and  all  measurements  when  the  buoy  heading 
was  within  a  heading  bin  were  grouped  together.  The  mean  measured  vertical  vorticity, 
standard  deviation,  and  standard  error  were  computed  for  each  buoy  heading  bin  and 
plotted  in  Fig.  3-28.  The  autocorrelation  of  vertical  vorticity  had  a  time  scale  of  5 
samples;  the  number  of  samples  used  to  compute  standard  error  was  divided  by  5  to 
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Fig.  3-28.  Vertical  vorticity  broken  up  into  buoy  heading  bins  for  a  20  minute 
record.  Heading  bin  means,  standard  deviations,  and  standard  errors  are  shown 
and  indicate  no  significant  bias  that  varies  with  heading. 

estimate  the  number  of  independent  samples.  The  sample  means  of  each  heading  bin  do 
not  show  any  statistically  significant  bias  that  varies  with  buoy  heading.  The  bias  in 
vertical  vorticity  does  not  change  with  heading  and  the  buoy  and  vorticity  sensors  have 
radial  symmetry  so  the  bias  of  vertical  vorticity  must  be  zero. 

The  presumed  zero  bias  in  down-wave  vorticity  was  checked  by  breaking 
measurements  of  horizontal  vorticity  up  into  heading  bins  and  assuming  the  actual 
vorticity  was  constant.  If  the  actual  vorticity  is  constant  and  the  buoy  rotates,  the  buoy 
referenced  vorticity  will  be  a  sine  wave  when  plotted  with  respect  to  heading.  In  Fig.  3- 
29.  the  heading  bin  means  and  standard  deviations  are  plotted.  Also  plotted,  as  dotted 
lines,  are  the  sine  wave  from  the  20  minute-mean,  earth-referenced  vorticity  rotated  into 
buoy  coordinates,  with  the  standard  error  added  and  subtracted  from  this  sine  wave.  The 
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Fig.  3-29.  Horizontal  vorticity  in  buoy  coordinates,  averaged  in  buoy  heading  bins, 
along  with  the  standard  deviations  of  vorticity  and  the  20  minute  mean  vorticity 
rotated  into  buoy  coordinates  for  each  heading  plus  and  minus  the  standard  error. 

horizontal  vorticity  had  an  autocorrelation  time  scale  of  one  sample.  The  bin  averaged 
vorticity  means  in  buoy  coordinates  do  not  show  any  statistically  significant  biases 
(differences  from  the  sine  wave  resulting  from  the  mean  earth-referenced  vorticity  rotated 
into  buoy  coordinates)  that  vary  with  heading.  No  bias  of  horizontal-vorticity  variation 
with  heading,  combined  with  the  radial  symmetry  of  the  buoy,  indicate  that  there  is  no 
significant  down-wave  vorticity  bias. 

No  inference  of  cross-wave  vorticity  bias  can  be  made  from  breaking  buoy 
referenced  vorticity  into  heading  bins.  Recall  that  cross-wave  vorticity  corresponds  to 
velocity  in  the  plane  of  orbital  wave  motion.  Vertical  asymmetry  of  the  buoy  and  vertical 
asymmetry  of  wave  velocity  (wave  attenuation  with  depth)  along  with  free-surface  effects 
prevent  an  inference  of  no  cross- wave  bias  that  is  independent  of  buoy  heading.  From  Fig. 
3-29.  we  can  infer  that  cross- wave  vorticity  bias  does  not  change  significantly  with  buoy 
heading  but  we  cannot  infer  that  there  is  no  cross-wave  vorticity  bias. 
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WAVE  SPECTRAL  REJECTION  OF  VORTICITY  SENSORS 
This  section  measures  the  wave  spectral  rejection  in  a  laboratory  wave  tank  and  in 
a  bottom  deployment  in  the  coastal  ocean.  Buoy  flexing  causes  a  measured  spectral  peak 
which  is  measured  and  explained. 


LABORA  TORY-WA  VE  SPECTRAL  REJECTION 

The  prime  motivation  to  measure  vorticity  was  to  filter  out  irrotational,  surface- 
gravity  waves.  This  section  measures  how  well  this  goal  has  been  met.  The  wave  spectral 
rejection  of  the  1 5-centimeter  path  vorticity  sensor  prototype  was  measured  in  a  small 
wave  tank.  Time  series  of  the  sensor  measurements  are  shown  in  Fig.  3-30.  The  top 
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Fig.  3-30.  Wave  tank  time  series  of  velocity,  shear,  and  circulation  over  path 
length. 
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graph  is  the  measured  velocity  of  one  path  of  a  circulation  path  square  and  had  a  standard 
deviation  of  2.06  cm/s.  The  middle  graph  is  the  measured  velocity  difference  between  two 
parallel  paths  of  the  circulation  path  square  and  had  a  standard  deviation  of  1 .2  cm/s.  the 
bottom  graph  is  the  circulation  divided  by  the  path  length  which  had  a  standard  deviation 
of  0.27  cm/s.  Low  frequency  internal  waves  and  turbulence  were  evident  from  the 
circulation  trace  that  required  low-pass  filtering  of  the  velocity  data  to  be  evident.  The 
same  data  is  plotted  in  the  frequency  domain  in  Fig.  3-3 1 .  The  sensor  has  a  better  than  20 
dB  rejection  of  wave  velocity.  Most  of  the  noise  floor  in  this  spectrum  is  residual 
turbulence  in  the  tank  from  waves,  the  wavemaker,  and  sensor  wakes.  The  electronic 
noise  floor  in  this  measurement  was  1.8  *  104  cm2/s2Hz.  The  45-centimeter  path  vorticity 


Fig.  3-31.  Vorticity  meter  spectral  rejection  of  waves,  avel,l  and  avel,2  are  two 
acoustic  velocity  paths  of  the  circulation  square. 


sensors  rejected  surface-wave  velocities  by  30  dB. 


122 


OCEAN-WA  VE  SPECTRAL  REJECTION 

A  bottom  deployment  of  the  1.5-meter  path  vorticity  sensor  shows  its  spectral 
wave  rejection.  In  May  of  1994,  the  instrument  was  deployed  near  Woods  Hole,  MA.,  in 
an  area  of  Vineyard  Sound  having  large  tidal  currents.  The  data  shown  in  Fig.  3-32  were 
measured  during  a  storm  with  surface  swell  that  reached  the  bottom. 


Fig.  3-32.  1.5  m  benthic  vorticity  meter  power  spectral  densities  of  acoustic  path 
water  velocities  and  circulation  divided  by  the  sum  of  the  path  lengths 

The  measuring  volume  was  centered  at  1.75  meters  above  the  bottom.  During  the  time 
period  from  which  these  spectra  were  computed,  the  mean  current  was  22.2  cm/s  and  the 
rms  horizontal  wave  velocity  at  the  sensor  was  5.5  cm/s.  The  power  spectral  densities  of 
the  velocities  of  the  individual  acoustic  paths  making  up  a  square  are  shown  as  dotted  lines 
and  the  power  spectral  density  of  circulation  is  shown  as  solid  lines.  The  power  spectral 
density  of  circulation  over  path  length  has  been  scaled  by  0.25  to  give  the  same  energy  at 
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high  frequencies  as  the  velocity  spectra.  These  data  were  reported  by  A.  J.  Williams, 
Terray,  Thwaites,  and  Trowbridge.  (1994).  This  instrument's  effective  surface  wave 
rejection  shows  the  utility  of  measuring  vorticity. 

The  spectra  shown  in  Fig.  3-32  are  from  time  series  and  cannot  be  converted  to  a 
wave  number  spectra  by  the  usual  frozen  turbulence  model.  The  usual  test  for  when  the 
frozen  turbulence  hypothesis  can  be  used,  of  the  velocity  standard  deviation  being  less 
than  half  the  mean,  is  not  met  (Stull,  1989).  These  sensors  measure  an  area  averaged 
vorticity,  not  point  vorticity.  When  the  frozen  turbulence  assumption  can  be  used  to 
transform  frequency  spectra  into  wavenumber  spectra,  for  eddies  with  a  spatial  scale  more 
than  twice  the  vorticity  meter  size,  the  spectral  transfer  function  should  be  one.  For 
eddies  with  spatial  scales  smaller  than  twice  the  sampling  volume,  line  averaging  and  area 
averaging  will  cause  attenuation.  For  eddies  with  spatial  scales  less  than  a  third  the  path 
length,  the  velocity  paths  will  be  uncorrelated  and  should  have  a  spectral  response  like  a 
travel-time  acoustic  velocity  meter.  There  is  a  tradeoff  between  spatial  resolution  of  a 
measurement  and  the  accuracy  of  that  measurement. 

BUOYFLEXING 

The  double-sensor,  vorticity  measuring  buoy  is  a  tall  flexible  structure  that  has  a 
bending  resonance  that  can  be  seen  in  deployment  data,  Fig.  3-33.  The  buoy  was  built  to 
minimize  drag  in  the  measurement  volumes  by  keeping  structural  elements  as  slim  as 
practical.  The  resonance  peak,  at  just  below  three  Hertz,  is  in  the  measured  horizontal 
velocity  but  not  the  vertical  velocity.  This  peak  is  enlarged  in  Fig.  3-34.  Compensating 
for  buoy  rotation  with  the  rate-gyro  in  computing  horizontal  velocity,  cannot  account  for 
buoy  flexing,  leaving  this  error  in  measured  earth  referenced  velocity.  This  flexing 
resonance  could  only  be  reduced  by  making  the  structure  larger,  which  would  cause  more 
drag  and  flow  disturbance.  However,  this  peak  is  small  compared  to  wave  energy  and  can 
be  ignored.  Deployments  of  the  single  sensor  buoy  did  not  show  this  peak.  The  single 
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Fig.  3-33.  Measured  power  spectral  density  of  east  and  vertical  velocity.  The  buoy 
flex  resonance  is  on  the  lower  right  side  of  the  graph. 


Fig.  3-34.  Enlarged  view  of  measured  spectral  peak. 


sensor  buoy  has  a  much  higher  flexing  natural  frequency,  and  at  this  higher  frequency  the 
wave  excitation  forces  are  smaller. 

This  chapter  modeled  and  measured  vorticity  meter  accuracy.  A  potential  flow 
model  was  developed  that  predicted  a  calibratable  reduction  in  gain  but  no  bias  or  noise. 
Wake  models  were  developed  to  predict  bias  in  constant  flow.  Electronic  noise  and  bias 
were  shown  to  be  much  smaller  than  flow  disturbance  related  error.  Bias  in  constant  flow 
was  found  be  less  than  1  x  10‘2  /s  with  the  45-centimeter  path  sensor  and  variance  in 
constant  flow  was  measured.  Bias  in  wave  flow  was  measured  to  be  less  than  1  x  10'2  /s 
in  the  down- wave  and  vertical  direction  and  about  1  *  10‘2  /s  in  the  cross- wave  direction 
in  open-ocean  conditions.  Rejection  of  surface  waves  was  measured  and  showed  the 
utility  of  measuring  vorticity  when  a  measurement  of  either  turbulence  or  shear  is  desired 
in  the  presence  of  waves.  The  vorticity  measuring  buoy  was  shown  to  be  able  to  make 
shear  measurements  in  the  upper  five  meters  of  the  ocean  over  finer  resolutions  than  has 
been  practical  before. 

The  major  vorticity  sensor  error  sources  are  reviewed  along  with  how  the  errors 
scale  with  sensor  size  and  whether  the  error  is  noise  or  bias,  in  Table  3-4.  In  this  table,  d 
refers  to  the  outer  diameter  of  the  pod  spheres,  D  refers  to  the  center  stalk  diameter,  and  L 
refers  to  an  acoustic  path  length.  The  transducer  pod  diameters  for  the  1 5  cm  path  and  45 
cm  path  vorticity  meters  are  7/8  inch  (2.22  cm)  and  1  1/8  inch  (2.86  cm)  and  the  center 
tube  diameters  are  Vz  inch  (1.27  cm)  and  1  Vz  inch  (3.81  cm)  respectively. 
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FLOW  ERROR 

VELOCITY 

VORTICITY 

NOTES 

SOURCES 

SCALE 

SCALE 

Electronic  noise 

1/L 

1/L2 

noise,  rms  0. 1  mm/s 

Zero  drift 

1/L 

1/L2 

bias,  1/3  cm/s  if  casual 
1  mm/s  if  careful 

Sound  speed 

(c,2  -  c 

U  d 

r2)/ c,2  xgain 

change  in  gain 

1°  C  -  0.6%  error  in 
scale,  can  be 
compensated 

Potential  flow 

oo 

21 

UD  tan  0 

cos  0 

UD  tan  6 

reduction  in  gain 

7.4  %  for  prototype 

0  =  0 

Wake  2D 

bias 

center  strut  wake 

L 

TT  (  d\\ 

L2 

ujd\  4 

bias 

Wake  3D 

U  —  3 

pod  wake 

l  Lj 

z  UJ 

Lift  related 

U„D/L 

U„D/L2 

bias 

Advection  of  wake 

function  of  wave  field  and 
sensor  motion 

bias  and  noise 

0.01  /s  in  28 
cm/s  relative 
velocity 

Vortex  shedding 

1.5XU„D/L 

1.5xU„D/L2 

noise 

Table  3-4  Major  sensor  error  sources 
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CHAPTER  4.  APPLICATIONS 


Chapter  four  describes  two  deployments  of  the  vorticity  measuring  buoy  to 
demonstrate  the  potential  of  the  instrument  to  measure  shear  and  vorticity  over  finer 
resolutions  than  before  practical.  The  first  application  measured  shear  in  the  upper-ocean 
boundary-layer  within  5  meters  of  the  surface.  The  second  application  measured  gradient 
Richardson  number  with  fine  resolution  in  a  thermocline.  An  application  of  measuring 
turbulence  with  a  vorticity  meter  in  the  bottom-boundary  layer  was  described  in  the 
Ocean-Wave  Spectral  Rejection  section. 


4.A.  UPPER-BOUNDARY-LAYER  SHEAR 

The  first  application  of  the  double  45-centimeter  path,  shear-measuring  buoy  to  be 
discussed  is  measurement  of  shear  in  the  upper-ocean,  boundary  layer.  Shear,  turbulence, 
and  mixing  in  unstratified,  turbulent,  boundary  layers  over  rigid  walls  have  been  well 
studied  and  are  well  understood.  By  comparing  the  shear  in  an  ocean  boundary  layer  to 
the  stress  over  this  layer,  the  effect  on  mixing  of  the  free  surface,  stratification,  and 
possible  organized  motions  like  Langmuir  circulation  can  be  measured.  The  shear 
expected  in  the  log  layer  next  to  a  rigid  wall  is  given  by  equation  (4-1),  where  dUtdz  is  the 

dU  u . 


shear,  a,  is  the  friction  velocity  (the  square  root  of  the  shear  stress  divided  by  density),  k  is 
von  Karman's  constant  and  is  usually  assumed  to  be  0.4,  and  z  is  the  distance  from  the 
wall. 


The  shear-measuring  buoy  was  deployed  in  Buzzards  Bay,  Massachusetts,  in  May 
1994.  The  measured  shears  are  compared  to  the  shear  one  would  expect  in  an 
unstratified,  rigid  wall,  turbulent  flow  of  the  same  shear  stress,  Fig.  4-1  (Monin  and 
Yaglom,  1982).  The  wind  stress  was  concurrently  measured  by  a  sonic  anemometer  and 


129 


calculated  using  the  inertial  dissipation  method  (Fairall  and  Larsen,  1986).  The  wind 
during  the  deployment  increased  in  strength  from  about  5.2  m/s  to  about  15.2  m/s,  causing 
the  water-side  friction  velocity  to  increase  from  0.62  cm/s  to  1.82  cm/s.  The  error  bars 


Fig.  4-1.  Buzzards  Bay  deployment  shears  measured  and  estimated  from  wind 
stress.  The  upper  graph  is  the  sensor  at  0.83  m  depth  and  the  lower  graph  is  the 
sensor  at  2.45  m  depth. 
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are  derived,  as  described  in  the  Bias  in  Wave  Flow  section,  from  an  open-ocean 
deployment  with  significant  swell  and  negligible  wind  stress,  and  are  scaled  with  rms 
velocity  relative  to  the  sensors.  The  sensor  at  the  shallower  depth  of  0.83  meters, 
measured  shear  that  coincided  with  the  log-layer  shear  prediction  for  most  of  the 
deployment.  Toward  the  end  of  the  deployment,  as  the  windstress,  waves,  and  wave 
breaking  increased,  the  shear  started  to  drop  and  diverge  from  the  wall  model.  The 
vertical  heave  velocity  spectra  at  the  start  and  end  of  the  deployment  are  shown  in  Fig.  4- 
2.  The  vertical  heave  velocity  variance  increased  from  100  cm2/s2  at  the  start  to  252 
cm2/s2  at  the  end. 


Fig.  4-2.  Vertical  velocity  spectra  at  the  start  and  end  of  the  Buzzards  Bay  Ma 
deployment.  The  heave  velocity  variance  increased  from  100  cm2/s2  to  252  cm2/s2. 
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The  sensor  at  the  deeper  average  depth  of  2.45  meters,  measured  about  three  times 
more  shear  than  would  be  predicted  by  a  log  layer.  Average  temperature  stratification 
measured  at  2.45  meters  depth  gave  a  Brunt-Vaisala  frequency  of  0.012  Is  (period  of 
about  nine  minutes),  and  the  measured  gradient  Richardson  number  at  this  depth  dropped 
from  0.20  to  0. 10.  Stratification  was  inferred  from  temperature  assuming  that  salinity  was 
well  mixed  in  the  upper  five  meters.  Surveys  of  stratification  done  with  a  conductivity, 
temperature,  depth  instrument  (CTD)  in  this  part  of  Buzzards  Bay  during  this  time  of  year 
have  shown  salinity  to  be  well  mixed  over  the  upper  five  meters  and  the  error  in 
stratification  resulting  from  not  measuring  salinity  should  be  less  than  twenty  percent.  The 
Richardson  number  time  series  are  plotted  in  Fig.  4-3.  These  are  reasonable  Richardson 


Fig.  4-3.  Gradient  Richardson  numbers  measured  during  the  Buzzards  Bay 
deployment. 

numbers  for  a  stressed  boundary  layer  with  turbulent  mixing.  A  more  detailed  vector 
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breakdown  of  these  results  into  directions  downwind  and  crosswind  is  in  Appendix  C. 

The  increased  shear  of  the  sensor  at  2.45  meters  as  compared  to  the  predicted  shear  of  an 
unstratified  boundary  layer  is  a  measurement  of  the  stratification's  inhibition  of  vertical 
turbulent  mixing  and  consequent  enhancement  of  shear. 

The  measured  data  from  this  deployment  will  now  be  replotted  in  vector 
components  along  with  shears  to  show  the  sensor-wave-correlation  bias.  As  the  Sensor- 
Wave-Correlation  Bias  section  described,  velocity  or  shear  derived  from  velocity 
measured  from  a  platform  that  moves  in  a  correlated  way  with  the  waves,  can  have  a  bias 
that  can  be  as  large  as  the  Eulerian  shear.  The  vorticities  and  shears  are  plotted  with 
respect  to  the  dominant  wave  direction.  The  measured  wave  velocity  directional  spectrum 
calculated  by  the  maximum  entropy  method  is  plotted  in  Fig.  4-4.  The  wind  was  from 


Fig.  4-4.  Wave  velocity  directional  spectrum  calculated  by  the  maximum  entropy 
method  for  the  Buzzards  Bay  deployment. 
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205°  but  the  bay  is  oriented  east-west  on  a  magnetic  compass  causing  the  longer  period 
waves  to  be  from  the  west.  The  shorter  period  waves  were  from  the  southwest.  The 
peakiness  of  the  directional  spectrum  may  be  an  artifact  of  the  maximum  entropy  method; 
both  the  heave  spectrum  and  the  Longuet-Higgins  spectrum  show  only  two  peaks. 

The  vorticity  and  shear  measured  by  the  shallow  sensor  at  0.83  meters  are  plotted 
in  Fig.  4-5  and  by  the  deeper  sensor  at  2.45  meters  in  Fig.  4-6.  These  are  plotted  in  east 
and  north  coordinates  that  correspond  to  the  longer  period  waves.  Because  all  the  waves 
were  not  from  the  same  direction,  it  was  not  possible  to  rotate  the  measurements  to 
correspond  with  all  the  waves.  The  shears  measured  were  from  the  velocity  paths  on  the 
top  and  bottom  of  the  vorticity  measurement  volumes  and  in  the  shallow  sensor  measure 


Fig.  4-5.  Vorticity  and  shear  corresponding  to  that  direction  of  vorticity,  measured 
by  the  sensor  0.83  m  deep  during  Buzzards  Bay  deployment. 
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significantly  more  shear  in  the  down-wave  direction.  The  shallow  sensor’s  shear 
corresponding  to  the  north  vorticity,  was  1.8  x  10'2  /s  greater  than  north  vorticity.  This 
corresponds  with  the  predominant  waves  radiating  east  and  is  a  measurement  of 
correlation-wave  bias.  Shear  and  vorticity  in  the  other  direction  and  measurements  in  the 
lower  sensor  have  differences  that  are  smaller  than  the  error  bars  but,  in  all  but  the  lower¬ 
sensor,  east  vorticity,  are  of  the  right  sign.  The  lower-sensor,  east  vorticity  only  differs 
from  its  corresponding  shear  by  0.2  *  10'2  /s.  Removing  correlation- wave  bias  is 
important  when  measuring  shear  from  buoys  in  the  upper  five  meters  of  the  ocean,  and 
this  instrument  does  so  intrinsically. 


Fig.  4-6.  Vorticity  and  shear  corresponding  to  that  direction  of  vorticity, 
measured  by  the  sensor  2.45  m  deep  during  Buzzards  Bay  deployment. 
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4.B.  THERMOCLINE  RICHARDSON  NUMBER 

The  second  application  of  the  shear  measuring  buoy  to  discuss  is  measurement  of 
shear  and  local  shear  instability  in  an  internal  boundary  layer.  In  August  of  1993,  the  shear 
measuring  buoy  was  deployed  with  a  single  sensor  in  the  thermocline  in  Massachusetts 
Bay.  Fig.  4-7  shows  a  typical  time  history  of  temperature,  magnitude  of  vertical  shear, 
and  gradient  Richardson  number. 

The  histogram  of  the  gradient  Richardson  number  is  shown  in  Fig.  4-8.  In  this 
figure,  the  histogram  of  the  Richardson  number  measured  over  a  1.0-meter  bin  size  is  also 
shown.  The  latter  measurements  were  made  the  same  day  using  an  acoustic  doppler 
current  profiler  (ADCP)  to  measure  shear  and  a  conductivity  temperature  depth 
instrument  (CTD)  to  measure  stratification.  The  averaging  time  for  the  1.2MHz  ADCP 
was  30  seconds  for  Fig.  4-8,  and  the  shear  for  both  Figs  4-7  and  4-8  was  lowpass  filtered 
with  a  25-second  period  filter.  There  is  a  greater  proportion  of  low  Richardson  number 
measurements  at  the  smaller  vorticity  meter  38.6  centimeter  resolution  than  the  ADCP 
1 .0  meter  resolution.  This  difference  is  in  qualitative  agreement  with  the  theory  by 
Desaubies  and  Smith  (1982)  that  predicts  more  small  mixing  events  than  large  mixing 
events  in  a  random,  internal- wave  field.  Twenty-five  percent  of  the  Richardson  numbers 
measured  over  the  38.6-centimeter  volume  were  unstable  while  12.7%  of  the  Richardson 
numbers  measured  over  the  one  meter  volume  were  unstable.  The  Richardson  number 
measurements  were  made  in  the  middle  of  a  dye  tracer  release  experiment  that  measured  a 
vertical  diffusivity  of  0.06  cm2/s  (Geyer  and  Ledwell,  1994).  This  deployment  shows  that 
the  instrument  can  measure  shear  and  local  shear  instability  in  ocean  internal-boundary 
layers. 
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Richardson  Number 


18 


Massachusetts  Bay  Temperature,  Shear,  and  Gradient  Richardson  Number 


Fig.  4-7.  Gradient  Richardson  number  measured  in  the  Massachusetts  Bay 
thermocline.  In  the  top  graph  of  temperature,  the  dot  dash  line  is  the  temperature 
at  top  of  the  measurement  volume  and  the  solid  line  is  at  the  bottom.  In  the 
Richardson  number  plot,  a  dotted  reference  line  is  drawn  at  0.25. 
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Normalized  Histogram 


Fig.  4-8.  Relative  frequency  of  Richardson  numbers  measured  by  the  vorticity 
meter  over  38.6  cm  and  by  an  ADCP  over  1.0  m. 
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CHAPTER  5.  SUMMARY 


5.A.  CONCLUDING  REMARKS 

The  major  goal  of  my  research  has  been  to  develop  an  instrument  system  capable 
of  measuring  shear  in  the  upper  five  meters  of  the  ocean.  The  data  measured  and 
presented  in  this  thesis  are  a  test  of  the  instrument.  This  is  the  first,  three-axis 
oceanographic  vorticity  instrument  deployed  in  the  ocean.  The  design  has  already  been 
used  in  a  subsequent  project  to  measure  coastal,  bottom-boundary-layer  turbulence  in  the 
presence  of  surface  swell. 

To  remove  platform  movement  from  flow  measurements,  I  used  a  low-power, 
strapdown,  inertial-measurement  unit  for  the  instrument  system.  Constraints  of  buoy 
motion  allowed  the  use  of  an  Euler-angle  algorithm  with  complementary-filtering  that 
gave  good  performance  results  from  slow  samples  of  noisy,  drifty,  and  low-power  inertial 
sensors. 

The  instrument  is  suitable  for  measuring  shear  in  the  upper  five  meters  of  the 
ocean,  measuring  fine  scale  shear  and  vorticity  in  internal-boundary  layers,  and  coastal, 
bottom-boundary-layer  turbulence  in  the  presence  of  surface  swell.  Measurement  errors 
are  dominated  by  asymmetries  between  sensor  wakes  and  circulation  paths.  This  analysis 
has  shown  the  tradeoff  between  spatial  resolution  and  accuracy  in  measuring  shear  or 
vorticity. 

In  the  first  chapter,  I  showed  the  need  for  measurements  of  shear  in  the  upper  five 
meters  of  the  ocean,  the  desirability  of  filtering  out  irrotational  water  motion  such  as 
surface  gravity  waves  when  measuring  shear,  and  the  need  to  remove  or  compensate  for 
sensor  motion  when  measuring  shear.  Sensor-wave-correlation  bias  of  velocity  and 
velocity-derived  shear  measurements  was  reviewed  and  non-intrusive  measurements  of 
vorticity  were  shown  to  not  have  this  bias.  The  only  reason  to  average  vorticity 
measurements  is  to  reduce  fluctuating  wake  error  effects. 
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In  the  second  chapter,  I  motivated  the  mechanical  design,  described  the  electrical 
system  and  inertial  measurement  unit,  reviewed  rotation  and  strapdown  inertial 
computation,  described  the  inertial  algorithm  used,  and  described  simulations  of  the 
inertial  system.  In  measurements  of  shear  or  vorticity  with  this  instrument,  error 
contribution  from  imperfect,  inertial  compensation  is  much  smaller  than  wake-related 
errors. 


The  third  chapter  developed  models  of  and  described  measurements  of  the 
instrument's  major  error  sources.  Electronic  noise  and  drifts  were  shown  to  be  much 
smaller  than  flow-related  errors.  The  simpler,  triangular  circulation  path  was  shown  to 
have  much  greater  bias  in  constant  flow  than  the  chosen  square  circulation  path.  Variance 
in  constant  flow  was  shown  to  result  from  vortex  shedding  off  bluff  sensor  parts.  Simple 
attempts  to  streamline  the  sensor  were  shown  to  be  capable  of  increasing  instrument 
biases.  Laboratory  measurements  that  simulated  waves  showed  down-wave  and  vertical 
biases  of  less  than  1  *  10'2  /s,  that  cross- wave  bias  varied  linearly  with  rms  relative 
velocity  and  did  not  vary  stongly  with  down-wave  drift  velocity.  Cross-wave  vorticity 
bias  in  wave  flow  was  measured  in  an  ocean  deployment  with  no  wind  stress  and 
significant  swell  to  have  a  bias  as  large  as  1  x  10"2  /s.  Open-ocean  measurements  showed 
that  vertical  and  horizontal  vorticity  bias  did  not  change  significantly  with  sensor  heading. 
The  sensor's  spectral  surface-gravity-wave  rejection  was  measured  in  both  wave  tanks  and 
in  the  coastal  ocean,  and  showed  the  utility  of  measuring  vorticity. 

The  fourth  chapter  described  two  ocean  deployments  of  the  instrument. 
Measurements  of  upper-boundary-layer  shear  were  compared  to  models  of  turbulent, 
boundary-layer  behavior.  Sensor-wave-correlation  bias  of  velocity-derived  shear  was 
measured  to  be  as  large  as  Eulerian  Shear.  Measurements  of  the  Richardson  number  in 
the  thermocline  over  small  scales  were  compared  to  Richardson  numbers  measured  over 
larger  scales  and  were  consistent  with  the  theory  that  small,  unstable  internal  wave  events 
are  more  common  than  large,  unstable  events.  Data  from  these  deployments  showed  the 
vorticity  meter’s  viability  in  measuring  small-scale,  ocean-boundary-layer  shear. 
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5.B.  FUTURE  WORK 

The  vorticity  meter  and  the  shear  measuring  buoy  it  assembles  into  is  complete, 
calibrated,  and  ready  to  be  redeployed.  Uses  of  the  system  in  future  deployments  include 
longer  deployments  measuring  shear  over  several  wind  shifts  and  measurements  of 
Langmuir  cells.  The  vorticity  sensor  requires  shorter  averaging  times  than  measurements 
of  shear  with  velocity  meters;  this  shorter  averaging  time  allows  slow  traverses  across 
Langmuir  cells.  The  sensor  would  have  to  be  scaled  close  to  the  size  of  the  cells, 
pathlengths  of  about  a  quarter  of  the  cell  depth  would  be  appropriate  for  adequate 
sampling  with  good  signal-to-noise  ratio. 

If  ocean  scientists  desired  the  ability  to  measure  low-frequency  wave  spectra, 
possible  future  system  improvements  would  include  installing  more  accurate  rate  gyros, 
sampling  faster,  and  or  adding  position  or  velocity  measurement  with  a  Kalman  filter. 
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APPENDIX  A. 


Mechanical  Drawings 

A-l  Shear  measuring  buoy  assembly 

A-2  Vorticity  sensor  assembly 

A-3  Pressure  vessel  assembly 

A-4  Lower  hub  detail 

A-5  Upper  hub  detail 

A-6  Custom  end  cap  detail 

A-7  Instrument  electronics  case  detail 

A-8  Transducer  mount  detail 

A-9  Central  stalk,  pod  arms,  and  braising  jig 

The  major  mechanical  drawings  for  the  vorticity  sensor  are  included  in  appendix  A. 
All  the  drawings  scales  have  been  reduced  to  fit  on  8  'A*  1 1  inch  paper.  A  vorticity  sensor 
of  another  path  length  could  be  readily  designed  by  scaling  these  drawings.  On  the 
assembly  drawings,  the  out  of  plane  pod  arms  are  not  shown.  The  lower  end  cap  is  a 
standard  BASS  end  cap.  The  braising  jig  rods  are  screwed  into  each  acoustic  path  for 
braising  with  screws  with  modified  ends  to  fit  in  the  tight  confines  of  the  pod  centers.  The 
dimensions  on  these  drawings  are  in  inches  which  was  required  to  get  reasonable 
machining  bid  prices. 
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Fig.  A-7.  Instrument  electronics  case  detail  drawing 
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APPENDIX  B. 


Appendix  B  lists  all  the  processing  files  for  the  shear  measuring  buoy,  the 
processing  is  done  in  Matlab  and  follows  the  block  diagram  in  Fig.  2-20.  The  main 
program  is  vortcasv.m  which  calls  the  various  subfunctions  and  scripts.  Remdfc.m  is  a 
script  file  that  removes  the  first  and  last  data  records  in  each  229KB  disk  file.  There  is  a 
short  time  gap  of  about  fifteen  seconds  of  unrecorded  data  when  the  data  logger  writes  to 
disk;  these  disk  writes  are  not  coordinated  with  the  start  or  end  of  a  data  record.  The 
unpacking  program  Bunt.exe  can  give  a  corrupted  first  or  last  record  of  a  disk  file  when 
given  incomplete  data  records.  The  first  and  last  minute  of  processed  data  of  each  diskfile 
should  be  ignored  to  avoid  filter  transients.  Processing  variables  are  listed  and  described 
below  in  table  B-l. 


Table  B-l  Data  Variables 


Data  variable 

description 

tim 

time  integers  hr:min:s:count 

podl,  pod2,  pod3. 

velocity  measurements  for  each  circulation  path 

pod4,  pod5,  pod6 

integers 

temp 

thermistor  output  in  A/D  integers 

imul 

accelerometer  output  in  A/D  integers 

imu2 

rate  gyro  output  in  A/D  integers 

imu3,  imu3n 

magnetometer  output  in  A/D  integers 

imulc 

calibrated  accelerometer  output  in  m/s2 

imu2c 

calibrated  rate  gyro  output  in  rad/s 

imu3c 

calibrated  magnetometer  output  in  mG 

tempc 

calibrated  thermistor  output  in  deg  C 

podlz,  pod2z,  pod3z 

velocity  path  zeros  in  A/D  integers 

pod4z,  pod5z,  pod6z 

avel,  bvel,  cvel, 

calibrated  velocity  of  each  acoustic  path  of  each 

dvel,  evel,  fVel 

circulation  path  in  cm/s 
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acirc,  bcirc,  ccirc, 
dcirc,  ecirc,  fcirc 
avort,  bvort,  cvort 
dvort,  evort,  fVort 
watvelrl,  watvelr2 
vortbl,  vortb2 
watvortl,  watvort2 
BANGLEF 
BANGLEFF 
BANGLES 

BANGLESF 

imu3d 

headslow 
heading 
eulerdot 
BANGLE 
VORT1,  VORT2 
watvell,  watvel2 

baccel 

Buoyvel 

Buoyvelf 

VEL1,  VEL2 


circulation  over  path  length  for  each  circulation  path 
in  cm/s 

area  average  vorticity  of  each  circulation  path 
in  1/s 

relative  sensor  velocity  in  buoy  x,y,z  coordinates  in  cm/s 
relative  sensor  vorticity  in  buoy  x,y,z  coordinates  in  1/s 
absolute  vorticity  in  buoy  x,y,z  coordinates  in  1/s 
integrated  rate  gyro  output  in  radians 
highpass  filtered  rate  gyro  output  in  radians 
slow  estimate  of  buoy  Euler  angles  from  accelerometers 
and  magnetometers  in  radians 
lowpass  filtered  slow  estimate  of  buoy  Euler  angles 
in  radians 

magnetometer  output  rotated  into  computational 
horizontal  frame  in  mG 
unwrapped  slow  estimate  of  heading 
unwrapped  estimate  of  heading 
Euler  angle  rates  in  radians  Is 
buoy  Euler  angles  in  radians 
earth  referenced  vorticity  of  each  sensor  in  1/s 
velocity  of  water  in  each  sensor  relative  to  accelerometer 
in  buoy  x,y,z  coordinates  in  cm/s 
buoy  acceleration  compensated  for  gravity  in  m/s2 
buoy  velocity  in  m/s  in  earth  coordinates 
high  pass  filtered  buoy  velocity  in  earth  coordinates 
in  m/s 

water  velocity  in  each  sensor  in  earth  coordinates  in  cm/s 
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Program  listing 


%  vortcasv.m  is  a  script  file  that  loads  the  raw  vorticity  meter 
%  data  with  vortload 

%  removes  change  of  tattletale  diskfile  change  wild  points 
%  switches  imu3,2  and  3  to  make  up  for  wiring  error 
%  converts  imu  data  into  physical  variables 
%  converts  thermistor  data  into  temperature  in  deg  C 
%  clears  many  of  the  no  longer  needed  variables 
%  last  revision  4/13/95 

% 

vortload 

remdfc 

imu3n=imu3  swtch(imu3); 
n2cal4 

clear  imul  imulz  imulg  imu2  imu2z  imu2g  imu3  imu3z  imu3g  dift  imu3n 

tempc=n2t4(temp); 

clear  temp  x 

svzer 

n2v 

clear  podl  pod2  pod3  podlz  pod2z  pod3z  pod4  pod4z  pod5  pod5z  pod6  pod6z 

%  Have  variables  avel  acirc  avort..  imulc  imu2c  imu3c 

%optional  line  to  clear  circulations  and  use  pack  command 

clear  acirc  bcirc  ccirc  dcirc  ecirc  fcirc 

pack 

[watvelr  1  ]=veltr45ff  avel,bvel,cvel); 

[  watvelr2]:=veltr4  5  f(dvel,  evel,  fvel) ; 
buoyangb 

%  optional  line  to  clear  buoyangl  intermediate  variables 
clear  BANGLES  BANGLEF  BANGLEFF 

headrate-[diff(BANGLESF(:,3))./.15;0]+filtfilt(Bfh,Afh,imu2c(:,3)); 

[vortb  1  ,watvort  1  ]=vortran2(avort,bvort,cvort,[imu2c(:,  1 :2)  headrate]); 
[vortb2,watvort2]=vortran2(dvort,evort,fVort,[imu2c(:,  1 :2)  headrate]); 

%  optional  line  to  clear  vel  and  sensor  vort  variables 
clear  avel  bvel  cvel  dvel  evel  avort  bvort  cvort  dvort  evort 
clear  BANGLESF 
VORT  1  =b2e(watvort  1  ,B  ANGLE); 

VORT2=b2e(watvort2,BANGLE); 

rel2abs4 

%********************************************************************* 


%  VORTLOAD  is  a  script  file  that  loads  the  .mat  files  from 
%  the  bunt  program,  and  converts  the  two's  complement  numbers 
%  from  the  BASS  fluid  velocities  into  regular  integers. 
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%  This  version  assumes  the  vorticity  measurements  are  in  pod  1 , 

%  pod2,  and  pod3  only. 

% 

load  time 

tim=x; 

load  pod  1 

podl=two2reg(x); 

load  pod2 

pod2=two2reg(x); 

load  pod3 

pod3=two2reg(x); 

load  pod4 

pod4=two2reg(x), 

load  pod5 

pod5=two2reg(x); 

load  pod6 

pod6=two2reg(x) ; 

load  temp 

temp=x; 

load  imul 

imul=x; 

load  imu2 

imu2=x; 

load  imu3 

imu3=x; 

remdfc 

********************************************************************** 
function  y  =  two2reg(x) 

%function  two2reg  takes  the  matlab  floating  point  number 
%  from  a  tattletale  BASS  current  meter  16  bit  two's 
%  complement  number  and  converts  it  to  signed  number. 

I=find(x>=2A15); 

y=x; 

y(I)=x(I)-2A16; 

********************************************************************** 

%  REMDFC  is  a  script  file  that  finds  the  time  gaps  from  the 
%  tim,3  and  tim,4  files  and  removes  the  corresponding  records  in 
%  all  the  variables.  The  file  version  assumes  that  all 
%  fluid  velocities  are  in  podl,  pod2,  and  pod3. 

%  This  version  will  not  work  for  sample  rates  slower  than  0.25  Hz. 

%  Last  revision  6/12/94 
% 

difs=diff(tim(:,3)); 


155 


difn=diff(tim(:,4)); 

I=find((difn~=l&difn~=-255)|difs>5|(difs>-55&difs<0)); 

[j,k]=size(tim);  %  to  remove  last  point 

I=[l;2;i;i+lj]; 

pod  1(1,  :)=[]; 

pod2(I,  :)=[]; 

pod3(I,  :)=[]; 

pod4(I,  :)=[]; 

pod5(I,  :)=[]; 

pod6(I,  :)=[]; 

temp(I,  :)=[]; 

imu  1(1,  :)=[]; 

imu2(I,  :)=[]; 

imu3  (I,  :)=[]; 

tim(I,  :)=(]; 

clear  difs  difn 

********************************************************************** 


function  imu3n  =  imu3swtch(imu3) 

%  function  imu3n  =  imu3swtch(imu3) 

%  switches  columns  2  and  3  of  imu3  to  correct  for  wiring  error 
%  after  this  function  imu3  will  be  in  x,y,z  like  the  other  imu's 

imu3  n=zeros(size(imu3 )); 
imu3  n(: ,  1  )=imu3  ( : ,  1 ); 
imu3n(:,2)-imu3(:,3); 
imu3n(:,3)=imu3(:,2); 

********************************************************************** 


%  n2cal.m  is  a  script  file  to  convert  vorticity  float  Bass  A/D  output  to 
%  physical  units  using  linear  calibrations  found  earlier. 

%  imulc  in  m/secA2 
%  imu2c  in  rad/sec 

%  imu3c  in  mG  and  assume  imu3,2  and  3  have  already  been  switched 

%establish  gain  and  zero  constants, 

%  these  calibration  constants  are  from  3/8/94 
imulg=[-0. 00030053  0.00030142  0.00030121]; 
imulz=[9.9501  -9.9478  -0.2539]; 
imu2g=[2.6848e-5  2.6965e-5  2.6691e-5]; 
imu2z=[-0.9006  -0.8951  -0.8767]; 
imu3g=[0.0 19056  0.018742  -0.017240]; 
imu3z=[-643.2  -581.4  533.3]; 

%  do  conversions 

imu  1  c=imu  1 .  *  (ones(size(imu  1  ( : ,  1  )))*imu  1  g)+ones(size(imu  1  (: ,  1 )))  *imu  1  z; 
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imu2c  imu2.*(ones(size(imu2(:,l)))*imu2g)+ones(size(imu2(:,l)))*imu2z; 
imu3c=imu3n.  *(ones(size(imu3n(:,  l)))*imu3g)+ones(size(imu3n(:,  l)))*imu3z; 
%optional  lines  to  clear  up  memory  space 
clear  imul  imu2  imu3  imu3n 

********************************************************************** 


function  [t]  =  n2t(n) 

%  function  t  =  n2t(bassint) 

%  This  function  computes  the  temperature  in  degrees  Centigrade 
%  from  the  Bass  recorded  integer  using  the  Steinhart  and  Hart 
%  equation. 

%  1/T(in  deg  K)  =  a  +  b*ln(r)  +  c*(ln(r))A3 
%  Coefficients  from  the  YSI  catalog  for  44030  thermistor 
a=0. 00 14051; 
b=0. 0002369; 
c=1.019e-7; 

% 

%calculate  thermistor  resistance  from  bass  recorded  A/D  number 
%using  bridge  and  amplifier  gains  from  vorticity  meter 
% 

[i,j]=size(n); 

r=(.0 1 6032*(2A  1 5  *ones(i,j)-n)+ 1 685  *ones(i  j))./(.3  87*ones(i  j)+  . . . 
5.83 14e-6*(n-2A15*ones(i,j))); 

% 

%compute  temperature  from  thermistor  resistance 
% 

invK=a*ones(i,j)+(b*ones(i,j)).*log(r)+(c*ones(i,j)).*((log(r)).A3); 
t=ones(i,j)./invK  -  273.15; 


%  svzer.m 

%  svzer.m  is  a  script  file  to  load  velocity  zeros  (in  signed  integers) 

%  for  the  double  sensor  deployment  south  of  Martha's  Vineyard  on  5/27/94. 
%  Sensor  1  into  cards  1,2  and  3  and  sensor  2  into  cards  4,5,  and  6. 

%  C,  1  diode  fixed. 

%  Sensor  1  zeros  from  svzer462-470  and  sensor  2  zeros  from  svzer333-450. 
%  last  revision  6/14/94 
% 

podlz=[5.8  15.9  -6.6  65.2]; 
pod2z=[-30.3  -17.7  -12.4  29.5]; 
pod3z=[124.3  -53.3  -39.5  -30.7]; 
pod4z=[23.9  -154.5  -29.3  256.8]; 
pod5z=[-33.1  -191.8  -109.1  214.6]; 
pod6z=[28.6  -173.5  34.2  265.5]; 
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%  n2v.m 

%  n2v.m  is  a  script  file  that  converts  45  cm  vorticity  meter  signed 
%  integer  velocity  numbers  and  calculates  corresponding  path  velocities 
%  in  cm/sec,  "circulation"  in  cm/sec,  and  vorticity  in  secA-l. 

%  This  file  assumes  that  appropriate  zeros(podlz,  pod2z,  etc.)  have 
%been  loaded  into  matlab. 

%  This  is  the  version  for  Matlab  4.0  and  two  sensors 
%  last  revision  6/13/94 
% 

avel=.  003 7 *  (pod  1  -ones(size(pod  1  (: ,  1 )))  *pod  1  z); 
acirc=avel( : ,  1  )+avel( : ,  2)+avel( : ,  3  )+avel( :  ,4) ; 
avort=acirc/47.34; 

bvel=. 003 7 *  (pod2-ones(size(pod2(  : ,  1 )))  *  pod2z) ; 
bcirc=bvel( : ,  1  )+bvel( :  ,2)+bvel( :  ,3  )+bvel( :  ,4); 
bvort=bcirc/47.34; 

cvel= . 003  7 * (pod3  -ones(size(pod3  (:,l)))*pod3z); 
ccirc=cvel( : ,  1  )+cvel( : ,  2)+cvel( : ,  3  )+cvel( :  ,4); 
cvort=ccirc/47 . 3 4 ; 

dvel= .0037*  (pod4-ones(size(pod4( : ,  1 )))  *  pod4z); 
dcirc=dvel( : ,  1  )+dvel( : ,  2)+dvel( : ,  3  )+dvel( :  ,4); 
dvort=dcirc/47.34; 

evel=.0037*(pod5-ones(size(pod5(:,l)))*pod5z); 
ecirc=evel( : ,  1  )+evel( : ,  2)+evel( : ,  3  )+evel( :  ,4); 
evort=ecirc/47. 34; 

fvel=. 003  7 *(pod6-ones(size(pod6( : ,  1 )))  *pod6z); 
fcirc=fvel( : ,  1  )+fVel( : ,  2)+fvel( : ,  3  )+fVel( : ,  4); 
fvort=fcirc/47 . 3 4 ; 

%  optional  line  to  clear  up  some  memory 
clear  podl  pod2  pod3  pod4  pod5  pod6 

********************************************************************** 

function  [watvelr]=veltr45f(avel,bvel,cvel) 

%  function  [watvelr]  =  veltr45f(avel,bvel,cvel) 

%  This  function  does  a  coordinate  transformation  for  measured  vorticity 
%  meter  velocities  from  the  sensor  paths  to  buoy  axis.  Since  there  are 
%  12  measured  velocities,  some  averaging  is  done. 

%  This  version  is  for  the  45  cm  path  vorticity  meter  after  the  chip  was 
%  changed,  ie,  for  the  Massbay  deployment  and  later. 

%  set  up  constants  to  save  later  computation 

s3o2=.25*sqrt(3)/2; 

o2s3=.25/(2*sqrt(3)); 

os3=.25/sqrt(3); 

s2o3=.25*sqrt(2/3); 

%  set  up  the  transformation  matrix 
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VT=[.125  -.125  -.125  .125  -.25  0  .25  0  .125  .125  -.125  -.125 

s3o2  o2s3  -s3o2  -o2s3  0  -os3  0  os3  -s3o2  o2s3  s3o2  -o2s3 
0  s2o3  0  -s2o3  0  s2o3  0  -s2o3  0  s2o3  0  -s2o3]; 

[k,l]=size(avel); 

watvelr=zeros(k,3);  %helps  matlab  interpreter 

%  do  matrix  multiplication 
watvelr=(  VT*  [avel  bvel  cvel]')'; 

********************************************************************** 


function  [vortb,watvort]=vortran2(avort,bvort,cvort,imu2c) 

%  function  [vortb,watvort]=vorttran2(avort,bvort,cvort,imu2c) 

%  This  function  transforms  the  coordinate  system 
%  from  sensor  vorticity  axis  to  buoy  axis. 

% 

vorts-[avort  bvort  cvort]; 

T=[l/sqrt(2)  0  -l/sqrt(2)  %set  up  rotation  transformation 

-l/sqrt(6)  sqrt(2/3)  -l/sqrt(6)  %matrix 

l/sqrt(3)  l/sqrt(3)  l/sqrt(3)]; 
vortb=(T*vorts')'; 

%  highpass  filter  pitch  and  roll  rate  to  take  out  rate  gyro  drift 
[vb,va]=butter(2,. 3/60, 'high'); 
imu2cf(:,  l)=filtfilt(vb,va,imu2c(:,  1)); 
imu2cf(:,2)=filtfilt(vb,va,imu2c(:,2)); 
imu2cf(:,3)=imu2c(:,3); 

watvort=vortb+2*imu2cf;  %vorticity  is  twice  the  rotation  rate 

********************************************************************** 


% buoyangb.m 
%buoyangb.m 

%  This  is  currently  a  script  file  to  compute  the  buoy 
%  Euler  angles  with  a  nonlinear  algorithm  and  a  complementary  filter. 

%  The  rate  gyros  are  integrated  and  highpass  filtered 
%  and  the  slow  response  is  from  the  linear  accelerometers 
%  and  magnetometer,  which  is  then  lowpass  filtered  and 
%  added  to  the  fast  angle.  The  slow  heading  estimate 
%  uses  an  intermediate  horizontal  reference  frame  to  calculate  heading. 

%  This  angle  is  then  used  to  recalculate  the  fast  estimate  of  Euler  angle 
%  with  the  true  nonlinear  F  inverse  matrix. 

%  This  version  uses  a  cutoff  period  of  30  seconds. 

%  This  file  is  for  a  150  msec,  sample  period,  and  assumes  a  321  Euler  Angle. 

%  Revised  10/27/94 
% 

BANGLEF=tpindint(imu2c,.  1 5); 
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% 

%compute  filter  coefficients 
% 

load  fill 530 

BANGLEFFO,  1  )=filtfilt(Bfh,  Afh,B  ANGLEFO,  1 )); 
BANGLEFF(:,2)=filtfilt(Bfh,Afh,BANGLEF(:,2)); 
BANGLEFF(:,3)=filtfilt(Bfh,Afh,BANGLEF(:,3)); 

% 

%compute  slow  angle  of  pitch 
BANGLES(:,l)=atan2(imulc(:,2),imulc(:,3)); 

% 

%  lowpass  filter  pitch 

BANGLESFO,  l)=filtfilt(Bfl,Afl,BANGLES(:,  1)); 

% 

%Add  together  slow  response  and  fast  response 
BANGLE(:,1)=BANGLEFF(:,1  )+BANGLESF(:,  1 ); 

% 

%Compute  roll 

BANGLES(:,2)=-atan2(imulc(:,l).*cos(BANGLE(:,l)),imulc(:,3)); 

BANGLESF(:,2)=filtfilt(Bfl,Afl,BANGLES(:,2)); 

BANGLE(:,2)=BANGLEFF(:,2)+BANGLESF(:,2); 

%Compute  magnetic  vector  in  intermediate  reference  frame 
imu3d=b2int(imu3c,BANGLE(:,l),BANGLE(:,2)); 

% 

%Compute  slow  estimate  of  heading,  filter  and  add  to  fast  estimate 
BANGLES(:,3)=atan2(imu3d(:,l),imu3d(:,2)); 

%unwrap  slow  heading 
headslow=unwrap(B  ANGLES( :  ,3 )); 

%Lowpass  filter  and  add  to  fast  estimate 

BANGLESF(:,3)=filtfilt(Bfl,Afl,headslow); 

heading=BANGLEFF(:,3)+BANGLESF(:,3); 

%Rewrap  heading 

B  ANGLE( : ,  3  )=heading-(floor((heading+pi)/ (2 *pi)))  *  2 *  pi; 

%  Recalculate  the  fast  estimate  of  Euler  Angles  with  the  nonlinear  update 
eulerdot=finv(imu2c,BANGLE); 

BANGLEF=tpindint(eulerdot,.  15); 

[i,j]=size(BANGLEF); 

BANGLEF=BANGLEF+ones(i,  1  )*BANGLE(  1 , 

BANGLEFF(:,  1  )=filtfilt(Bfh,Afh,BANGLEF(:,  1)); 

BANGLEFFO, 2)=filtfilt(Bfh,Afh,BANGLEF(:, 2)); 
BANGLEFF(:,3)=filtfilt(Bfh,Afh,BANGLEF(:,3)); 

%  add  slow  and  fast  estimates  of  Euler  angles 
BANGLE(:,  1)=BANGLEFF(:,  1)+BANGLESF(:,  1); 
BANGLE(:,2)=BANGLEFF(:,2)+BANGLESF(:,2); 
imu3  d=b2int(imu3  c,B  ANGLE( : ,  1  ),B  ANGLE(  :,2)); 
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% 

%Recompute  slow  estimate  of  heading,  filter  and  add  to  fast  estimate 
BANGLES(:,3)=atan2(imu3d(:,l),imu3d(:,2)); 

%unwrap  slow  heading 
headslow=unwrap(BANGLES(:,3)); 

%Lowpass  filter  and  add  to  fast  estimate 

BANGLESF(:,3)=filtfilt(Bfl,Afl,headslow); 

heading=BANGLEFF(:,3)+BANGLESF(:,3); 

%Rewrap  heading 

BANGLE(:,3)=heading-(floor((heading+pi)/(2*pi)))*2*pi; 

%Optional  line  to  delete  intermediate  variables 
%clear  headslow  heading 

********************************************************************** 

function  [area]=tpindint(func,delx) 

%fimction  [area]  =  tpindint(func,delx) 

%  This  function  uses  the  trapezoidal  rule  to 
%  calculate  the  running  area  under  data. 

%  The  error  should  be  proportional  to  the 
%  deltax  squared  times  the  second  derivative 
%  of  the  function.  If  fiinc  is  a  matrix,  [area]  will  be 
%  a  matrix  with  the  integration  carried  for  each  column 
% 

[m,n]=size(func); 

crude=cumsum(func); 

%  for  vector  case 
ifm=l|n=l, 

area=(crude- .  5  *(func(  1 )  *  ones(size(func))+func))*delx; 

%  for  matrix  case 
else 

area=(crude-.  5  *(ones(m,  1  )*func(  1 ,  :)+fimc))*delx; 

end; 

********************************************************************** 


function  [imu3  d]=b2int(imu3  c,psi,thet) 

%  function  imu3d  =  b2int(imu3c,psi,thet) 

%  This  function  does  the  coordinate  transformation  from  buoy  axis 
%  to  an  intermediate  computational  reference  frame  with  a  vertical 
%  z  axis.  It  is  used  by  the  buoyang2  algorithm  to  improve  it's 
%  heading  estimate. 

%  Revised  6/24/94 

%  The  following  vectors  are  components  of  the  rotation  transformation 
%  matrix  A21 

%  1 1  2  3| 
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|4  5  6| 
|7  8  9| 


% 

% 

Al=cos(thet); 

A2=sin(psi).  *sin(thet); 

A3=cos(psi).  *sin(thet); 

A4=zeros(size(psi)); 

A5=cos(psi); 

A6=-sin(psi); 

A7=-sin(thet); 

A8=sin(psi).  *cos(thet); 

A9=cos(psi).  *cos(thet); 

%size  the  imu3d  matrix  to  save  compiler  time 
imu3  d=zeros(size(imu3  c)); 

imu3d(:,l)=Al.*imu3c(:,l)  +  A2.*imu3c(:,2)  +  A3.*imu3c(:,3); 
imu3d(:,2)=A4.*imu3c(:,l)  +  A5.*imu3c(:,2)  +  A6.*imu3c(:,3); 
imu3d(:,3)=A7.*imu3c(:,l)  +  A8.*imu3c(:,2)  +  A9.*imu3c(:,3); 

function  eulerdot=finv(imu2c, BANGLE) 

%  function  eulerdot  =  finv(imu2c,BANGLE) 

%  This  function  multiplies  buoy  oriented  angle  rates  by 
%  the  F  inverse  matrix  for  the  nonlinear  Euler  Angle 
%  update. 

%  revised  10/27/94 

%  make  vectors  for  less  typing 

psi=BANGLE(:,l); 
thet=B  ANGLE( :  ,2); 
phi=BANGLE(:,3); 

%  matrix  F  inverse  |  1  2  3  | 

%  I  4  5  6  | 

%  |  7  8  9  j 

fl=ones(size(psi)); 

f2=tan(thet).  *sin(psi); 

f3=tan(thet).  *cos(psi); 

f4=zeros(size(psi)); 

f5=cos(psi); 

f6=-sin(psi); 

f7=f4; 

f8=sin(psi)./cos(thet); 

f9=cos(psi)./cos(thet); 
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eulerdot=ones(size(imu2c)); 


eulerdot(:,l)=fl.*imu2c(:,l)+f2.*imu2c(:,2)+f3.*imu2c(:,3); 

euIerdot(:,2)=f4.*imu2c(:,l)+f5.*imu2c(:,2)+f6.*imu2c(:,3); 

eulerdot( : ,  3  )=f7 .  *  imu2c( : ,  1  )+f8 .  *  imu2c( : ,  2)+f9 .  *  imu2c( : ,  3 ); 
********************************************************************** 


function  [EVect]=b2e(bvect, BANGLE) 

%fiinction  EVect=b2e(bvect,BANGLE)  6/7/94 

%  This  function  does  the  coordinate  transformation  from  buoy 
%  axis  to  earth  axis 
% 

%make  vectors  for  less  typing 
psi=BANGLE(:,l); 
thet=BANGLE(:,2); 
phi=BANGLE(:,3); 

%The  following  vectors  are  components  of  the  rotation  transformation 
%matrix  A  1 1  2  3| 

%  |4  5  6| 

%  |7  8  9| 

A 1  =cos(thet).  *  cos(phi); 

A2=sin(psi).*sin(thet).*cos(phi)-cos(psi).*sin(phi); 

A3  =cos(psi) .  *  sin(thet) .  *cos(phi)+sin(psi) .  *  sin(phi); 

A4=cos(thet).*sin(phi); 

A5=sin(psi).*sin(thet).*sin(phi)+cos(psi).*cos(phi); 

A6=cos(psi).*sin(thet).*sin(phi)-sin(psi).*cos(phi); 

A7=-sin(thet); 

A8=cos(thet).  *sin(psi); 

A9=cos(thet).  *cos(psi); 

%size  the  EVect  matrix  to  save  the  compiler  time 
EVect=ones(size(bvect)); 

E  Vect( : ,  1  )= A 1 .  *  bvect( : ,  1 )+ A2 .  *bvect( : ,  2)+ A3 .  *bvect( : ,  3 ); 

EVect(:,2)=A4.  *bvect(:,  1  )+A5.  *bvect(:,2)+A6.  *bvect(:,3)i 

E  Vect( :  ,3  )= A7 .  *bvect( : ,  1  )+A8 .  *bvect( :  ,2)+ A9 .  *bvect( :  ,3); 
********************************************************************** 

%  rel2abs4.m 
%  rel2abs4.m 

%  This  is  a  script  file  that  computes  water  velocity  in  earth 
%  coordinates.  This  algorithm  adds  the  omega  cross  r  term  in 
%  the  buoy  reference  and  computes  and  adds  the  buoy  velocity 
%  in  earth  coordinates. 

%  This  version  is  for .  1 5  sec  sampling  and  buoy  velocity  is 
%  highpass  filtered  in  buoy  coordinates  and  in  earth  coords 
%  revised  2/27/95 
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rl=157;  %  radius  in  centimeters 

r2=318.7; 

g=9.8;  %  m/secA2 

[n,  m]=size(watvelr  1 ); 

%  Add  in  omega  cross  r  term 

watvell=watvelrl+[imu2c(:,2)*rl  -imu2c(:,l)*rl  zeros(n,l)]; 
watvel2=rwatvelr2+[imu2c(:,2)*r2  -imu2c(:,l)*r2  zeros(n,l)]; 

%  Compensate  buoy  acceleration  for  gravity 

baccel=imu  1  c+[g* sin(B ANGLE( : ,2))  -g*sin(BANGLE(:,l)).*cos(BANGLE(:,2)) ... 
-g*cos(BANGLE(:,2)).*cos(BANGLE(:,l))]; 

% 

%  compute  and  use  highpass  filter  coefficients 
[bfh,afh]=butter(2,. 3/30, 'high'); 

Buoyvel=tpindint(b2e(baccel, BANGLE),.  1 5);  %  Velocity  in  m/sec 

%  highpass  filter  a  second  time  now  in  earth  coords 
Buoyvelf=zero  s(n,  m) ; 

Buoyvelf(:,  1  )=filtfilt(bfh,afh,Buoyvel(:,  1 )); 
Buoyvelf(:,2)=filtfilt(bfh,afh,BuoyveI(:,2)); 

Buoyvelf(:  ,3  )=filtfilt(bfh,  afh,Buoyvel( : ,  3 )); 

VEL 1  =b2e(watvel  1  ,B ANGLE)+ 1 00*Buoyvelf;  % VEL  in  cm/sec 

VEL2=rb2e(watvel2,BANGLE)+100*Buoyvelf;  %watvel  in  cm/sec 

%  optional  line  to  clear  intermediate  variables 
%  clear  watvell  watvel2  baccel  Buoyvel  Buoyvelf  bfh  afh 
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APPENDIX  C. 

This  appendix  catalogs  the  rotating  arm  bias  measurements  made  at  the  David 
Taylor  Model  Basin,  shows  the  measured  bias  in  constant  flow  of  the  15-centimeter  path 
vorticity  sensor,  and  shows  the  measured  buoy  response  functions.  In  the  rotating  arm 
tests,  the  45  cm  path  vorticity  sensor  was  mounted  on  an  apparatus  which  moved  the 
sensor  in  a  circular  motion  keeping  all  axes  constant.  All  points  of  the  sensor  traversed  in 
the  same  circular  motion  relative  to  the  tow  tank  carriage,  as  the  carriage  traveled  down  a 
large  tank  of  still  water.  In  these  tests:  the  y  sensor  direction  was  parallel  to  the  tank 
length  and  carriage  motion,  the  z  direction  was  vertical,  and  the  x  direction  was  horizontal 
and  across  the  tank  width.  The  apparatus  traversed  the  sensor  in  a  circular  motion  who’s 
axis  was  the  x  direction. 

The  measured  vorticity  means  in  constant  flow  (Fig.  3-10)  corrected  for  center 
strut  wake  bias,  is  repeated  as  Fig.  C-l.  The  standard  error  shown  is  the  standard 
deviation  of  each  sample  time  series  divided  by  the  square  root  of  the  number  of  samples 
of  each  mean,  and  is  a  measure  of  how  much  the  sample  mean  is  expected  to  fluctuate. 

The  measured  vorticity  means  with  the  sensor  rotated  at  ten  second,  seven  second,  five 
second  and  3.8  second  periods,  are  shown  in  Figs.  C-2,  C-3,  C-4,  and  C-5. 

The  measured  uncorrected  biases  of  the  15-centimeter  path  vorticity  sensor  towed 
through  still  water  are  shown  in  Fig  C-6  with  the  axes  labeled. 

The  measured  buoy  response  functions  are  plotted  in  Fig.  C-l.  The  heave,  surge, 
pitch,  and  pitch  rate  response  functions  were  calculated  from  measured  cross  spectra 
between  a  velocity  measured  at  0.83  meters  depth  and  the  buoy  response  variable.  For 
each  cross  spectrum,  the  amplitude,  phase,  and  coherence  are  plotted.  For  frequencies 
where  the  coherence  is  below  0.7  the  amplitude  and  phase  may  not  be  meaningful.  The 
heave  response  function  is  with  respect  to  vertical  velocity  at  0.83  meters  depth.  The 
North  surge  response  is  with  respect  to  north  velocity.  The  east  surge  response  was  the 
same,  as  one  would  expect  and  is  not  shown.  The  pitch  and  pitch  rate  response  functions 
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are  with  respect  to  horizontal  velocity  in  the  pitch  plane.  The  roll  and  roll  rate  response 
functions  are  the  same  respectively,  as  one  would  expect  from  symmetry. 


Fig.  C-l.  Vorticity  means  in  constant  flow  using  the  center  stalk  bias  correction 
described  in  the  text. 
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ig.  C-2.  Vorticity  means  with  the  sensor  rotated  at  a  ten  second  period.  The  x- 
xis  vorticity  is  corrected  for  the  center  stalk  wake  bias  as  described  in  the  text. 


Fig.  C-3.  Vorticity  means  with  the  sensor  rotated  at  a  seven  second  period.  The  x- 
axis  vorticity  is  corrected  for  the  center  stalk  wake  bias  as  described  in  the  text. 
This  figure  is  a  repeat  of  Fig.  3-15. 
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Fig.  C-4.  Vorticity  means  with  the  sensor  rotated  at  a  five  second  period.  The  x- 
axis  vorticity  is  corrected  for  the  center  stalk  wake  bias  as  described  in  the  text. 


Fig.  C-5.  Vorticity  means  with  the  sensor  rotated  at  a  3.8  second  period.  The  x- 
axis  vorticity  is  corrected  for  the  center  stalk  wake  bias  as  described  in  the  text. 
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Fig.  C-6  Measured  biases  of  15-cm  path  vorticity  sensor  towed  through  still  water. 
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Fig.  C-7.  Buoy  response  functions  for  heave,  surge,  pitch,  and  pitch  rate  relative  to 
velocity  at  0.83  meters  depth. 


APPENDIX  D. 

This  appendix  rotates  the  data  from  the  upper-boundary-layer  shear  deployment 
into  down-wind  and  cross-wind  components,  and  compares  these  to  two  models  of 
boundary-layer  behavior.  The  wind  coordinate  measurements  are  compared  to  two 
models  of  boundary  layer  behavior,  the  unstratified,  turbulent,  wall  layer  and  the  model  by 
Santala  (Santala,  1991).  The  unstratified-wall-layer  shear  is  given  by  (D-l)  and  assumes 


du  _  u, 
dz  K z 


(D-l) 


no  stratification.  In  this  equation,  u  is  the  downwind  drift  velocity,  u .  is  the  friction 
velocity  which  is  the  square  root  of  the  windstress  over  the  water  density,  and  k  is  von 
Karman's  constant  usually  assumed  to  be  0.4. 


The  second  model  that  the  data  were  compared  to  was  developed  by  Santala  to 
describe  ocean-upper-boundary-layer  behavior  measured  in  the  winter  off  the  northern 
California  coast.  In  this  section,  data  was  rotated  into  the  wind  coordinates  shown  in  Fig. 
D-l  with  z  positive  upwards.  Recall  that  the  direction  of  vorticity  is  defined  to  be 
perpendicular  to  the  velocities  and  directions  over  which  the  velocities  vary,  that  the 
vorticity  is  derived  from.  Defining  the  cross-wind  direction  for  shear  opposite  from  the 


downwind  shear  and  vorticity 
z  positive  up 


crosswind  shear 


Fig.  D-l.  Windstress  oriented  coordinate  system  used  in  this  appendix. 
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cross-wind  direction  for  vorticity  at  first  seems  odd,  but  with  this  convention  positive 
vorticity  corresponds  with  positive  shear.  In  this  coordinate  system,  the  Santala  model  is 
given  by  table  D-l.  In  this  table  vc  is  crosswind  drift  current  and  vd  is  the  downwind  drift 


Crosswind  shear 


dv  1.41m 

c  * 


dz 


for 


>  0.2  *105 

2 

w 


Downwind  shear 


dz 


u 


kz*1.64*105 


0 

/  \ 

^-1.24*105 

2 

U 

* 

KZ 


0  <  <  1.2*1 05 
2 
U _ 


■r  1.24*105  <  <  2.88*105 

V  2 

u . 

2.88*105  <  ¥- 


Table  D-l.  Santala  upper-boundary-layer  shear  model  in  coordinates  defined  by 
Fig.  D-l  -  . 


current.  This  model  predicts  large  crosswind  shear  and  a  zero  downwind  shear  layer  over 
a  log  layer  with  a  transition  layer  in  between.  As  windstress  increases,  the  inter-layer 
transitions  deepen,  and  at  some  depths,  the  shear  can  be  reduced  due  to  much  more 
effective  mixing  while  the  shear  stress  increases.  In  this  model,  Santala  does  not  mention 
a  development  time. 

The  measurements  for  this  comparison  were  from  the  upper-boundary-layer 
deployment  described  in  the  first  half  of  the  applications  section.  The  shear  measuring 
buoy  was  deployed  in  Buzzards  Bay,  Massachusetts  during  a  period  of  increasing 
windstress.  Windstress  was  measured  with  a  sonic  anemometer  and  calculated  with  the 
inertial  dissipation  method  (Fairall  and  Larson,  1986).  No  evidence  of  Langmuir  cells  was 
seen  during  the  deployment.  As  the  windstress  increased,  the  Santala  model  transition 
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depths  deepened,  Fig.  D-2.  When  the  shallower  transition  depth  sinks  below  a  sensor,  the 
model  predicts  zero  downwind  shear. 


Fig.  D-2.  Santala  model  layer  transition  depths  predicted  from  the  measured 
windstress  during  the  Buzzards  Bay  deployment 

The  buoy  data  were  rotated  into  directions  downwind  and  crosswind,  and  the 
vorticity  predicted  shear  were  compared  to  an  unstratified  wall  layer  and  the  Santala 
model.  The  downwind  shear  from  vorticity  in  the  0.83  m  depth  sensor  is  shown  plotted 
with  both  the  Santala  model  and  wall  layer  model  in  Fig.  D-3.  In  the  first  half  of  the 
deployment,  the  shear  follows  closely  the  wall  model.  But,  as  the  windstress  continued  to 
increase,  the  shear  diverged  from  the  model  and  decreased  as  the  windstress  increased. 
With  enough  time,  the  measured  shear  could  have  approached  the  Santala  model. 

The  downwind  shear  from  vorticity  in  the  2.45  m  depth  sensor  is  shown  plotted 
with  both  models  in  Fig.  D-4.  The  measured  vorticity  was  much  larger  than  either  the 
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wall  model  or  Santala  model.  The  measured  stratification  over  the  lower  sensor  was 
larger  than  the  stratification  between  sensors,  Fig.  D-5.  The  top  thermistor  failed  leaving 
only  thermistors  at  the  bottom  of  the  shallow  measurement  volume,  and  at  the  top  and 


Fig.  D-5.  Stratification  in  lower  sensor  and  between  sensors 


bottom  of  the  deeper  measurement  volume.  The  stratification  over  the  shallow  sensor  was 
assumed  to  be  the  same  as  the  intersensor  stratification  for  the  Richardson  number  plot  in 
the  applications  section.  This  larger  stratification  can  explain  the  less  effective  vertical 
turbulent  mixing  implied  by  the  larger  measured  shear. 

Crosswind  shears  from  measured  vorticity  and  Santala  model  shear  are  plotted  for 
both  sensors  in  Fig.  D-6.  The  measured  crosswind  shear  in  both  sensors  was  in  the  same 
direction  as  predicted  by  Santala.  The  measured  crosswind  shear  in  the  shallow  sensor 
starts  at  about  zero  and  does  not  approach  the  Santala  model  until  near  the  end  of  the 
deployment  when  the  windstress  was  strong.  The  measured  crosswind  shear  over  the 
lower  sensor  was  larger  than  the  Santala  model,  which  can  again  be  explained  by  the 
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larger  stratification  over  the  lower  sensor.  The  wall  layer  model  predicts  no  crosswind 
shear. 

In  most  of  the  measurements,  the  wall  layer  model,  when  compensated  for 
stratification,  was  close  to  the  measured  results  for  low  windstress.  As  the  windstress 
grew,  the  measured  results  started  to  approach  the  Santala  model.  The  shears  measured 
during  this  deployment  may  have  been  influenced  by  the  shape  of  Buzzards  Bay  as  the 


Fig.  D-6.  Crosswind  shear  from  vorticity  and  windstress  for  both  sensors. 

wave  directional  spectrum  was.  The  shear  measuring  buoy  is  a  useful  tool  to  measure 
shear  in  the  upper  five  meters  of  the  ocean,  and  with  further  and  longer  deployments 
should  help  ocean  scientists  to  better  understand  this  i  mportant  part  of  the  ocean. 
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